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Summary. — In this paper we consider a simple model of an electron- 
photon cascade in which the transformation probabilities are functions 
of absorber depth. The mean and variance of the number of electrons 
in the cascade as a function of absorber depth are obtained and their 
properties diseussed. An equation which yields the probability that the 
cascade will eventually terminate is also given. 


4. — Introduction. 


In a previous paper (1) we have presented an introduction to the Bellman- 
Harris theory of age-dependent stochastic processes (?), and considered a simple 
model of a cascade process based on this theory. L. JANossy (*) has employed - 
an approach similar to that of Bellman and Harris in the derivation of the 
G-equations of cascade theory; and the same technique has been used by 
J. M. BLATT (4) in developing a stochastic theory of tr ‘acks in nuclear emul- 
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sions. For a discussion of the Bellman-Harris theory, and related theories, 
we refer to reference (5). 

In this paper we wish to consider a model for an electron-photon cascade 
in which the transformation probabilities are functions of absorber depth. 
In order to formulate this model it is necessary to consider a generalization 
of the Bellman-Harris theory, since these authors considered only constant 
transformation probabilities. This generalization has been obtained by W. A. 
O’N. WAUGH (5). The Waugh process can be described (in the nomenclature 
of cascade theory) as follows: Let X(t) denote the number of electrons at 
absorber depth t in a one-dimensional electron-photon cascade, and let G(t) 
denote the probability distribution of the absorber depth, say 7, travelled 
by an electron from its formation until it is transformed. Hence G(T), 
0 <t<oo, is the integral distribution for all paths of lengths less than or 
equal to 7. At the end of its path of travel the electron is transformed into 
n electrons (n= 0,1, 2,...) with probabilities q,(t), with each new electron 
having the same distribution G(r) for the depth of the absorber it will travel 
before being transformed. For example: q(t) is the probability of absorp- 
tion; gr) is the probability that the electron is absorbed and one new 
electron formed, or that the electron continues to penetrate the absorber; 
and q.(t) is the probability that two new electrons will be formed, the original 
(parent) electron being absorbed, etc. In each case q,(t) isa prescribed 


function of 7; and these probabilities satisfy the condition > q,(t)=1 for 


n=0 
all TETE co: 

Let X(t) be an integer-valued random variable representing the number 
of electrons at depth t; hence the stochastic process {X(t), t >0} represents 
the development of the cascade. Now let p(a, t) = Pr (X(t) =), & — 0, 1, 2, ...; 
and let 
(1) m(s, t) = Y pla, t)s®, |s|<1 

x=0 
denote the generating function of the probabilities p(x, #), starting with one 
incident electron at absorber depth zero. When the transformation proba- 
bilities depend on absorber depth, the generating function x (s, t) satisfies the 
functional equation 
¢ . 
(2) a(s,t) ={h[a(s,t— t), t]dG(t) + s[1 — G(t)],. 
0 
(6) A. T. BHARUCHA-REID: Elements of the Theory of Markov Processes and their 


Applications, (New York, forthcoming). 
(5) W. A. O’N. WavuGu: Biometrika, 42, 291 (1955). 
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where 


(3) h(s, t) = x Qn(T)8" , Is <1 


n=0 


that is, h(s, 7) is the generating function of the transformation probabilities 
In(T). 

Differentiation of (2) with respect to s yields integral equations of the 
renewal type for the moments of Y(t). For example, the mean or expected 
value of X(t) is given by the solution of the integral equation 


(4) EX(t) = mit) = fatom te vjdG(t) EL Gi), 
where : 

oh 
(5) k,(t) = ae 2 Lx 


and the second derivative of z(s,t), which we denote by m,(t), satisfies the 
integral equation 


(6) malt) = [Em — IF a(x) + Jisceimst— aac, 
where 
(7) = | 


Since the variance of X(t) is given by 
(8) Var X(t) = ms(t) + mi (0) —[m,(t)]", 


the solutions of (4) and (6) enable us to obtain the mean and variance of X(t). 
Higher moments can be obtained by continued differentiation. 


9. — An electron-photon cascade model. 


In this section we consider a simple model of an electron-photon cascade. 
In order to consider a particular process of the Waugh-type, it is necessary 
to specify the functions G(t) and h(s, t) which characterize the process. We 
now assume 


(9) G(t) = 1—exp[— À], A>0, 


TA 
Le) 
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and 
| go(t) = 1—f exp[—az], x >0,, dep 
(10) | amava], 0<y<1 
| Qu(t) = 1 — (qo(t) + ge(7)) = (6 — y) expl— ar]. 


The functions (9) and (10) which characterize the cascade can be described 
as follows: An electron at absorber depth # has probability ZAt+0‘At) of 
being transformed (experiencing a collision) in the interval (t, t+At). When 
transformation takes place the electron can: (i) be absorbed with probability 
q(t), Where q(t) is an increasing function of 7, equal to 1—f when 7 —0, 
and equal to unity when t = oo; (ii) produce a new electron and continue to 
penetrate the absorber, or be absorbed and produce two new electrons with 
probability q(t), where gs(7) is a decreasing function of 7, equal to y wien 
t=0, and equal to zero when 7= x; or (iii) continue to penetrate with 
probability q,(t), where q,(t) is also a decreasing function of 7, equal to 
(8 —y) when t= 0, and equal to zero when t= co. In order for q,(t) to be 
non-negative we must have f > y. 

From the above, we see that the possibilities when transformation takes 
place are that the number of electrons in the cascade can be reduced by one 
with probability q(t), unchanged with probability q,(t), or increased by one 
with probability q(t). In the definition of the transformation probabilities 
the variable t represents the depth of absorber penetrated by an electron from 
its formation until it is transformed. We also remark that the negative-expo- 
nential form of G(t) expresses the Markov property of the stochastic cascade. 

Before proceeding we add a few remarks concerning the relationship of 
the above model to some of the stochastic models considered by other workers. 
Firstly, we note that for parameter values &— 0, 6 =y—=1, the above model 
is identical with the pure birth model considered by W. H. Furry (7). Secondly, 
for «= 0, and f #y, we obtain a model of the birth-and-death type (1, 8). 
Finally, our model is a generalization of a model considered by N. ARLEY (8) 
in which the probability of absorption in the interval (t, At) was assumed 
to be ut A+ o(At). For the Arley model, and our model, the total cross-section 
is a function of +. 

While energy considerations are not explicitly taken into account in our 
model, the exponential form of the transformation probabilities provides a 
simple approximation to more realistic cascade processes if one assumes that 


(€) W. H. Furry: Phys. Rev., 52, 569 (1957). 


(8) N. ARLEY: On the Theory of Stochastic Processes and Their Applications to the 
Theory of Cosmic Radiation (New York, 1949). 
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the energy dissipated by a particle is proportional to the thickness of absorber 
traversed since its formation. 


From (2), (9) and (10) we obtain the following integral equation for the 
generating function 2(s, t): 


t 


(11) as, t) =|{ galt) + nas, t — 7) + galt) x*(s, tm}: 
‘Aexp[— At]dt + s exp[— At] = 
= fe — Bexp[— at] +[(8 — y) exp [— ronan t—t) + 


+ y exp [— at]z*(s, t— 1)}4 exp[— Ar]dr +s exp[— Aue 


The above integral equation can be reduced to the differential equation 


(12) = yln*(s, t) —[« +A — (B — y)) |a(s, t) + 


+ œ + A(1— p) + a(s —1) exp [— AM]. 


A solution to this non-linear and non-homogeneous equation which admits 
a series expansion in s is difficult to obtain. Hence we concern ourselves 
with the study of the mean and variance of the number of electrons in the 
cascade. 


3. - Integral equations for the moments. 


From (4), (9) and (10), we see that the mean number of electrons in the 
cascade is given by the solution of the integral equation 


(13) m,(t) = ap + y)m,(t — t) exp [— (a + A)t] dt + exp [— At]. 


This integral equation can be reduced to the differential equation 


dm, (t) I 


+ —{a +1 (B+y)]}m() +e exp[- A], 


(14) 


20 
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the solution of which is 


(15) malt) = Die exp [— At] — 


—H(B + y) exp [—[a + (1 — (6 + y))]t}}- 


1 
Lx — A(B + 


In solving (14) we have used the initial condition m,(0)= 1, è.e. the cascade 
is initiated by one electron. Should the cascade be initiated by N electrons, 
which are statistically independent and produce independent cascades, the 
mean number of electrons at absorber depth ¢ is obtained by multiplying (15) 
Dy aN. 

In connection with the mean number of electrons in the caseade, it is of 
interest to consider two problems. The first is to determine the mean number 
of electrons for large absorber depth, i.e. as t > co. The second is to find 
the value of ¢ which maximizes the mean, and then obtain the associaed mean 
number of electrons for various parameter values. For the first problem it 
is easy to show that as t + © the following limiting values, depending on 
certain relations between the parameters, obtain: 


| 0, for f+<1 
(16) Hesse à 2, tor Gex, Paypal, 


od, Or eA fp 


For the second problem, it is equally easy to show that the value of t 
which maximizes m,(t) is | 


iin Sth atop bre AR 
(17) HE pa] cs nr] - lE 


Since the above expression for #,,, can assume a large number of values, 
depending on the various parameter values, it is of interest to calculate pa 
for the special cases considered in the determination of mf. For B+y—1 
we see that ¢,,.=0; hence the maximum mean value is one (or N in the 


case X(0)= N), and the curve of the mean number of electrons decreases 
to zero as t +00. For f=y 


E (2ACa + A(1 — 28)] 
(18) lune = rag 18 | : i 


In order for ¢,,, to be positive and finite we must have 


2pLa + A(1 — 2B)] 


x 


>1 and a—248>0; 


3 
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hence we must have À <æ(2B)-1. In this case m,(¢) attains its maximum at 
tax given by (18), and then decreases to zero. For a= À, B=y=1, we find 


that #,, = co; hence m,(t) does not achieve a maximum for finite absorber 
depths. For t= co, we have from (16) m= 2. Finally, for «<A, B=y=1 


1 {2a — 4) 
(19) Lu = ——, 1 . 
a— 24 ae aes 
Since we assume x <A, t,,, for this case does not exist. 
In general, if we put ¢=¢,,, in (15) we obtain 


A 


| S x [CB +r) Le + 20 — (8 + y)) oes) 
M (tmax) [a — Ap + y)] | is | 14 
ARMOR re RAR pe py rer 
fe — A8 +1 à 


as the mean number of electrons at the maximum absorber depth t,..: 
We now consider the variance of the number of electrons in the cascade. 
From (6), (9) and (10) we have 


(21) m,(t) = 2ya fm —t)P exp[— (a + A)r]dr + 


+ A(B + [mat —t)exp[— (« + A)t]dt. 


0 
This integral equation can be reduced to the differential equation 


(22) POL as fort A(1 — (B + y)) ]ma(t) + 2yALm,(t)}? . 


Since m,(t) is a known function, the above equation can be solved to give m,(t). 
The solution of (22) is 


(23) malt) = 2yAB-* exp [At] {a*(2B + A)" exp [(B + A)t] — 
— 2a(B—«)A-1 exp [— At] + (B— a)2A-! exp [At]} + Cexp[At], 


where A——[x+À(1—(8+7))], B=a—A(B+y), and © is the constant of 
integration. Using the initial condition m,(0) — 0 we find 


CAB AREA 2a(B — a)" + (B— x) AT}. 
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From (8), (15) and (23), we obtain the rather complicated expression for the 
variance of the number of electrons in the cascade: 


(24) Var X(t) = B-? {(2yA(2B + A)-1— 1]o? exp [— 22] + B exp [— At] — 


— 2a(2y +1)(B — a) exp[(2A + B)t] + (2yZ24-1— 1)(B— a)? exp [2At] + 
+ [B(B—«) — 2yA[a?(2B + A)-1-+ 20(B—a«)s-4 + (B — «)?A-]] exp [At]} . 


If we consider the asymptotic variance for large absorber depths, we find 


0, for f+y<1 
(25) lim Var X(t) = i pa 


| co, for x<A, P=y=1. 


For the case x=4, and B=y=1, A=0 and B<0; hence in order to 
obtain the variance in this case it is necessary to determine m,(t) and m,(t) 
when the above relations between the parameters obtain. From equation (22) 
we obtain 


(26) amt) _ atm, 
where 
(27) mat) = 2— exp [— At}. 


Solving (26) with the initial condition m,(0) — 0, we obtain 


(28) malt) = 8(At + exp [— At}) — exp [— 24t] — 7. 
Hence 
(29) lim Var X(t)=.c6, 


COT =A pie 


4. — Coefficient of variation or relative fluctuation. 


We now add a few remarks on the coefficient of variation (or relative 
fluctuation) of the number of electrons in the cascade. The importance of 
this function has been stressed by several authors (8). The coefficient of va- 
riation is defined as : 


(Var X(t))# 


(30) V(t) = EX (i) 
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From (8) we see that 


Mo(t) il 


(31) ROSE ES 


Because the expression for m.(t) is rather complicated, the function V(t) is 
difficult to discuss for finite absorber depths. However the asymptotic coef- 
ficient of variation V*= lim V(t) can be obtained from (16), (25) and (29). 
From a study of certain special cases we observe that the introduction of 
transformation probabilities depending on absorber depth increases the flue- 
tuation. 


5. — Probability of cascade termination. 


Since we have introduced a probability of absorption, q(t), it is of interest 
to determine the probability, say P(t), that the cascade will terminate at 
absorber depth #, t< 0, i.e. P(t)=Pr (X(t)=0). It has been shown (6) that 
P(t) is a solution of the integral equation, 


(32) Pit) = [fan + gle) PU 7) + walt) Pt 1} AG (2) = 


= fi pente +8 — rexel PO— 9) + 
+ y exp[— ar] P2(t — 7)}A exp [— it|dr. 


Like equations (11), (13) and (21), the above integral equation can be reduced 
to à differential equation and solved to give P(t). However, in order to deter- 
mine the probability that the cascade will eventually terminate, we have only 
to consider an algebric equation. 

From the theory of discrete branching processes (°) it is known that the 
probability of eventual termination is given by the smallest non-negative 
solution of the functional equation n(s)=s. For age-dependent processes it 
has been shown (2) that this probability, say Q — lim P(t), is given by the 
unique non-negative root, less than unity, of h(Q)= Q. For the Waugh pro- 
cess it is necessary to consider, in place of the above equation, the functional 


equation 


(33) h(Q, tT) = Q. 
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Hence, for our particular model, we seek the solution of 


Cha h(Q, t) =A4+BQ+00=0, 
where 
2 fui mar = fa1— Bexp[— arl] exp[— Àt]dr = seem 
co Le] - 
8 | 1(T)g(T) At fe» ) exp[— (x + A)t] dr = Oi 
J 2)t]d 2 
cfa dr =|yZexp[— (x + À)r] Ae 


0 


with A+B+C=1. Therefore, the probability that the cascade will even- 
tually terminate is given by the smallest non-negative solution, less than 
unity, of 

(35) yaQ? —[a —A(lL— (B+y))10 + x +2A(1— 8) = 0. 

For values such that 


[a — A(1— (B+y))] > 4yA(x + A1 — 8) 


two real roots, say è, and &,, of (35) exist; and Q — min (¢,, ¢.). 


RIASS-ULN TO=(*) 


Nel presente lavoro si considera un modello semplice di una cascata di elettroni 
e fotoni in cui le probabilità di trasformazione sono funzioni dello spessore dell’assor- 
bitore. Si ricavano la media e la varianza del numero degli elettroni nella cascata in 
funzione dello spessore dell’assorbitore e se ne discutono le proprietà. Si dà anche 
un'equazione che fornisce la probabilità che la cascata si arresti. 


(*) Traduzione a cura della Redazione. 
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Analytic Properties of Scattering Amplitudes 
as Functions of Momentum Transfer. 


H. LEHMANN (*) 


Institute for Advanced Study Princeton - New Jersey 


(ricevuto il 31 Luglio 1958) 


Summary. — Scattering amplitudes are shown to have analytic pro- 
perties as functions of momentum transfer. The partial wave expansions 
which define physical scattering amplitudes continue to converge for 
complex values of the scattering angle, and define uniquely the ampli- 
tudes appearing in the unphysical region of non-forward dispersion rela- 
tions. The expansions converge for all values of momentum transfer 
for which dispersion relations have been proved. 


4. — Introduction. 


The purpose of this note is to derive some properties of scattering ampli- 
tudes which follow from causality in relativistic quantum theory. It will be 
shown that a scattering amplitude has—for fixed energy—analytic properties 
as a function of scattering angle or momentum transfer. This consequence 
of causality is distinct from the existence of dispersion relations (4) which 
express analytic properties of a scattering amplitude as a function of energy 
for fixed momentum transfer. However, our results are of interest mainly in 
connection with dispersion relations for non-forward scattering. They imply 
that for all values of momentum transfer for which dispersion relations have 


(*) Now at the University of Hamburg, Germany. 

(1) M. L. GOLDBERGER: Phys. Rev., 99, 979 (1955). 

(2) K. SYMANZIK: Phys. Rev., 100, 743 (1957). 

(3) N. BOGOLIUBOV, B. Mepvepev and M. PoLivanov: lecture notes. Translated 
at the Institute for Advanced Study, (Princeton, 1957). 

(4) H. J. BREMERMANN, R. Orume and J. G. TAYLOR: Phys. Rev., 109, 2178 (1958). 
These papers contain numerous other references. 
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been established (?)—in x-N scattering this is the case for momentum transfer 
< (2:(2m—+u)/(2m — u))} Au —the so-called non-physical region is COR PIE 
determined by the physical phase shifts via the partial wave expansion. This 
possibility has been suggested previously (°). 

We disregard spins and consider the scattering of a charged particle of 
mass m (nucleon) with initial momentum p, final momentum p’, and a neutral 
particle of mass u < 2m (meson) with initial momentum k, final momentum ki, 
We assume in addition that these particles are not coupled to other particles 
of charge zero and mass < 2 or of charge one and mass < mu. Our results 
are valid also for the scattering of equal particles (m = u) under corresponding 
restrictions on the mass spectrum. 

Ag is well known, several equivalent expressions for the scattering matrix 
in terms of Heisenberg operators may be given. With the notation 


<p'k' out|pk in) = <p'k' in| pk in) + (2x) d(p +k — p'—k')T 


the amplitude 7 may be written as (5) 


(1) Ti — [are exp Le 7 <p'|R'A (5)4 (- 3) Ip> 5 
or 
(2) i = — fare exp ela {0|R'A (5) v(— JP ind : 


A(«) and y(x) are the meson and nucleon field operators. R’ denotes a retarded 
commutator. For example, 


R' A(x) p(y) = — i(le — p?)(Cy — m?) 6(a@ — y)[A(x), ply)] . 


The state vectors refer to incoming or outgoing particles with definite momenta 
as indicated. 

Eq. (1) is used in the derivation of dispersion relations. Eq. (2) on the 
other hand yields directly information about the scattering amplitude as a 
function of momentum transfer, since—in the center of mass system—this 
variable appears only in the exponential. However, more information is ob- 
tained in both cases by observing that from either (1) or (2) the following 


(5) M. L. GOLDBERGER: Proceedings of the Sixth Annual Rochester Conference, 
(New York, 1956). 


(5) H. LEHMANN, K. SYMANZIK and W. ZIMMERMANN, Nuovo Cimento, 6, 319 (1957). 
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expression for the imaginary part of the amplitude may be derived, 


(3) ImT= af adr dr, exp È = a, i(k phas| 


| 4 di) x 2 72 
sora) (Step +m rie (rs 


> denotes a sum over all states with total four-momentum p+k. Hq. (3) is 
“4 


(on the energy shell) equivalent to the unitary requirement for the scattering 
amplitude. It was first used by BoGoLIUBOV et al. (*). A simple proof of this 
relation is given in Sect. 3. 

The general method used to obtain explicit consequences of these expres- 
sions may be described as follows: 
In each case the scattering amplitude—or its imaginary part—appears as the 
Fourier transform of a retarded commutator or of a sum over products of 
such commutators. Therefore it is simply related to the Fourier-transform 
of the corresponding unretarded commutator which is in the case of (2) or (3) 
given by 

x 


4) ig) = [are exp tigeicol[i (5) f(—5)]IP +47 


where j(x) = (D —y?)A(x); f(a) = (A —m?)y(x). We know about F(q): 


(4a) it is the Fourier transform of a function that 
vanishes for space-like # ; 


F(q) vanishes unless 


te HED, 
DE ko E go 0 and as +f i > mi 
(4b) or 
k 2 
Porta g,>0 and (PE — gh > mt. 


The latter statement follows directly if a sum over intermediate states is 
introduced in (4). ™ and m, are the masses of the lowest intermediate states 
which contribute to the two terms of the commutator. In the m-N case 


M, = Ils Ma = mM: 
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Dyson (7) has solved the problem of finding a representation of all functions 
satisfying these conditions. His result is: For F(q) to satisfy (4a) and (4b) 
it is necessary and sufficient that it can be represented as 


(5) F(a) = | du | drte( qo — o)O[(g — te)? — x2} plu, 22). 


g(u, x?) is arbitrary if the vectors (p+k)/2+u and (p+k)/2—wu both lie in 
the forward light-cone and 


x >Max Lo: M — V (e -_ DE Mz — VE — o} È 


g vanishes outside this region. It depends, of course, also on the quantum 
numbers y and on p+k. All our results will be based on applications of 
Dyson’s theorem. 

For the Fourier transform of the retarded commutator which appears in (2) 
or (3) we have the relation (q'= (q,, q)) 


dqoF 
F,(g) = SEO ; Img>0, 
0 


if F, is sufficiently bounded. 
Inserting (5) gives 


dx? (u, x? 
6 ue . 
(6) F,(g) = [a uf 


In general we cannot expect F, to be bounded enough for (6) to hold in this 
form. The necessary modification (*) does not alter the analytic properties 
we are interested in. It is therefore sufficient to discuss (6). We shall use 
it both for Eq. (2) and Eq. (3). 


2. — Momentum transfer properties of scattering amplitudes. 


By inserting (6) into Eq. (2) we obtain 


(7) alc (liu diga DR) 
~ Da ((k'— p')/2— wu)? — x2’ 


(7) F. J. Dyson: Phys. Rev., 110, 1460 (1958). 
(§) R. Jost and H. LEHMANN: Nuovo Cimento, 5, 1598 (1957), Eq. (4-5). 
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y is an invariant function of the vectors u, p, k. The integration extends over 


the region given in (5). We choose the center of mass system to evaluate (7) 
and introduce the variables 


i) 
wate, se, 
(8) or 
pw Um + WW — (MA, ogg g = 1 — 24 
K? = 4W? , SPE fe 


g depends then only on u?, %, uk, 2, W. It vanishes outside 


[ 0 <u <W/2; — W/2 Lu Lu <W/2—4; 
ls, | x>Max {0; m,—V(W/2 + Wo) — U5 m,— VW[2— uo)? — wt. 


Introducing polar co-ordinates in u-space, (7) becomes : 


(10) T(W, cos Ÿ) = = hy af [affa 


P(Uo, U?, COS € sin B, x?, W) 


R°+ w+ ne (vo + (m? — p?)/2W)? — cos Sn 
2Ku sin f 
or 
© 27 

(a, cos «, W) 

(11) T(W, cos d) — [sofa] Ci (3 "ni a)’ 
aw) 0 

with 


1 
p(w, cosa, W) = gie fan fan far for 


5 Ket ut + 28 — (uo + (m° — p°) 
È “i 2Ku sin p 


2 M. Elo Us cos a sin B, x, W). 


The lower limit x,(W) is determined by 


æo(W) = 


[K2 + ut + xt — (Wot (m° ee 
in ? 
| 2Ku 
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if u,, u, x vary over the region (9). The minimum can be calculated in an ele- 


mentary manner. The result is 


(mi p?)(m;— me) |} 
K*[W? — (m, — m,)*] 


(12) æ(W) = |1 + 


In (11) the scattering angle cos Ÿ appears only in the denominator. We may 
therefore consider cos 3 as a complex variable and the scattering amplitude 
as an analytic function of cos Ÿ. . Moreover, this can be done separately for 
the real and imaginary parts of the amplitude. 

Singularities of these functions can occur only if the denominator on the 
right hand side of (11) vanishes. That is for 


cos 9 = x-cosa + iva? — 1 sine. 
We have therefore the following result: 


(13) The real part and the imaginary part of the scattering 
amplitude are analytic functions of cos Ÿ, regular inside 
an ellipse in the cos Ÿ-plane with center at the origin 
and with axes x, V#?—1. 


We shall see presently—making use of Eq. (3)—that the imaginary part of 
the amplitude is regular in a larger domain, namely: 


(14) Im T(W, cos d) is regular in cos Ÿ inside an ellipse with 
center at the origin and with axes 2% — 1; 2% /, xa — 1 
æ% is given by (12). 


Using cos Ÿ =1— (2A?/K?) we can, of course, re-express (13) and (14) as analytic 
properties of the scattering amplitude as a funetion of momentum transfer. 

These results (we defer the proof of (14)) lead—using well-known mathe- 
matical theorems (*)—to the following properties of the partial wave expansion 
of the scattering amplitude. Let 


T(W, cos 8) = az, i (21 + 1)0,(W) P,(cos 8) , 


(15) n K 
with C(W),=—- DW d cos UE cos ee (COs d) . 


Er 


(*) E. T. Wurrraker and G. N. Warson. A course of modern analysis, 4th Ed., 
{Cambridge 1940), p. 322; G. Szeaô, Orthogonal Polynomials, (New York, 1939), p. 238. 
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The Legendre series converges inside the domain of regularity of the re- 
presented functions; i.e. for cos inside the ellipses (13) or (14) for Re T or 
Im T respectively. Also 


ni 1 
13a lim | Re €,(W) |!" < — = 
( ) lo 1 Lo ci Va — 1’ 
(14a) lim |Im 0 (W)|™" < 


(am) + Vari — D 
Taking into account the unitarity relation 
(16) Im 0,(W) > [Re C,(W)? + [im CAW)F, 


we may note that (13a) is actually a consequence of (14a); i.e. if Im T(W, cos #) 
is regular in the domain (14) it follows immediately that Re T(W, cos d) is 
regular in (13). 

We cannot conclude, of course, that the amplitude T(W, cos?) actually 
has singularities on the boundary of the domains (13) or (14). Using more 
physical information it may well be possible to improve these results. 

To discuss the connection of the above statements with the non-physical 
region of dispersion relations, let us consider Im T(K?, A?), the imaginary 
part of the amplitude, as a function of c.m. momentum and momentum 
transfer—the physical region is given by K? > A2. However, in the dispersion 
relation Im T(K?, 4?) is needed for all K? > 0. The expansion 


LAVE 2 A2 
(50) Im, A) = TE (2 + 1) Im OAW) Br (1-58): 


defines a continuation of Im 7(K?, A?) into the non-physical region. The 
series converges if 


At < Kei = K+, 


i.e. it converges for all K?>0 provided 
A? < Min {Ka} . 


This leads to the restriction 


- Il Nuovo Cimento. 
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in the x-N case. (A? < 2u? for equal particle scattering with M, = M = 2p). 
We have still to show that the continuation given by (15a) is indeed the 
correct definition of the non-physical region in the dispersion relation. 


3. — Connection with dispersion relations. 


To identify Im T(K?, A?) as given by (15a) with the dispersion relation 
integrand it is convenient to consider (1) (which defines a function 7 for 
arbitrary real vectors k, k’) not only on the energy shell (k? = k? = pu?) but 
for the more general case ka=k°=t; keeping p? =p = m?. T can then 
be considered as a function of 


o TEV +P) palpa pa PBI dos È, 


2V(p + p')° = 


We derive first Eq. (3). 
(1) leads directly to 


(AT) Sime TIw,e, 217) 3 fato exp i is = Si x 


ali (5).i(-3)]l® = 


il 
=5{M(0, €, 4°) — M(- 0, 6, 4%}, 


with 
eh ate abe ie EN INR ei 
18) ato, &, 4") =[arwexp |i ("> *) a] ali (5) (5) = 


= (2x) 2 <p'|j(0) |p + k, y> <p +k, y|9(0)|p> - 


Let pin(p') denote the annihilation operator for an incoming nucleon with 
momentum p’. Then 


<p'|9(0) |p +k, y> = <O| pin (p')j(0) |p +k, y> = <O|[p in (p'), j(0)]|p +k, y> 


if (p+k—p'P=k?=C < 9? 


since <0|j(0)gin(p')|p+%,y)=0 in this case. With the relation (5) 


PR Nr a : 
in JON = pri [Ate exp Lip’ IR’ Ayla) , 
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and an analogous treatment of the second factor in (18) we have 

(19) HM (0°C, 43) = 2a fato, dt, exp [ip'x, — ip] > <O| R'A(0)y(2.)|p + k, y>- 
"i 

nali Seale) 

RA (7) yt (-5) Lo». 


We note that the imaginary part of the physical scattering amplitude is given by 


<p +k, y|R’A(O)y*(a,)|0> = 22 | du, dur, exp 


> <0 | Ri A (3) (2)? +k, y> <p + ky y 


(20) Im T =} M(o, u?, A?). 


The term M(—a, u?, A?) does not contribute since it vanishes for © > V/A? +. 
Hence (19) is on the energy shell equivalent to (3). 


To find analytic properties of M as given by (19) we use again the integral 
representation (6). This leads to 


st dxid'u,Adx3D(U,,%1) Ur, %2, P + k) 
21 A? 
(21) M(o, b; L= ie 12 — un)? — x1][((k —p) 12 — u)? — |? 
where 


Du, #15 Un %2, PYk) = > Plu, #1) p+k) p, (ur, #2, p+k) 


is built up from the weight functions y, corresponding to the individual terms 
on the right hand side of (19). © is a real, invariant function which satisfies 
the support conditions (6) in each pair of variables u, x separately. We choose 
the center of mass system and replace w by 


W2= 20V A? + m? + 242 + m + È. 


Then 
02) Mo, £4) ==: 
d'u, dx} dtu, dx P(r, ui, Xi, Uzo » Ua x, Uy Uy/U,U2, W) 
[ (m? — c)/2W tire U10)? — x — (k'— u;)? | [((m? = C)/2W + Uso) ?— x (k —u,)?| 


With polar co-ordinates 


(23) M(W, 6, A? ) a KAD rg mom ifn ww 


Du o, UÈ, 22, cos « sin f sin B, + COS f, COS Ba, W) 
[a(£) — cos (9 — x)|[*2(2) — cos (x — ot) | i 
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K+(6) + ui + mi — ((m? —C)/2W + Uio)®. 


sr È 2K(t)u, sin f; E 
W2 (pra Cet 4 2W2 
rio = OMR, 


We note that 


27 


1 . 1 fi 
Ju [æ, — cos (0 — x) | [Lx — cos (ye x) | 


(a |(Wae? — 1)) + (æ (Var — 1) 


#7 : 
ta (2,0, + Va? — Vai — 1 — cos (8 — a) 


From (23) analytic properties in € and A? follow. On the energy shell, 
i.e. î= u?, we can introduce y = 7,7,+Vx?—1V23—1 as a new integration 
variable. Only integrations over y and « remain; the other integrations result 
only in a new weight function ®(y, cos a, W). 

The minimum value of y is 


Yo = Min {a,9, + Vel =1V ai —1} = 203 — 1. 


Therefore 


27 = 


M(W?, A) =2Im 7 = [av fas Dy, cos a, W) 


y — cos (I — a) 
2a°—1 0 


This proves the statement (14). 

It can be seen now that the analytic continuation of Im T defined by 
(15a) yields the correct non-physical part of the dispersion relation integrand. 
In the proofs of these relations (*) it is shown first that—as a consequence 
of Eq. (1)—a dispersion relation in © holds if È is taken real and €< — A?. 
The absorptive part in this relation is M(W, ¢, 4?) as given by (18) and (23). 
The dispersion relation for the physical value ¢=y? is then obtained by 
analytic continuation in ¢, provided M(W, È, A?) is an analytic function of È 
regular for Re € < u? in a neighborhood of the real axis. It follows from (23) 
that this condition is satisfied if 


At < Mina: 


This is also the condition for the convergence of the Legendre series. The 
absorptive part of the dispersion relation is then given by (23) with ¢€ = p?; 
i.e. the non-physical region is obtained by analytic continuation in A? which 
can be carried out by the Legendre expression (15a). 
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The possibility of evaluating the non-physical region in this manner has 
been discussed earlier (5). While no proofs were given, it was believed on the 
basis of threshold arguments that such a procedure could be valid only if 
A? < u?, due to a branch point of the scattering amplitude as a function of A?. 
We have shown that the expansion converges also for higher values of A?; 
the limit being A? = 2u? in the case of equal particle scattering. We believe 
that this is due to the fact that the real and imaginary part of the amplitude 
are separately analytic functions of A? and have different properties. For 
the dispersion relation only the imaginary part is needed and it has a larger 
domain of regularity. The mentioned branch point is likely to be present in 
the real part. 

While we have no good reason to believe that our results are best possible, 
the expected appearance of a singularity in the real part gives us—via the 
unitarity relation (16)—an upper limit to the values of A? for which the Le- 
gendre expansion for the imaginary part might converge. In the case of equal 
particle scattering the expected branch point of the real part at 4?=W?/4 
leads to the limitation A? < 8u? for the Legendre expansion of the imaginary 
part. 


* * * 


The author would like to thank Professor J. R. OPPENHEIMER for the 
hospitality of the Institute for Advanced Study, and to acknowledge financial 
support from the National Science Foundation. He is indebted to several 
physicists at this Institute for valuable discussions. 


RIASSUNTO (*) 


Si dimostra che le ampiezze di scattering hanno proprietà analitiche come funzioni 
del trasferimento dei momenti. Gli sviluppi parziali delle funzioni d’onda che definiscono 
le ampiezze fisiche di scattering continuano a convergere per valori complessi dell’an- 
golo di scattering e definiscono unicamente le ampiezze che compaiono nella regione 
non fisica delle relazioni di dispersione non in avanti. Gli sviluppi convergono per tutti 
i valori del trasferimento dei momenti per cui sono state dimostrate esatte le relazioni 


di dispersione. 


(*) Traduzione a cura della Redazione. 
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Summary. — The photographic plate method has been used to detect 
3He particles emitted from helium gas exposed to bremsstrahlung of 
maximum energy 70 MeV. The *He(y, n)*He cross-section is found to 
vary between values of (0.63-+0.10) mb at 40 MeV and (0.25+0.04) mb 
at 60 MeV, the integrated cross-section being 7.5 MeV mb between these 
limits. The angular distribution of the *He particles in the C.M. system 
has a peak in the forward direction and may be described by the ex- 
pression f(0) = B sin?0 (1+ C'cos 0) with © approximately equal to 0.5. 


1. — Introduction. 


The photodisintegration of ‘He nuclei by bremsstrahlung of high energy 
was observed in helium-filled cloud chambers by GAERTTNER and YEATER (1) 
and by NICOLAI and GOLDWASSER (?). Their results indicated that the inte- 
grated cross-sections for the (Y, n) process are rather greater than the cor- 
responding values for the (y,p) process, but detailed excitation functions 
were not obtained. The photoneutrons from helium were detected by FER- 
GUSON, HALPERN, NATHANS and YERGIN (*), who showed that the (y, n) cross- 


(*) E. R. GAERTINER and M. L. YEATER: Phys. Rev., 88, 146 (1951). 

(*) V. O. NicoLar and E. L. GOLDWASSER: Phys. Rev., 94, 755A (1954). 

(3) G. A. Ferguson, J. HALPERN, R. NATHANS and P. F. YERGIN: Phys. Rev., 
95, 776 (1954). 
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section increases from zero at the threshold energy of 20.6 MeV to a value 
of about 1.3 millibarns near 25 MeV, higher energies being unattainable. Ex- 
periments by FULLER (*) with photographic plates arranged to record photo- 
protons emitted from helium gas yielded a (y, p) cross-section curve rising to 
a maximum value of 1.8 millibarns near 26 MeV and falling sharply to less 
than 0.1 millibarn at 40 MeV. A curve of similar shape for the same reaction 
was obtained by REID, SWINBANK and ATKINSON (°) using a cloud chamber. 
DE SAUSSURE and OSBORNE (*) used photoplates to record the *He particles 
emitted in the (y n) reaction and found a (Y,n) cross-section falling slowly 
from 0.5 millibarn at 40 MeV to 0.1 millibarn at 120 MeV. Despite the fact 
that the (y, p) and (y, n) processes are essentially similar in ‘He there appeared 
to be wide discrepancies between the two cross-sections, especially near 40 MeV. 

More recently a cloud chamber investigation of the photodisintegration 
reactions in ‘He has been reported by GORBUNOV and SPIRIDONOV (°) in a paper 
giving detailed results for the (y, p) reaction. Before this work became known, 
the present authors had begun to scan photoplates which had been exposed 
in helium irradiated by bremsstrahlung from the 70 MeV beta-synchrotron at 
Queen’s University, Kingston, Ontario. The preliminary results gave infor- 
mation about the (y,n) cross-section above 40 MeV which was reported to 
the Conference on Photonuclear Reactions at Washington, D.C., in April 1958. 
More detailed studies have been pursued in order to compare the data with 
the (y, p) cross-section curve published by GORBUNOV and SPIRIDONOV. Par- 
ticular attention has been paid to the (y,n) cross-section between 40 MeV 
and 60 MeV and to the angular distribution of the *He particles in the centre- 
of-mass co-ordinate system. 


2. — Experimental method. 


In general the method closely followed that used in experiments on the 
photoprotons emitted from oxygen, carbon and nitrogen (8). The narrow 
bremsstrahlung beam entered a high-pressure helium chamber through a lead 
collimator 35 cm thick and the charged particles emitted from the gas were 
detected in photoplates disposed symmetrically about the beam. In preli- 
minary experiments it was found that most of the particles recorded were 


(4) E. G. FULLER: Phys. Rev., 96, 1306 (1954). 
J. M. Rep, P. SWINBANK ahd J. R. ATKINSON: Physica, 22, 1142A (1956). 
J. M. RED: private communication (1958). 
(5) G. DE SAUSSURE and L. S. OsBorNE: Phys. Kev., 99, 843 (1955). 
(7) A. N. Gorsunov and V.M. SPIRIDONOV: Zu. Eksper. Teor. Fiz., 38, 21 (1957). 
(8) D. L. Livesey: Can. Journ. Phys., 34, 1022 (1956); 35, 987 (1957). 
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protons produced in one or other of the reactions 


‘He + hv = 2H + *H —19.8 MeV, 
‘He + hy = 2H + ™+*H — 26.0 MeV, 


4He + hy = 21H + 21n — 28.2 MeV. 


There were however several dense tracks possibly due to *He produced in 
the (y, n) reaction 


4He + hy = in + He — 20.6 MeV. 


One of the plates used (hereafter called Plate A) was of Ilford type E-1, 
and grain-counts were carried out on all tracks longer than 30 um in a limited 
area of this plate. Most of the tracks were shown to be protons by compa- 
rison with proton tracks produced within the emulsion by the fast neutron 
background. A second group consisted of short tracks as dense as those of 
radioactive «-particles, and there was also an intermediate group ascribed to 
deuterons and tritons. Many of the supposed *He tracks were too short and 
too dense to be suitable for grain-counting so a different method had to be 
used for their identification. 

In subsequent experiments the photoplates were reversed in their holders 
immediately after exposure and a weak source of Actinium A, B and O was 
placed at the centre of the chamber which was immediately evacuated. This 
further exposure (of 30 minutes) introduced «-particles of ranges up to 30 um 
into the plates from the «wrong » direction. The «-particle tracks occurred 
in most microscope fields and greatly facilitated identification of the *He tracks. 
Since the ionization density in a short *He track is very nearly equal to that 
produced by ‘He, a criterion was adopted whereby a supposed *He track was 
accepted only if it was at least as dense as a nearby «-particle track. As a 
result a few genuine *He tracks may have been rejected, but elimination of 
triton tracks, numerous in these plates, should have been ensured. 

A check on the reliability of the method was provided by variations in the 
processing of the plates. In the series of Ilford C-2 plates (B, C, D and E) the 
development was successively reduced, so that while plate B received over- 
night immersion in cold HQ developer (*) plus 30 minutes in HQ at 20 °C, 
plate E received the cold immersion only. These plates showed considerable 
variation in the density of proton tracks, but investigation of similar areas 
revealed that the proportion of ‘He tracks accepted was (12 -+1)% of all 


() W. K. Dawson and D. L. Livesey: Can. Journ. Phys., 34, 241 (1956). 
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tracks exceeding 10 um in length within the limits of experimental error. This 
ratio did not depend on the local abundance of Actinium a-particles or on the 
area selected. 

A possible source of error lay in the occurrence of recoiling ‘He nuclei 
ejected from the gas by fast neutrons. This was checked experimentally by 
measurements of the fast neutron flux, both within the beam and in the region 
of the plates, and by counting the number of helium tracks entering plate A 
in the « wrong » direction. The results showed that less than 5% of the re- 
corded tracks could be due to this process. Radioactive contamination conld 
also cause background etfects, but no definite a-particles could be detected 
in an exposure without gas in the chamber. A survey of possible sources of 
error indicated a maximum systematic error of 10%, but this is in the oppo- 
site sense to the errors caused by rejecting too many sHe tracks. 

In the final examination of the plates the two observers worked indepen- 
dently, scanning considerable areas at the edges nearest the beam, and also 
at the far edges to check the expected variation of yield with distance from 
the beam. The photon energy and angle of emission in the center-of-mass 
co-ordinates were calculated for each *He track accepted. The 3He energy 
varies rapidly with laboratory angle for a given photon energy; the lowest 
sHe energy recorded (4.5 MeV) corresponded to a photon of 34 MeV at the 
smallest useful angle of 30° and to a photon of 45 MeV at the maximum angle 
of 150°. Accordingly a correction was applied for the loss of tracks at back- 
ward angles for photon energies below 45 MeV. An additional correction for 
the events occurring outside the useful range of angles amounted to not more 
than 5% above 45 MeV. 

The integrated (y, n) cross-sections for various ranges of photon energy 
were calculated from the expression: 


n = (Far?) N(y/R?)(*/27) 5 


where n = number of events recorded in a certain energy range per unit 
area of plate, 


F = effective photon flux per unit area of the beam and per unit 
energy range, 
y = radius of the collimated beam, 
= number of helium nuclei per unit volume of gas, 


y = perpendicular distance of the beam axis from the plane of the 
photoplate, 


R = mean distance from the beam axis to the scanned area 


and x = integrated cross-section over the chosen energy range. 
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Corrections were worked out for the effects of beam polarization, finite 
beam cross-section and finite scanned area, but in no case did these exceed 2 Ge 
The value of F was calculated from the X-ray dosage in nominal roentgens 
indicated by the synchrotron monitor, which was calibrated in a separate 
experiment on the activation of copper foils by the reaction ®Cu(y, n)°Cu. 
The integrated cross-section for this reaction was taken to be 0.60 MeV barns 
from the threshold up to 70 MeV (1°). This figure led to a direct estimate of 
the photon flux at 18 MeV in the bremsstrahlung spectrum. At all other photon 
energies the flux F was calculated by normalizing a Schiff curve to the point 
fixed at 18 MeV. It should be pointed out that the maximum synchrotron 
energy was not stabilised and the procedure may not be reliable above 60 MeV; 
at lower energies the values of F should be accurate to + 10%. 


3. — Results and discussion. 


The detailed results obtained for the integrated cross-sections in different 
energy regions are shown ih Table I, where the standard errors are based 
on statistical fluctuations only. The different plates and observers gave quite 
consistent results and the final mean values were weighted according to the 
areas of plate scanned. The total integrated cross-section from 40 MeV to 
60 MeV amounts to 7.5 MeV mb with a standard error of + 0.9 MeV mb 


TABLE I. 
Nn: Ob- ATEN Mes Integrated cross-sections (MeV mb) 
ME SS em |(40+45)MeV(*)| (45-250) MeV | (50-60) MeV 
A DLL 1.70 SRO 2.4 +0.5 1.9 +0.5 2.3 +0.6 
B DLL 1.50 dro 2.6 +0.5 3.1 +0.6 2.3 +0.6 
B DLL 1.50 4.4 3.4 +0.8 2.85 +0.75 2.3 +0.8 
C IGM 4.00 990 2.2 °+0.5 1.75 40.4 2.4 +0.5 
D DLL 2.00 303 2.9 +0.6 3.2 +0.7 3.2+0.8 
E DLL 4.00 3.55 2.8 +0.45 2.25+0.45 2.8+0.6 
E IGM 2.80 4.4 2.2 +0.7 1.5 +0.5 3.2 +0.9 
Weighted means: 2.57 +0.25 2.21 +0.25 21-03 
Total cross-section: 7.5+0.5 MeV mb 
(*) Corrected for losses at angles exceeding 110°. 


(9) A. I. Berman and K. L. Brown: Phys. Rev., 96, 83 (1954). 
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after full allowance has been made for all errors arising in the experiment. 
In comparison, NICOLAI and GOLDWASSER found a value of 28 MeV mb for 
the entire energy range up to 135 MeV, 
and GoRBUNOV and SPIRIDONOV obtained 
about 6.8 MeV mb for the (y,p) reaction 
between 40 MeV and 60 MeV. 

The shape of the cross-section curve was 
found by plotting the data at 2.5 MeV 
energy intervals, as shown in Fig. 1. The 
results were everywhere consistent with a 
cross-section falling smoothly from an esti- 
mate of (0.63+0.10) mb at 40 MeV to 
(0.25 0.04) mb at 60 MeV, with a best value 
of (0.374-0.04) mb at 50 MeV. This curve. 
agrees closely with the (y, p) data of GOR- 
BUNOY and SPIRIDONOV but falls rather more 
rapidly than the earlier (y,n) curve of 
DE SAUSSURE and OSBORNE. The new results 
should reduce slightly the bremsstrahlung- 
weighted cross-section as estimated by RUSTGI Salta lobbed “against: photon 
and LEVINGER (1), but their theoretical con- energy. Upper: estimated (+, n) 
clusions and those of FOLDY (13) are not cross-section as a function of 
materially affected. photon energy; the smooth curve 

In order to obtain the final angular di- has been fitted to pire ane best 
stribution of Fig. 2, the distribution of *He da a ee Re 
tracks in centre-of-mass co-ordinates was 
corrected for the variation of the effective 
solid angle of collection with angle in the laboratory system. The difference 
between the centre-of-mass angle and laboratory angle is very nearly con- 
stant over a wide energy range, so the correction was applied by dividins 
each ordinate by the sine of the corresponding mean laboratory angle. The 
distribution for all events above 45 MeV shows a distinct forward peak 
near 70° and the data above 50 MeV are also asymmetrical about 90°. Attempts 
were made to represent the distributions by functions of the form {(0) = 
— A+B sin?0 (1+0 così). In all cases the value of A/B was not significantly 
different from zero. The best fit was obtained with 0 approximately equal 
to 0.5, but experimental errors at the extreme angles precluded a real test 
of the function f(0). The earlier work of DESAUSSURE and OSBORNE gave al- 


di 


u01)235 SS04I PajOWIISZ 


Number of* events x 2.5 MeV interval 


06 50 60 7 
Photon energy MeV 


Fig. 1. — Lower: number of (+, n) 


(11) M. L. RUSTGI and J. S. LEvINGER: Phys. Rev., 106, 530 (1957). 
(2) L. L. FOLDY: Phys. Rev., 107, 1303 (1957). 
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$0 most a pure sin?@ distribution for the 
energy range 36 MeV to 84 MeV. The 


results of GORBUNOV and SPIRIDONOV 


NO) on the (y, p) reaction are also in agree- 
: pe il ment with a very low value A/B; 
their value of C is 1.05 + 0.16 for all 
a events between 30 MeV and 170 MeV. 
x A test of the effects of the cosine 
* term is provided by calculations of the 
3 ratio of events at»forward angles to 
S hv50 MeV 


those at backward angles. In our work 
this is 14 +0.2 for energies above 
45 MeV, whereas GORBUNOY and SPIRI- 
DONOV obtained 2.4 for photoprotons 
in the same energy region. The two 
cases are not strictly comparable because 
3 : the (y, n) reaction involves the emis- 
ere sion of only one charged particle. The 
nates, after correction for the variation È ; ioe 
of solid angle of collection with angle forward predominance in *He emission 
in the laboratory system. Lower: all implies backward emission of neutrons, 
events above 50 MeV. Upper: all events such as occurs in the photodisintegration 
above 45 MeV, the curve representing of the deuteron at high energies. Apart 
Ba erro oh ER Ad, from differences in the values of C, 
the (y,n) and (y,p) reactions appear 
to have similar properties in the energy region investigated. 


Relative 
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Fig. 2. — Angular distribution of *He 
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RIASSUNTO (*) 


Il metodo delle lastre fotografiche è stato usato per rivelare particelle di He emesse 
da elio gassoso esposto a bremsstrahlung di energia massima 70 MeV. La sezione d’urto 
‘He(y, n)*He si trova variare tra valori di (0,63 +0,10) mb a 40 MeV e (0,25 +0,04) mb 
a 60 MeV, la sezione d’urto integrata essendo 7.5 MeV mb tra questi limiti. La distri- 
buzione angolare delle particelle di *He nel sistema del centro di massa ha un picco in 
direzione anteriore e si può descrivere con l’espressione 1(0) = B sin?6 (1+0C cos 6) 
con C approssimativamente uguale a 0.5. 


(*) Traduzione a cura della Redazione. 
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On the Bound State Problem in Quantum Field Theory (*). 


W. ZIMMERMANN (*) 


Institute for Advanced Study - Princeton, N. J. 
Physics Department, University of California - Berkeley, Cal. 


(ricevuto il 1° Agosto 1958) 


Summary. — A causal and invariant scalar field involving a stable bound 
state is investigated. A formula for the S-matrix is derived and it is 
shown that the bound state can be described by a local and invariant 
field operator. For simplicity only the case of spin zero particles and 
bound states is considered; however, the extension to other cases is 
possible. 


Introduction. 


For a relativistic quantum mechanical system with an energy momentum 
operator P, the one particle states are defined as eigen states of — Pi with 
a discrete non-vanishing rest mass (*°). In general, there may be several kinds 
of particles, each characterized by a certain value of rest mass, spin,charge, etc. 
Usually a division is made between elementary and composite particles. But 
it seems to be hard to define this distinction in a convincing manner. In the 
conventional formulation of quantum field theory each of the elementary 


(*) This research was supported in part by the United States Air Force under 
contract no. AF 49(638)-327 monitored by the AF-Office of Scientific Research of the 
Air Research and Development Command. 

(+) On leave of absence from the Max-Planck Institut für Physik, Gottingen, 
Germany. | 

(1) For a discussion of the particle aspect in quantum field theory compare R. HAAG: 


Dan. Mat. Fys. Medd., 29, 12 (1955), especially chap. no: 
(2) We exclude the case of particles with zero rest mass. 
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particles is described by a basic field operator whereas the composite particles 
(like the deuteron in the ground state, etc.) appear as the stable bound states. 
of the system. But this definition of a composite particle depends, of course, 
on the formalism. The same particle which is regarded as an elementary par- 
ticle in one formalism may be a composite particle in another (**). 

In this paper we want to investigate whether the principle of microscopic 
causality sheds any new light on this question. We consider a model of a 
causal and invariant scalar field A(x) describing just two kinds of particles, 
an elementary particle of mass m and a composite particle of mass M, both 
of spin zero. It will be shown that it is possible to define a field operator B(x) 
for the composite particle (explicitly expressed in terms of the original field 
A(x)) which is in some respect analogous to the field operator of an elementary 
particle. The invariant operator B(x) satisfies the requirement of microscopic 
causality, furthermore the S-matrix can be expressed by the vacuum t-fune- 
tions (5) of the field operator A(x) and B(x) exactly in the same way as in 
the case of elementary particles only. Therefore the principle of microscopic 
causality offers no possibility of distinguishing between elementary and com- 
posite particles. With respect to the S-matrix the elementary particles as 
well as the composite particles of the model are described by invariant field 
operators which satisfy the requirement of microscopic causality (8). 

The formalism developed in this paper is closely related to other recent 
investigations. Starting from similar requirements one can derive the same 
results by applying the method of strong operator convergence in the form 
developed by HAAG (*°). An equivalent formalism was also obtained from 
the recursion formulae for retarded functions recently proposed by NISHIJIMA 
as an axiomatic formulation of quantum field theory (1911). 


(?) Compare E. FERMI and C. N. YANG: Phys. Rev., 76, 1739 (1949). In this paper 
the hypothesis that z-mesons may be composite particles is discussed. 

(*) In this connection an interesting example is Giirsey’s model of a theory of 
elementary particles (F. GiRsEY: to be published). There, following a suggestion of 
Heisenberg, all particles appear as composite and no elementary particles correspond 
to the basic fields. i 

(5) The vacuum t-functions are defined as the vacuum expectation values of mul- 
tiple, time ordered, operator products. 

(5) This possibility was first mentioned by N. BocoLJuBOVv: unpublished lecture notes. 

(°) The S-matrix of the model considered is causal according to the definition given 
in H. LEHMANN, K. Symanzik and W. ZIMMERMANN: Nuovo Cimento, 6, 319 (1957). 

(*) R. HAAG: Proc. of the Lille Conference on Mathematical problems of the quantum 
theory of fields (1957), in print. 

(9) R. HaaG: preprint, to be published. 

(1°) K. NISHIJSIMA: Progr. Theor. Phys., 17, 765 (1957). 

(1!) K. NISHIJIMA: preprint, to be published. 


3960 


ON THE BOUND STATE PROBLEM IN QUANTUM FIELD THEORY 599 


1. — General conditions. 


We consider the model of a neutral scalar field described by a hermitian 
operator A(x) assumed to be invariant under the inhomogeneous Lorentz group. 
We assume that the principle of microscopic causality 


(1) [A(æ), A(y)] = 0 for (e — y} > 0 


holds and that no negative eigenvalues appear in the energy and rest mass 
spectrum. The operators A(#) will be supposed to form an irreducible ope- 
rator ring. For simplicity we assume that there are just two non-vanishing 
discrete eigenvalues m? and M? of — P, and that 


(Q, A(a)®) #0 if PD m? D 
but 
(Q, A(x)P) = 0 


if — a en) EN UM 
(Q, A(x) A(y)P) #0 
In addition we suppose that the states ® and Y have spin zero. 

This situation may occur for example in the case of A‘-coupling if there 
is an elementary particle of mass m and a stable two particle bound state 
of the mass M. 

In order to describe the bound states we introduce the operator (!) 


(2) B(w, €) = TA(w + €)A(w — È) 
and define (1%) incoming fields, according to WIGHTMAN (1), by 
AL) = Ale) +[Analms een, 


(3) 
| B,,(a, €) = Bla, €) + | Anal M, 0a — 2')j(@", È) da’, 


(2) £ may be taken as spacelike or timelike 4-vector, but it is supposed that & #0 
(3) The integral expressions in eq. (3) are to be understood in the sense of weak 
operator convergence which means that 


(®, Ain (a) ¥) = (D, A(x)P) + [ane (m, ca a')(®, j(a')P) da’, 
between any normalizable state vectors ® and W. 


(4) The following definition of the asymptotic fields A,,(x), Aou(x) was first sug- 
gested by WIGHTMAN: private communication. 
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with the current operators 
ja) = (EH; —m?)A(x), 


02 
(3) O: SES 
j(æ, €) = (D, — M?)B(e, &), Cm, 


The outgoing fields A,,,(v), Bu (1, £) are correspondingly defined with the 
‘help of the advanced functions A,,,(m,#) and 4,,,(M, x). The general con- 
ditions under which the so defined operators Aow(x) exist have been given and 
investigated in details by GREENBERG and can easily be extended to the case 
.of the operators Bout CSST EN 

Finally we demand that the incoming field operators A,, (2), B,,(#, £) to- 
gether form an irreducible operator ring and correspondingly the outgoing 
fields A, (0), Bou (4, €). 

So far we have listed the general requirements which we need in the work 
which follows. We conclude this section by deriving some simple properties 
.of the operators Au (x), Bout (x, È). AS a consequence of definition (3) they are 
solutions of the Klein-Gordon equations for the masses m and M, respectively: 


(O, — m?) Au (x) =0, 


(4) | 
(CO, — M°)Bou (a, E) = 0. 


We have the invariance properties 


0A out (x) 


. OBout x, 
2 = —iPy; Aow(a)], quela, 6) 


OL, 


(5) Us i[P,, Bout (a, &)] ’ 


de, 


Agu (La) = U(L)Apw(#)U(L)-1, Bow (La, LE) = U(L)Bou(t, &)U(L)", 


in 


for an arbitrary Lorentz transformation L. (U(L) denotes the unitary operator 
transforming A(x) into A(Læ).) These invariance properties are easy to prove 
if (3) is written in momentum space. 

As a consequence of (4), (5) the vacuum expectation values of the incoming 
fields vanish: 
(6) (Q, A out (0) Q) =) , (Q, Bout (2, È) Q) = 0. 


For example, 
02 
M?(Q, B,, (x, 8) A 25 (2, B,(a, &)Q) Un 
u 
because P= 0. 


(5) O. W. GREENBERG and A. &. WIGHIMAN: preprint, to be published. 


(*) In Sect. 2 we will use somewhat more restrictive conditions than Greenberg, 
for details see reference (2°), 
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Furthermore Acut(x) and Bout (x, €) satisfy the asymptotic conditions 


lim (®, A,(t)W) = (D, Aow, W), 


(7) mar 
Rr (D, Fr (t, È) P} = (®, Bout, (È) d} ’ 
with 
| 0 0 
Bag == Fe M à x a ®lm) __ * (y) == Bila | 
rit, È) if a be È) gi (0) — F*(x) Bie Bia, DE 
(7°) CRETE 
Bontp (€) = — ifase {rate 15) ce F*(e) Ge F*(x) te Bout (x, al, 


(correspondingly the definition of A, and Aout) for any normalizable solution 
f(z), F(x) of the Klein-Gordon equation 


(D — m?) f (v7) = 0, 


O M*) F(a) = 0. 


Eq. (7) can be proved by forming the integral expression (7’) for both sides 
of (3) and taking the limit t + + oo. 


2. — Commutation relations for the asymptotic fields. 


21. Elementary particles. — Our first aim is to derive commutation re- 
lations for the incoming (or outgoing) fields defined in Sect. 1. We begin with 
the case of elementary particles and prove the relation (1?) 


(8) [A,,(@), Ain (y)] = [Am (0), Aone (¥)] = À A(m, & — y) 


if A(#) is normalized in the following manner. Let ®, be a one-particle state 
with energy momentum eigenvalue k, 


PD = k, Dr ’ 
(Dr, Dr) a 2ko ds (k — k') ’ 


—K =m, kr=|Vk+m|. 


(7) This was already shown by GREENBERG, reference (15) under the assumption 
that the equal time commutator [4(2) A (y) 21, 18 à e-number. Here we use a different 
method which can be extended to the case of bound states. 


E 39 - II Nuovo Cimento. 
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Then it follows from translation invariance that 


(Q, A(a)®,) = € CPE 


(2x)? 


with the constant 


¢ = (2x) (9, A(0) Di) . 


Now A(z) shall be normalized by the condition c = 1. 
It is easy to determine the matrix elements of A,, (x) between the vacuum 
state and an arbitrary state vector. From the definition (3) it follows that 


exp [tka] 


for the one-particle state ®,. On the other hand if ® is an eigenstate of 
— Pi with a rest mass x? 4 m? we have 


(Q, A,, (2) D) =0 
because 


(m? — x2)(9, A,, (x) D) a 


follows from (4) and (5). 
With this result we can calculate the vacuum expectation value of the 
commutator (8). We see that 


(CA (AO) =] $e (2, Au (D) (Pas A (O) = à Am, 2 — 0) 


(9) ($2, [Aout (ar) Aout (y)] Q) = à A(m, x — y). 


Now we turn to the operator form (8) and using the conditions of Sect. 1 
we shall show that 


(i) the commutators of the incoming and the outgoing fields coincide: 


(10) [A out (0) REA y)] = - [ A: in A,, (y)] 
and that 


(ii) the commutators of the incoming field is a e-number: 


(11) LA, (0); An (y)] = (9, [An (0), 4, (y)]Q) . 
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Relation (8) is, of course, a consequence of the statements (i), (ii) and Eq. (9). 
In order to prove (i) and (ii) it is convenient to expand the field operators 
with respect to a complete orthonormal system {f,(«)} of the positive frequency 
solutions of (LJ — m?) f(x) = 0: 
A(t) =D f(a) Az (a) + > f.(&) A; (&), 
Aout (LA Dit h(a a) Acut, + >. mie x) Aout, x 


For the proof of statement (i) we start with the identity (15) 
(12) far favi (©) f;(y)Kz Ky PA (a) A(y) = ay [ae fi (@)fe(y) Ke Ky TA(x)A(y) , 
Ke i 
= ——; — 2 n 
2 Oat, ‘i 


This relation is not self evident because there are simple examples of patho- 
logical fields which do not satisfy eq. (12). But using the fact that as a con- 
sequence of causality and spectrum conditions the vacuum expectation values 


(Q, A(m)--A(a,) F(A(@), A(y)) As) +A (Ym) 2) 


are boundary values of Wightman’s analytical functions (1°) a more detailed 
investigation (2°) justifies the interchanging of the x- and y-integration in 


| da fa ff) KLET(Q, A (ay)... Ap) T(A(a), A(y)) A) «+ A(Ym)2Q) - 


Then relation (12) holds for every matrix element of TA(x) A(y) because, ac- 
cording to the irreducibility of the operators A(x), any state vector can be 
written as a linear superposition of vectors A(x,)... A(#,)Q. 

Rearranging the integral on the left hand side with the help of Green’s 
theorem we get 


i| dyfs(y) Kr TA(æ)A(y = — ifaw =. Jura 400 2 he 
a A(x) Ag re, A surg A(a) ’ 
8 29(2) fla 
(ie = aa, ? g(a) = jt) Ho g(&) de. 


(18) In the sense of weak operator convergence (compare footnote (28), correspond - 
ingly for all following operator expressions containing time integration. 
fo} 
(19) A. S. WiGHTMAN: Phys. Rev., 101, 860 (1956). | 
(20) W. ZIMMERMANN: Order of integrations in reduction formulae, unpublished 
manuscript. 
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Carrying out the integration over # in the same way, we obtain 

far faurreo (y)K™ KY T'A(x)A(y) = AcupA: tin — A outg Aion À cute Ainp+ Aina À ing 
For the integral on the right hand side of (12) we have; correspondingly 

| dy DAT KF TA A 0) = A Asay — Ans Aig Asp Ais + Asia 


Inserting these expressions into (12), we get 


[as Anal [AE Asp) = 0 


With the corresponding relations between Ae, and Av, or Am, and Aout, 
1 in in 
statement (i) follows: 


LA (x), Api (y)] wr FA (x), A (y)] . 


The second statement (ii) may be derived from the identity 
(18) fa faurien RE RELA) A) AG) = 
= av Jace: (w)fa(y) Kr Ky TA(x) A(y)A(e) , 


which can again be rearranged with the help of Green’s theorem. The final 
result is 


A(z NASA inp) = A outa A vp ]A(e) 


Using statement (i), we get 


[En Arg) A(e)] = 0, 


and this shows that [4,,,, 4,,,] is a c-number because we have assumed that 
the A(x) form an irreducible operator ring. With the corresponding results 
for [Aj,, Ait], [Ass Aingl; We have proved the statement (ii). 


2°2. Bound states. — In the next step we want to derive commutation 
relations for the field operators Bout(+, £). We begin with the proof of the fol- 
lowing statements: 
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(i) The commutators of the incoming and outgoing fields coincide: 
(14) [Bone (5 È) Bou (Y PI = [Bn (®, | È) By, n]. 
(ii) The commutator of the incoming field is a e-number: 
(15) [B,, (a, €) Bn (y, 1)] = (2; CB, (@, €) B. (x, n)19) . 


We expand the operators B and Bet with respect to a complete orthonormal 
system {F (x)} of the positive frequency solutions of (O — M?) F(x) =0: 


Bia, &) = F(a) BY (0, €) + > Py (@) By (0 5); 


Bout (x £) => FA x) Bou, E) + > A a) Bot, aie) N 
The coefficients are determined by (?!) 
Bt(a,, €) = F i414,” Bw E) sE Fi) — FF) et Bix Al 
x 0) 3 Il , C2 x = 0X5 , f ’ 
Fi(x) = FP, (2), Fy (@) = Fz (x), 


(and similarly for Bos , (&)). 
As a consequence ‘of causality and spectrum conditions we have the identity 


(16) Jarfaur (x) F;(y) EX KXT A(« + €) A(x — É)A(y + Aly — 1) = 


: | dy | dn F* (0) Fy(y)KAKYT Aw + 6) Ale —HAY + mA — n) - 


After integration over 2, and y, this relation yields 


[By But (n)] — CB; (§), Bit ()1=0, 


out A 


and similarly for [B*,,(), Bé,(m)l. So we have proved that 


[Bou (0, €), Bout (Ys n)] = CB, (€, &), BY; n)] - 


(21) If & is spacelike we have 
Bi (ates E) = By (Xo. §)* - 


This relation is no longer true if £ is a timelike vector. 
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Statement (ii) may be derived from 
(17) | da | dy F*(0) Fy EM KITA + 8)A(e — AY + MA — MAG) = 
= (ay | do P (0) Fyly) KEKE TA(E + #)A(e — A (y + mAG— MAG); 


which yields 


[Bio (€) Bus (M)] A(2) = ALB. (€) Bre (M)] 
or 


[[B; (©) Bag (MIA@)] = 0. 


Therefore the commutator [B,,(x, £) B,,(Y, 7)]is a c-number (statement ii). 
Now we can determine the commutator (15) by calculating the vacuum 
expectation value. For this purpose we determine first the matrix elements 
of B, (x, £) between the vacuum state and an arbitrary state vector. 
If ®, is an eigenstate of P, belonging to the eigenvalue k, and the rest 
mass — ki = M? we have 


_ exp [tke] 
(2, TA(x + E) A(x — E)D,) = Co = Fx(É), 
with 
(18) F,(£) = (2x) (2, TA(E) A(— €)®,) 
and 


(0, + M*)(Q, TA(a +&)A(a—£)G,) = 0. 


From this and the definition (3) of B,,(#, €) it follows that 


in 


ik 
(19) (2, B, (a, ED) = (Q, Bla, E) Dx) = eae F,(é). 


On the other hand if ® is an eigenstate belonging to the eigenvalue &, with 
the rest mass — ki M? we get from (5): 


(ki a M?)(Q, B,, (x, £)®) = à (Eee M?) (9, B,, (x, £)®) — 0 
hence 
(2) Boia E)D) = 0. 
Therefore we have 
ak 


(9, Bin (æ, Ë) By (y; n)Q) =| se 


(Q, B;, (a, E)®,) (D, PB; (y, 7) Q) = 


le fede . 
a oa | ae F,(É)F,(n) exp Like — y)], 
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where we have used the invariance under space time reflection for 
(20) (Dr, TA(n) A(—9)Q) = (2, TA(y) A(—0)®,) = Fam). 


Introducing the Fourier transform B,, (ky &), we have 


a 
Balt, £) = yoy [dak exp Likelo( + MBA, 8) 
and defining creation and annihilation operators 


Bi (k, €) = B,(k, &), for ky =+|Vie+ M|, 
B;(k, €) =B,(—k,é), for ky = —|Ve + MI, 


we get from (19) and (15) the final commutation relations 


(21) [Bx (k, €) By (p, m1] = Fa (6) F,(n)2ko Os (KE — p) ; 
(22) [B#(k, £)B=(p, n)] =0. 


Finally we only mention that the same methods used in this Section to 
prove the commutation relations (8) and (21) may be applied in order to derive 


(23) LA, (2), By (49) = 0. 


3. — Derivation of an S-matrix formula. 


The operators Bou (4, È) depend on the relative co-ordinate € of the bound 
state. For the S-matrix we are only interested in the center of mass motion 
of the bound states and want to carry out the limit £ > 0. Therefore we 
define operators Boxt(x) by 

Bout (x, é) 


(24) Bel) = Foe) 


F(E) = (220)! (Q, TA(E)A a E)D,) ’ 


1 


where ®, is the bound state at rest. 
In order to prove the existence of the limits € > 0 in (24) we divide (21) 


by F,(£) and get 
+ 
Bek, €) ,B 


= F,(n)63(k — p) - 
F, pe in(P, n) (7) ( P) 
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Differentiation with respect to £, yields 


0 Ba) 
e er, BnlP, )| = 9, 
LE x (€) i (P n 


and correspondingly from (22), (23) 


D Bak, €) ul 
pi AIAR 25 == 0 
a P,(£) » Bi(P; n| ’ 


a BEE x 
ES ae ap) SIE 


Thus (0/06,)(B; (k, E)/(F,(€)) commutes with all operators A,, (x) and B,, (y, 7). 


Since we have assumed that the operators A,,(#), B,,(y, 7) together form an 


irreducible operator ring (0/0£,)(B; (k, &)/F,,(€)) is a e-number: 


0 Bi(k,é) _ 0 (2, Bulk, 6)2) 
OF, Fate) > QE, F,(§) 


= 0, 


because of (6). Therefore the expression 


Bi (hy €) 


Tree BEE 
FE) Si 
is independent of £ and 
. Bin(k E) 
25 TRE k == lim = 2824 
( ) ( ) E>0 F, (È) ed 


does exist because 


Fg) _ P(RE, ©) 
PE) 0) 


is in the limit éÉ +0 independent of k: 


d RIO 
Dica 


From (21), (22), (23) we get: 
[Bk (k), B;,(p)] a 2ko ds (K Pp); 


[B;.(k), Bi (p)| = 0, 
[Bi (E), A, (y)] = 0. 
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Hence 
(26) [Bow (w), Box (y)] = iA(M, x—y), 
(27) [Bout (x), Aout(y)] = 0. 


Further we have 


(0 — M?) Bowi(4) AI) 
and 


exp [aka] 


(9, Bou (v) Dx) = = CEE ? 


k+ M°=0, 


if ®, is the bound state with four momentum k,. Bou(x) is given in terms 
of A(x) by 


(28) Bow(4) = = im ALZI, + É)A(AC — &) + da’ Asay (M, æ— x')J(x', 8! 


o Fy (€) 19 


According to (25) B=(k, &) is a multiple of B,,(k). From this fact we con- 
clude that the operators A,,(x), B,,(y) together already form an irreducible 
operator ring. Therefore the whole Hilbert space Æ can be built up by the 
creation operators 


Agu(k), = Baw(k) 


of incoming or outgoing elementary particles or bound states of momentum k. 
The state vectors 


(29) Dirt = Ata (hy). Aj, (hr) By, (By 1) ++ Bin (Bn) 2 
as well as 
a sia a, = Ani (k, Je Ai (ki) Boat (ki41) ie Bla (kn) Q, 


(30) 
(k2=—m? for i<l, ki =— M? for i>l) 


form a complete orthonormal system of Æ. The S-matrix is defined as the 
operator transforming the incoming into the outgoing states: 


(31) CA Spr fl (Dr Pr pu” ©) ; 


Since ‘both systems (29) and (30) form a complete basis of the Hilbert space 
the S-matrix is unitary. 
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It is possible to express the S-matrix by the vacuum expectation values 
of T-products only. To show this we derive from (7) the following reduction 
formula (?*) 


(32) [ST(2, … ox), Ad] = + à | SE (MET ST (0, Mi) ; 


(33) [ST (#,... dx), Bia(0)] = + ifasrs (2)K*OT (a, rs 2@—C, 2+ 6), 
with 
fi (e) = f, (@), 


T (ay...) = TA(a,) -.- Alta), ey, 


If we now insert plane waves instead of the wave packets f(x) we get 


ie(k ° m 
(34) [ST(c,... dr), Ain (k)] = — aot de exp [ike]K7 ST (a, ... U2), ° 


(35) [ST (a, ... 0), Bi (k)] = 


a & tke) KS Dori 
earn de exp [ike] HST (a, ... te, 2— 6,2 HE) 


Since the state vectors (29) form a complete basis S can be expanded with 
respect to the incoming fields: 


© _ 4@\m+n 
136). 8 ==> Sera dE Ae idl well peck lege lido 


Apolo 10) 
. LL O(k; + m°) i O(k? + M?) : Ain (K,) + Ain (Km) Bin(h) -- 9 3196 LA ces 
t= j=1 
where the coefficients are 


i QU'A QUO Em; a eee ln) == 


= (hy) ... e(En)e(h) … &(In)(Q, Le LS, AX(K2)] … A (Km) | BE (l)]... BAU) 12) = 


1 
~ Coin, (G) … ai di, eee dyn de... de, exp. [t(D eye + dle) > 


Re. RE RE KS (OD, TY Yn ya Gyo Gas Ba La a eee 


(22) All reduction formulae in H. LEHMANN, K. SYMANZIK and W. ZIMMERMANN: 
Nuovo Cimento, 1, 425 (1955), Sect. 2 are contained in the single formula (32) which 
was given by K. SyMANZIK: unpublished. 
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The last line of (37) follows by iterating (34), (35) and taking the vacuum 
expectation value. Transforming (37) into co-ordinate space we get (°°) 


= x (= ue 4 1! 7m >M Tr M 
S 2. ai lim FASI di 07, eg de KR 


CO Tu dr Cc Cy a 28 + Sy Su Gos Sn — €)Q): 
a Ain (Y1) eee Ain (Ym) Bin (21) VOL Bin(@n): 2 


(38) 


The function F,(£) which enters in the expansions (37) and (38) can easily 
be expressed by vacuum t-functions (i.e., vacuum expectation values of 
T-products) if we take the vacuum expectation value of reduction formula (35) 
Tor ba 29) == Fm, > — È: 


(39) FE) Fey) = — à | dy exp [iky]E® (0, TA(OA(—2)Aly + mA — 2) - 
Putting £=»n we get 
(40) FE =— 4 | dr exp [ike]K¥ (0, TA(E)A(—£) Ale + Al — 90). 


Inserting (40) in (37) or (38) the S-matrix is completely given by the vacuum 
t-functions. 


4. — Local field operators for bound states. 
In this section we want to take the limit & + 0 of the operator 
B(a; €) = TA(x +£) A(x — 6) 


itself. We assume the existence of 


| TA@+ HAW — 8) — (0, TAA 82) 
Chi BOLT Cm (0, TA()A(- £)®) 


(23) One gets corresponding expansions of the T-product replacing S by ST(a, … 7) 
and 
TY, ve Yms a + Li» . Sn + Cn: cia C1 sea €n Ca) ’ 
by 
TR oxo Lys Us ne Uns À + Caofsa3 Ont ons ey (ROTA ee pre 


in the expansions (37) and (38). 


ns 
19 
Del 


bui 
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which in a formal sense can be written as (24) 

(42) B(a) = a + A(z)? -- a ?b 

with the two (probably divergent) renormalization constants a and b 


| a =—ilexp[—iMx,](Q, TA(0)?A(a)?Q) dx, , 
(43) | 


(according to (40)) 
b 2, A(0)Q) 
The operator B(w) is, of course, invariant under the inhomogeneous Lorentz 
group and commutes with B(y) and A(y) if (e—y)?>0. Furthermore we 
get from (7) or (28) the asymptotic conditions 


(44) lim (®, B,(t)¥) == (®, Bout, P) 
td +o 
with 
. 5 
B,(t) = — i[d3xB(2) ta F*(x), 


xy=t 
for any normalizable solution of (LJ, — M?) F(x) = 0. 

Therefore the field operators A(#) and B(x) together satisfy the three prin- 
ciples of invariance, causality and asymptotic conditions exactly in the same 
formulation which was given in an earlier paper for elementary particles 
only (2°). From these principles we get in the usual way an expansion of the 
S-matrix: 


~ CA Anes 71 >M 7e 7 
(45) Se > RESI dy, ese dYm dz, eee dz, Ky; "00 ca à eee Ki - 
mn=0 he 


TICA TA(y;) ..» A(Ym)B(é1) «.. B(2n)2) An (Ys) An (Yn) Bis CD Beal oo 


which may also be obtained directly from (38) taking the limit È > 0. 

We remark that there remains at least one formal difference between the 
case of two independent elementary particles and the case of one elementary 
and one composite particle. If B(x) belongs to a bound state composed of 


elementary particles described by a field A(z), it is possible to represent B(x) 


(24) For the case of A4-coupling at least in perturbation theory the matrix elements 
of B(x) can be made finite to any order of the coupling constant if a and b are deter- 
mined by (43). 

(5) H. LEHMANN, K. SyMANZIK and W. ZIMMERMANN: Nuovo Cimento, 1, 425 (1955). 
Although there only the case of one scalar field was considered the generalization to 
several fields is obvious. 
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—roughly speaking—as a polynomial in A(x) (in our model Eq. (42)). We 
may impose this relationship as an additional condition beyond the principles 
of invariance, causality and the asymptotic condition in order to exclude the 
case of two elementary particles. 


5. — Examples. 


In this Section we shall give simple examples for the expansions of the 
S-matrix which we have derived in the last sections. Instead of the model 
of one scalar field we consider the more interesting case of a charged scalar 
«nucleon » field y(x) interacting with a neutral scalar «meson » field A(x). 
We assume that there are just three non-vanishing discrete eigenvalues 7°, 
x? and M? of — ie m? belongs to the one nucleon states, x? to the one meson 
states and M? to a stable bound state Ÿ of charge two and spin zero (« deu- 
teron » state). As a consequence of charge conservation 


(Q, y(æ)) = ( 


vanishes identically for a deuteron state Ÿ. Under the corresponding assump- 
tions, as formulated in Sect. 1 for the case of one scalar field, we can apply 
the methods developed in the last sections. As an example for the S-matrix 
expansions (36), (38) and (45) we consider the special scattering process with 
one deuteron and one meson of momentum q, resp k, in the initial state and 
two outgoing nucleons with momentum p, and p, in the final state: 


1 : 
(36/) oe; Spi) = apne |i dx, dy dz exp [i(qz + ky — prt, — Pata)] 
qa 


KEKE KG KY (Q, Ty(ar)y (2) A(y)p(e + pe — 09), 


i a 
= lim ———_—_ | dr, da, dy dz ex 2 + ky — PX, — Pot) |” 
(38" Ln ro l r, da, dy dz exp [i(q2 + ky — pits — Pots) 
KP KM REA (Q, Ty(asp(@a) Ay) ple + pe — 2) 2) , 
= | da, dx, dy dz exp [i(qz + ky — Pas — Pato)] 
(45°) 


MEA Re dt, (AA Ty (x1) p (#2) A(y) B(2) 2) 1 


pi=pi=—m, g=—M?, k?=—x*. 
The function F,(¢) is given by 


E(t) = (2x) (2, Ty) y(— OD) = 27 (G,, Te (- È) 2) 
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and satisfies — 
oC)? = im [a exp [ige] FK" (9, Tp(S)p(— Sy + Cyp(e—2)2) . 


The deuteron field operator B(z) is defined by 


_ Ty(e + C)p(z— 6) 
Ru F,(t) 


where F,(£) belongs to the deuteron state W, at rest. In a formal sense B(z) 
may be written as 
B(e) = e ple)? 


with the renormalization constant 


c= (2x) (9, y(0)2%)? = — tf exp [— iM] (Q, Ty(0}y(x)?Q) dx . 


Finally we remark that expressions (36’), (38’) and (45°) hold also if the 
model contains additional stable bound states. In this case the incoming 
nucleon, meson and deuteron fields are, of course, not irreducible in the whole 
Hilbert space, and the S-matrix expansions contain in addition the incoming 
fields belonging to the higher bound states. But if the state vectors OS 
and @”;”: contain only nucleon, meson and deuteron states, expressions for 
the S-matrix elements (@";", Os") can be derived from the asymptotic 
properties of Aow, Wout and Bout in a similar way as in Sect. 3. 


* * * 


I would like to thank Professor OPPENHEIMER for the kind hospitality of 
the Institute for Advanced Study, and I am grateful to the International 
Co-operation Administration, Washington, for a grant. I am indebted to many 
physicists in Gòttingen, Princeton and Berkeley for helpful discussions. 


RIASSUNTO (*) 


Si esamina un campo scalare causale e invariante coinvolgente uno stato stabile 
legato. Si deriva una formula per la matrice S e si dimostra che lo stato legato si può 
descrivere per mezzo di un operatore di campo locale e invariante. Per semplicità si 
considera solo il caso di particelle di spin nullo e di stati legati; l’estensione ad altri 
casi è tuttavia possibile. 


(*) Traduzione a cura della Redazione. 
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The Eigenvalues of a Randomly Distributed Matrix. 


R. ENGLMAN (*) 
H. H. Wills Physics Laboratory - University of Bristol - Bristol 


(ricevuto il 4 Agosto 1958) 


Summary. — The elements of a symmetric matrix are supposed to be inde- 
pendent random variables, with a distribution law that depends on the 
distance of the element from the diagonal. What is the distribution of 
the eigenvalues of this matrix, in the limit as its size becomes infinite? 
A method, based on the resolution of the symmetric matrix into two 
triangular ones, is developed, which reduces the problem to a form which 
has the promise of being soluble by numerical methods. 


1. — In 1953 F. J. Dyson (1) gave a method for the calculation of the distri- 
bution function of the eigen-frequencies of a disordered linear chain, in the 
limit when the chain becomes infinitely long. By a disordered chain is meant 
«a chain of one-dimensional harmonic oscillators, each coupled to its nearest 
neighbours by harmonic forces, the inertia of each oscillator and the strength 
of each coupling being a random variable with a known statistical law. » 

An indirect effect of Dyson’s method was the catalysis of a spate of papers 
purporting to generalize his results to three dimensions, mainly with a view 
to obtain the frequency spectrum of an imperfect crystal. The most recent 
list of these papers is, to my best knowledge, given by MARADUDIN, MAZUR, 
MoNTROLL and WEISS (?). Valuable as these works are they are not accu- 
rately described as « generalizations » of Dyson’s solution, since this was exact 
whereas most of the more recent works are effectively perturbation treatments. 

The work presented here also grew out from a study of Dyson’s paper. 
It is a generalization of that in the sense that an exact «solution » is given 


(*) Present address: Physics Department, The Technion, Haifa; Israel. 

(1) F. J. Dyson: Phys. Rev., 92, 1333 (1953). 

(2) A. A. MARADUDIN, P. Mazur, E. W. MONTROLL and G. H. Weiss: Rev. Mod. 
Phys., 30, 175 (1958). 
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for the frequency spectrum of a chain, whose constituents are coupled to any 
number of neighbours according to known statistical laws. Since a finite piece 
of three dimensional lattice can always be «mapped » onto a one dimensional 


one, this can also be fitted into this scheme provided we choose the law of 
coupling suitably. Nor is our interest confined to the calculation of the fre- 


quency spectrum, but rather we give a programme for calculating the eigen- 
values of any symmetric matrix whose elements take values according to some 
statistical law. 

In this connection we have in mind, among other things, the distribution 


of electronic energies in a disordered crystal, a problem in which, in spite of 


its importance, only limited progress appears to have been made. 
In the above paragraph prudence and honesty compelled us to write so- 


lution in quotation marks. Stated more exactly the claims of this paper are as 


follows. 

When each of the constituents of a one-dimensional chain is coupled to 
a small number of its neighbours, the distribution of eigenfrequencies may 
be derived exactly from the solution of a system of integral equations. The 
number and complexity of these integral equations grows rapidly with the 
number of neighbours coupled to each constituent of the chain. These equa- 
tions are derived by a very simple and general method based on matrix theory. 
(There is a cryptic remark in Dyson’s paper bearing on the possibility of such 
a method. DEAN (*) comments on this remark). In particular, the reader 
will note how readily some of Dyson’s results are obtained in the simplest case 


‘of nearest neighbour interaction. 


On the other hand, when each constituent is coupled to many others 
(possibly to all), or in matrix language, when elements far off the diagonal 
are non-zero, then the analytical solution of the system of very unpleasant 
integral equations is clearly out of question. Indeed, in this case (and three 
dimensional problems are of this kind) we had better not speak about a solu- 
tion, but rather a reduction of the problem to a form amenable to numerica- 
work. We are, however, quite serious about this aspect of the problem and 
propose to carry out with God’s help detailed and extensive computations as 


soon as this is feasible. In the meantime a method of procedure is suggested 
tentatively. 


2°1. — We consider the eigenvalues H, of the Nx N symmetric matrix A 
(elements nn — @nm), or equivalently, the roots E, of the determinantal equa- 
tion 

re 
(1) 0 = det A4'(E)=][(E,— E), 


i 
ei 


(3) P. DEAN: Proc. Cambr. Phil. Soc., 52, 752 (1956). 
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where 


A'=A—EI. 


In most physical situations that are of interest to us, it is sufficient for 
the purpose of finding the energy eigenvalues, to consider a real symmetric 
Hamiltonian matrix, A. In the most general case A would be complex and 
Hermitian. Then our treatment would proceed in much the same way as for 
a real matrix but at the expense of typographical complications, which would 
inconvenience printer and reader alike. 

As N — co and the eigenvalues E, are expected to become densely distri- 
buted we define from (1) the functions R(E) and D(E) of the variable £, 


(ce) 


(2) R(B) = 5 log det. AB) = + D 10e ( E) = tox (8 E')D(E')d£'. 


D(E) thus defined is the density of eigenvalues in the neighbourhood of £. 
It was shown by Dyson that equation (2) may be inverted so that provided 
R(E) is known as a function of the complex variable E, D(E) may be found 
too. This means that we need only apply ourselves to the task of finding 
R(E) and this is in fact what we propose to do now. 


99. — The NxN matrix A’ will now be resolved into a product of an 
(Nx N) lower triangular T and its transpose T", 


(3) API. 
+ + + + + + + + + + + 
+ + + + + + + + + + + 
+ + + + + = + + + + + + 
+ + + + + + + + . + + 
+ + + + + + + + + + + 


The elements t,, (Æ0 for r<s) of T are determined by the set of 
4N(N +1) equations 


(4) Gn LE A 
r=1 
(5) Ann = > tartme 9 (n << m). 
ral 


40 - Il Nuovo Cimento. 
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The determination of t,, as a function of E will be seen to be unique by 
starting at the left upper corner of A and working one’s way downwards (and 
inwards) using at each stage 


n-1l 


(4.a) Inntmn = Amn — > tar tne 
rT=1 
or 
m= 
(5.a) ? = Gun —E—Dt,- 
r= 


In this line take Ret,,(E)> 0 for E > — 00.] 
Now by (3) 


det A’ = (TI be} 


so that in (2) 


(6) R(E) = (2/N) > logt,,(E) . 


Tf we now consider a matrix A whose elements a,, take values according 
to some distribution law, say f.,(a,.) (r<s) and g(a,,), then ¢,,, will also 
follow a statistical distribution law, say Fnn(tmn) and Gyn(tun). When the size 
of the matrix increases indefinitely then, as n > oo, G,,(t) is expected to 
converge to a value G(t) (*). In this case it is permissible to rewrite (6) as 


— 
1 
— 


LU) 2 frog 160 a : 


The integration here, as everywhere in the sequel, is over the full range of the 

distribution function. 
The reader is reminded that F,,,, Gn, and therefore also G are all func- 

tions of E through equations (4a) and (5a). In order to avoid crowding of 


symbols we suppress this dependence of E, which should, however, be borne 
in mind. 


3. — The determination of G(t) will complete, through the use of (2) and (7), 
our task of finding D(E). First we find G,,, Fam. The quantities on the 


(*) This convergence seems natural on physical grounds. Let it be stated, though, 
that it has not been proved directly, from the equations defining G,n: 


5 


THE EIGENVALUES OF A RANDOMLY DISTRIBUTED MATRIX 619 


right hand side of (4a), (5a) are statistically independent. Therefore 


RL 


n—1 
Eyn(e) = [fa y Ge (y) fa (2y — > di à) II dx dy, Le, ) Pg Ged ; 
j=1 r=1 


n-1 n—1 
Gorle) = 2f.. ses + B+ > x) IT dx, F,.(x,) . 
j=1 r= 
When N and »-> oo, then as noted earlier, we assume that Gi, > & 
and that F,, depends only on the horizontal distance r—s from the dia- 
gonal. Also we relabel the distributions F,, of t., (which have so far been 
denoted with the upper left corner of the matrix as the origin) as F7% with 
the point nn as origin of reference. Implicit in this notation is the assumption 
that as we move away from the diagonal F,n(y) will approach the d-function 
ò(y). The rapidity of this convergence will be dictated by the speed with which 
the off-diagonal matrix elements a,, will tend to a 6-function distribution as 
we move away from the diagonal. 
The above two equations can then be rewritten as 


(8) F(z) "il w | 


and 


… dy-yf,.(2y + > eiyi) Gy) Il da, dy, FS +” (x) (ys) (r=1,2,...) 
j=1 s= 


(9) G(2) = fof eg(e? +H + > x) TI Feo (as) de, . 


4. — We take a linear chain first. It is understood, we repeat, that f, will 
tend, with r > co, to a 5-function and in consequence so will F™, too. 

Two cases will now be considered: 4) when the constituent interact with 
the first nearest neighbours and b) when the coupling is with the first and 
second nearest neighbours. 


a) This is Dyson’s problem, put in a more abstract forms. 
Now 


f(x) = (x) for r > 2. 
Solutions of the system of equations (8) and (9) are 
Fc) = (x) for r > 2. 


Using these solutions we can eliminate PF from (8) and obtain from (9) an 
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integral equation for G(z), 
ae al | gt? + E + a*)f(ey)y (y) dy der. 


This is effectively the integral equation (44) or (40) given by Dyson for par- 
ticular choices of g or fi. 


b) The solution of (8) and (9) are 
F”(y) = Ô(y) for r > 3. 


together with the coupled integral equations: 


“ne =| i f | fi(ey + a) folary’) E(y) G(y') F® (ys) yy’ dx, dy, dy dy’ , 


G(z) = 2e] | | g(22 + E + ar + 208) f(ay") G(y') FO (x, )yy' da, da, dy’ . 


It is clear that the problem rapidly becomes unmanageable by analytical means 
except for the simplest choice of f,. 


5. — In describing the interactions in a three-dimensional structure (of M? 
constituents) by a matrix, a possible way of ordering is to label the rows or 
columns corresponding to the lit-mj-+nk constituent (i, j, k unit vectors; 
I, m, n integers: 1, 2,..., M) by 


r=l1+mM+nM®. 


Let us adopt the same scheme for r of F and f, occurring in (8) and (9). 
Then the meaning of these equations is clear. Their solution is another matter. 

It is fairly clear that as any one of J, m, n, tends to M (which is itself very 
large), F° will tend to a delta distribution, since f, does so ex-hypothesis. 
On the other hand, the topology of the coupling in any dimension higher than 
one effectively excludes the possibility of quasi-analytical solutions cf the 
form a) or b). What is to be done in this situation? 

Evidently, the question of numerical solution must be raised. Now, I 
think that some adaptation of the Monte-Carlo method to (8) and (9) is pro- 
bably the best approach, and in fact, intend to explore this avenue further. 
For the time being I must leave the problem in its present formal, program- 
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matic form hoping that the method given here will attract the attention of 
those interested in this field. 


* * * 


. The criticisms of Professor M. H. L. PRYCE, F.R.S. and Dr. D. BOHM were 
of great value in developing this paper. 


RIASSUNTO (*) 


Si suppone che gli elementi di una matrice simmetrica siano delle variabili casuali 
indipendenti con una legge di distribuzione dipendente dalla distanza dell’elemento 
dalla diagonale. Quale è la distribuzione degli autovalori di questa matrice al limite 
in cui la sua dimensione diventa infinita? Si sviluppa un metodo, basato sulla risoluzione 
della matrice simmetrica in due matrici triangolari, che riduce il problema a una forma 
che promette di essere di risolvibile con metodi numerici. 


(*) Traduzione a cura della Redazione. 
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On the Helium Film in Non-Static Conditions (*). 


S. FRANCESTTI 


Istituto di Fisica dell’ Universita di Firenze - Firenze 


(ricevuto il 4 Agosto 1958) 


Summary. — The problem of the helium film under non-static conditions 
is treated in this paper. Conclusions regard essentially the effect of the 
(subcritical) flow and that of temperature inhomogeneities on the form 
of the (saturated) film. The existence of a kind of thermomechanical 
effect of the film is demonstrated. The theory is a substantially phenom- 
enological one, based on the assumption of reversibility. This assumption 
is briefly discussed in Sect. 4, while the irreversibile flow is taken into 
consideration in Sect. 2. 


Introduction. 


The purpose of the present paper is to develop a simple, essentially phenom- 
enological, theory of the non-static film (including of course the static film 
as a particular case). 


1. — The isothermal, reversibly flowing film. 


The fundamental assumption is that the subcritical flow in helium II is 
a reversible phenomenon. This point has attracted the attention of various 
authors, F. LONDON in particular (*), and finds its empirical support in various 


(*) Part of the present work was the object of a contributed paper at the Kamer- 
lingh Onnes Conference in Leiden, June 1958. 
(1) See for instance F. Lonpon: Superfluids, vol. 2, p. 149-150. 
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experiments, among which the double-beaker experiment by DAuNT and MEN- 
DELSSOHN (2) is a particularly telling one (*). 

The assumption of reversibility means that a suberitically flowing film is 
a system in equilibrium, somewhat as a current-carrying superconductor is. 
It also implies the introduction of an extra intensive variable describing the 
state of flow and as such one may choose the flow ® per cm width of the film 
(or any equivalent quantity). 

The form of the film can be determined by writing the condition that the 
chemical potential has to be the same everywhere in the system. In particular, 
it must have the same value at any point in the film as in the bulk liquid. 

To find an expression for the chemical potential u, let us consider the 
Gibbs function G,,,, referring to an element of the film of unit area and thick- 
ness f. at some height z. The volume V of this system is ¢ itself, its mass is 
m = ot (where o does not differ appreciably from the bulk liquid density, at 
least for usual, saturated films). 

Calling F,,, the Helmholtz free energy, one has 


film 


(1) G = Pan pV= Ta + pt, 


film 
where p is the pressure in the vapor phase. 

F,,,, is meant to include the contribution otgz of the gravitational potential. 
On the contrary, it is convenient to write separately the corresponding term 
in the expression of G,,, by which symbol we mean the Gibbs function of 
an equal quantity of helium as considered above but placed at some height 
2 inside the bulk liquid. By definition 


(2) Csi = Es + otgz + Pat a By n (09% + Pat 1 


where p, is the hydrostatic pressure in the liquid. Clearly, eg?+ pn is a con- 
stant, and taking, as usual, the origin of 2 at the free surface of the bulk liquid, 
we have 0g2+p,= p for 2 = 0. Therefore 


(3) Giai = Eve Sr pi ° 


Of course, Fun is simply proportional to the mass of of helium and can there- 
fore be written 


(3’) Fonte = te ot , 


where f,n no longer depends on t. 


(2) G. Daunr and K. MenDELssonn: Nature, 157, 839 (1946). 
(3) The same experiment has been recently performed with special accuracy by 
J. F. ALLEN (communication to the Kamerlingh Onnes Conference). 
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To get the chemical potentials 4, and Myx» We have only to resort to the 
general formula 


ag 
: “(Fla 


that is, we have to take the derivatives of expressions (1)vand (3) with respect 
to m= ot, holding the intensive variables constant. We get 


à 1/0 Fin p 
(5) Lim 0 at À, n ’ 
1(0Fru p 
(5') Mou | a Nés È cs four x 0 
By equating we have 
ie 
(6) 0 pri (Fim Da Fu) = 0 , 


which, apart from the factor 1/0 (4), is just the equation employed in discus- 
sing the profile of the static film in some previous papers (5). 

The same equation must hold for the reversibly flowing film, provided the 
proper changes are made in the expression of F’,,,,. These changes amount 
to the addition of a kinetic energy term 


(7) Ke, = 40,051. 


(Of course, it is just the character of reversibility which ensures that the 
kinetic energy of the flow makes a purely additive contribution to the free 
energy or the Gibbs function). 

This positive term in the expression of G,,, gives rise—however—to a 
negative one in the expression of w,,,. This comes about because by taking 
the derivative of (7) with respect to ot, we have—according to definition (4)— 
to keep the flow ® constant (°). Assuming, as usual, that v, is the same throug- 


(*) The reason for not dropping the factor 1/g is to preserve for the two terms 
in eq. (6) their homogeneity with the chemical potentials, from which they differ only 
by the additive term p/o. 

(5) S. FRANCHETTI: Nuovo Cimento, 4, 1504 (1956); 5, 183, 1266 (1957), henceforth 
referred to as I, II, III. Although not explicitly stated, it is easy to see from the 
arguments employed (especially in I, p. 1507-8) that Fim includes the gravitational 
contribution, while Fu does not, as assumed above. 

(°) Alternatively, one can say that the equilibrium situation is characterized by 
the minimum of the Gibbs function (for the entire system) under given conditions 


of p, T and flow. The final equations are of course the same as with the procedure 
followed here. 
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hout any (horizontal) section of the film, this means a condition 


(8) vst = const 0 ; 


Os 


(expressing © in em*(cm s)-'). Eliminating v, between (7) and (8), one gets 


! ï | 2 
(7) En = 5 = Dai, 
Us 
from which a contribution 
A il 
(9) im it (Uri) stat ==. FP 9 n) EE 5 


to expression (5) and to the left hand side of eq. (6) is immediately derived (7). 
As is well known, the subcritical flow ® cannot exceed a maximum value 

D... Which is a function of the temperature. We therefore have 

(10) D = HDs 0<q <1. 

On the other hand, the experimental evidence seems to be compatible with 

the assumption that the maximal flow D... is given by 


(11) Da = a Din 4 


where ©,,, is a limiting value depending only on the kind of the wall. 


Through (10) and (11) we get an alternative expression for [gin — (Him)stat > 
namely 


NIH 


(9') Pim — (iim) stat as > 1? Dist sy 


We have therefore as equation determining t in the presence of a flow 


s ; Ais: 
re Dim t 3 0 DI (Fim = Fyunc) stat = 0 . 


ple 


(12) oF 


(7) We have entered these details, which seem fairly plain, because they seem to 
have been overlooked by some author. In particular one should carefully distinguish 
between the kinetic contribution to and the specific kinetic energy, the former being 
just the negative of the latter. 
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The diffrence F,,,-- ya, in the static case has been carefully investigated 
in the papers quoted in ref. (5). It is due on the one hand to the fields acting 
on the film (Van der Waals field of the wall, gravitational field), while on 
the other hand there are terms due to the smallness of the transverse dimen- 
gion of the film. The latter are two: a temperature independent term relating 
to the zero point energy oi the liquid, regardless of «its state of excitation 
(« Atkins’ term ») and a temperature dependent one due to statistic effects on 
the excitations themselves. These «limitation » terms are the only ones re- 
quiring atomistic considerations to be worked out. The choice of the model 
employed should however have only secondary importance. 

The limitation terms in F,,,— Fy, contribute to the left hand side of (12) 
(or to the expression of y,,,,) terms in ¢*, while the Van der Waals term is 
one in t-* and the gravitational potential is of course simply ge. 

For details the reader is referred to the papers quoted in ref. (°). However, 
the Van der Waals term can immediately be derived—in a somewhat sim- 
plified way—by noting that the Van der Waals energy per cm* of the film 
will have the form 


t 


—furvar. 


Lo 


where x is the distance from the wall, x, some (inessential) minimal distance 
of approach and a is a positive constant depending on the nature of the wall. 
The contribution to &,, due to the above term is its derivative with respect 
to ot, that is 


film 
(13) Mar 
Q 


If A is the constant in the expression — Ax-* for the potential energy of a 
helium atom at a distance x from the (plane) wall, a/o turns out to be 
(1/m)(A — 2A,,,) with m the mass of the helium atom and A,, characterizing the 
action of a fictitious «helium wall» (having the density of liquid helium). 
In other words 


(13') ee = (A= 245, = 2:19-10®(A/347), 


Some approximate values for A calculated according to SCHIFF (8) are 
listed in Table I (units are c.g.s.). | 


(8) L. I. Scurr: Phys. Rev., 59, 839 (1941). 
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Referring for instance to eq. (19) of paper IIT, one finds 


(14) 


TOT 


lo DI 
A (Finn — Four) stat = — o (at-* + bi?) + gz, 


with a given by eq. (13) and 


(14') b =[1.87 + R(7)]-10° (c.g.s. units). 
TABLE Je 
o ===" =" SS 
Bai apreres (scat ie hs pores Us. à ANS 10 ES | 
| Classe eat tito 43 » | 
Miata WE LR, ee | 1.7 0» | 
M a A, 0.97 » | 
Solid hydrogen or neon . . . | 0.5 » 
jd Liquid helium» . . . . . . As = 017) 


The detailed form of F(Z) is rather laborious to get and can be deduced 
in two instances from paper III. Table IT shows some approximate values 
for a copper wall. 


TABLE II. 
SP | R(T) T R(T) 
<1:0 negligible 2.00 0.82 
1.20 | 0.017 T, 1.49 
1.40 | 0.064 | 2.20 1.52 
| 1.60 0.17 2.50 1.56 
1.80 0.40 3.00 1.47 


Inserting expression (14) into eq. (12) one has 


al 
1% gi, tl (at +00) + 92 = 0, 
20 0 


or else 


(15) a + (b+ on Pit = 094°. 


Fig. 1 suggests à graphical way of solving this equation with respect to t 
(by given 2). It also shows an approached but quick evaluation of the flow 
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effect. Indeed one has 
RO 
—~ À 3ogzt?, Ay © 5 21 Din € , 


and therefore. 


At Os 2Dîm 
(16) rep Fag a ee 
t 6Goget? 


Visibly, the effect is an increase 
of the thickness for the flowing film. 
Intuitively, this is comprehensible, 
because while thickening the film 
decreases its free energy by reducing 

Pig A. the amount of kinetic energy asso- 

ciated with the given flow. (If at- 

tention is not paid to the considerations in footnote (7), the opposite result 
obtains.) 

Since 21° is an increasing function of 2, the effect is larger at lower heights. 
It is also larger at lower temperatures, because of the factor 0,/o which, more- 
over, makes the effect vanish at the A-point. 

In most experiments the flow has its «saturation » value, that is 7=1, 
while the temperature is sufficiently low to have 0/0 & 1. A typical value 
for D,,, on a clean metal or glass surface is ~ 10-4. With et? — 2-10-!! which 
is also a typical value in ordinary conditions, one sees from (16) that an effect 
of the order of 8% can be expected. 


2. — On the irreversibly flowing film. 


In all rigour, the treatment of the foregoing 
Section refers to a truly reversible flow, such as takes 
place on the walls of the inner beaker in the double- 
beaker experiment (flow without level difference). 
The situation is however somewhat different when 
the film is flowing from a higher to a lower level, 
as Sketched in Fig. 2. 

Here the chemical potential of the liquid to the 
left is higher than that to the right by an amount 
g(è; —2,). This means that the chemical potential 
cannot have the same value all along the film. 


© 
œ 
1® 
De) 
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If 2,—<, is sufficiently small, the loss of potential energy can be counter- 
balanced by a gain of kinetic energy without exceeding the limits of the sub- 
critical flow: this is what happens in oscillating film experiments. Excluding 
this case and considering instead a stationary flow in the conditions just de- 
scribed, we have clearly to do with an irreversible phenomenon in which the 
potential energy g(2,—2,) per unit mass of the flowing liquid must eventually 
be turned into heat. (This does not mean that the flow is non-isothermal, 
since—presumably—that small amount of heat will get absorbed by the wall). 

One feels, however, that at least in regions away from point B, the form 
of the film should be determined by the amount of flow, independently from 
any other circumstance accompanying it. 

Thus, one is led to think that the irreversible aspect of the flow is probably 
confined to a «transition zone» around point B, where the change of y (which 
will not be definable inside this region) and the « matching » of the different 
thicknesses of the ascending and descending film are to take place. 

Indeed, in that region there will be à « squeezing» of the flow, that is an 
acceleration of the fluid. This might well be the cause for friction setting in 
and causing the conversion of mechanical energy into heat. (It is hard to see 
how this could happen in the lower regions, where the thickness of the film 
is continually increasing). 

A proof that the squeezing is an important factor may be seen in the 
instability that arises when the difference in the thickness on the two sides 
is comparatively larger, that is when the upper level lies sufficiently near to 
the rim, leading to drop formation and increasing the flow above its normal 
value (°). 


3. — The non-isothermal (resting or flowing) film. 


In the non-isothermal case too, we shall require the chemical potential to 
be the same everywhere and hence 


(17) du = 0,, 


all along the film, which means that we shall treat the film as being in equi- 
librium. (This point is briefly analized in the next section). 

It is hardly necessary to recall that the same assumption is the basis of 
a well known method of deriving the fundamental equation for the thermo- 
mechanical effect between two vessels connected by a superleak (1°). The dif- 


(9) L. C. JACKSON and D. G. HensHaw: Phil. Mag., 41, 1081 (1950). 
(19) See F. LONDON: Superfluids, vol. 2, p. 70. Or else, L. Tisza: Phys. Rev., 72, 
845 (1947). 
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ference is that in the latter case a pressure gradient is originated as a conse- 
quence of a temperature gradient. In the film this is impossible, and in place 
of a pressure gradient we shall find a deformation of the film, which may also 
be accompanied by a flow and constitutes a sort of isobaric thermomechanical 
effect of the film. 

The chemical potential # at any point of the film can be thought to be a 
function of the variables t, 2, p, T and O=7/(0,/0)®,,,. We shall consider 
the (external) pressure p to be a constant, and for simplicity’s sake we shall 
assume the same to hold for the flow ®. (Although in the presence of some 
evaporation along the film this may not be the case). 

Eq. (17) becomes therefore 


Ou Ou Ou 
ata de + ap dT =0. 


The third differential coefficient can be replaced by the negative of the average. 
specific entropy —s, getting 


Ou Ou 
LT = — de — = 
(17') = dt + >> de —sdT=0. 


The appropriate expression for ~ can be obtained from eq. (5), combining 


RI 


it with eqs. (5') and (14) and adding the flow contribution (9’). The result is 


il Lise el x 
(18) notes o at? AU + 2 0-1" Pi) + gz + + foun » 
The last two terms give no contribution to eq. (17') which becomes 
1 NE 
(19) — Da at + (o aL 5 0.1" Pin) ci dt + gd: — sd Tr =0. 


By integrating with respect to variables # and 2 from 2= 0 (where t= 00) 
up to any height +, one has 


1 i ili : OT 
——at-*——|b + on Di, | -? + gz —] s — de = 
0 Al g 057 i) = Jz pira Om 


0 


As the temperature differences will always be very small, the last term may 
conveniently be written s(7 — 7,) with s an average value and 7, the tem- 
perature at <= 
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Written in a form analogous to (15), the above equation becomes therefore 


4 


; I 
(20) a+ [> +3 anti) = 09 ei eta, 


which reduces to eq. (15) when 7=7, troughout. 
It follows from this equation, comparing with eq. (15), that a situation 
with 7 > 7; corresponds to have + changed by an amount 


8 


(20') ie wee 
g 


CEE 


meaning an upward sliding of the film. If instead one has 7 < Ty a cor- 
responding downward sliding takes place, both displacements being in accor- 
dance with the character of a thermomechanical effect (displacement towards 
the warmer region). 

The kinship of the present effect with the ordinary thermomechanical 
effect is also apparent when one considers the well known experiment by 
DAUNT and MENDELSSOHN (1!) in which two vessels containing liquid helium 
at two slightly different temperatures are connected through the film. Were 
it not for the distillation which tends to counteract the effect (see next section) 
a situation of equilibrium would be reached with the liquid at two different. 
levels in the vessels. Eq. (19) would apply. By integrating all along the film 
from the one vessel to the other, the t-containing terms vanish and we find 


g(22— 41) = 8(T,— Ti). 


Since, obviously, 09(% — 4) gives the « pressure head» AP, we have actually 
found the well known equation of the fountain effect AP = os AT. 

It follows also from eq. (20) that t would become infinity and then nega- 
tive if the square bracket in the right hand side should turn respectively zero 
or negative. This means that to have a stationary film, the condition 


(21) ST — Ti) < 92 
or else 
! ar at 


(11) J. G. DAUNT and K. MenpeLssoun: Proc. Phys. Soc., 63, 1305 (1950). 


3593 


632 S. FRANCHETTI 


must be fulfilled. If this does not happen, that is if the (upward) temperature 
gradient is too large, instability will set in. Instability in the presence of a 
temperature gradient was reported (and correctly interpreted as a sort of 
thermomechanical effect) by HAM and JACKSON (7). 

Directly from eq. (19) one may derive an expression for dt/de: 


(22) = =— Bar + 2 ( + 5 Ou Di) = i (9 — 8 =) 7 

We can see from here and from (21') that dt/dz is always negative, that 
is, we cannot have a «bulging» (stationary) film. Eq. (22) shows, however, 
that dt/dz can reach (negative) values as small as we please when s(d7/dz)(> 0) 
is made to approach «from below» the value of g. 

This might account for the large values found in some experiments (°°) 
for the parameter q in the expression t= const 2! for the thickness of the 
film, if we note that 


(22') Es 
q 


The observation following eq. (20), that the thickness of the non-isothermal 
film is the same as for the isothermal film at a point shifted vertically by an 
amount Az given by eq. (20') suggests the following approached expression for 
the thickness change At due to temperature inhomogeneities 


bay Cth is di\ s 
Ai a(S) “ae =—(S) 2 
| ). = sl g ba To) 6 


Combining this with (22') one gets 


(23) At PO (7). : 


(23’) Fe (At < 0). 


If At is positive (7 > T;), |At/t| is also bounded because of condition (21) 


(12) A. C. Ham and L. C. JACKSON: Proc. Roy. Soc., 240, 243 (1957). 
(13) See for instance K. R. Arxins: Proc. Roy. Soc., 203, 123 (1950). 
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which imposes 

At al 
(23”) EL SS (=4 e 

| t na q iso CES 


Since, normally, in the isothermal case q ~ 2.5 + 3.0, we see that the rela- 
tive thickness increase cannot exceed (30—40) %. 

The average temperature gradient (7 — 7,)/ associated with a given re- 
lative change At/t follows immediately from (23). Im Table III are listed 
some values for At/t= 0.1. They show the great sensitiveness of the film to 
temperature disturbances, in the upper part of the He II region (1). 


TABLE III. 

Ms | grad T ug grad T 
(°K) | (°K/em) (°K) (°K/em) 
0.0 | co 1.5 1.4: 10-41 
0.5 | 2.7-10-2 2.0 2.9-10-5 
1.0 | 1.6-10-3 Ta 1.7-10-5 


Suppose now that in some part of the film the maximum relative increase 
of the thickness has been reached. This means, according to (22) and (23), 
di/d: = 0, that is the film surface is nearly parallel to the wall. Suppose more- 
over that in a neighbouring (higher) region At falls sharply. Then a step-like 
profile results for the film. (A temperature distribution giving an effect of 
this kind would, for instance, be a « peaked » one, such as could be originated 
by some heat radiation falling on a narrow strip across the film). 

Steps in the profile of the film have been reported by DILLINGER (1) and 
by Ham and Jackson (ref. (!°)). 

If At is to fall from its maximum value to zero in a range of the order of 
one em, it is easily found from eq. (23) and the known values of s, that d 7/de 
must take values ranging from 107? to 10-4 °K/em (depending on the values 
of T and 2). 


4, On the nature of the equilibrium in a stationary film. 


The foregoing treatment has been based on the assumption that a stationary 
film is a system in equilibrium. That in certain instances this has to be con- 
sidered as a convenient schematization rather than a rigorous description, is 


(14) The entropy values utilized in the computation have been taken from 
€. F. SQUIRE: Low Temperature Physies (New York, 1953), p. 64. The very small dif- 
ference between Spun AN Sim is of course of no significance in the present connection. 

(15) J. R. DILLINGER: Conf. de Physique des Basses Températures (Paris, 1955), p. 97. 
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shown by DAUNT and MENDELSSOHN’s experiment, ref. (11), where the two 
masses of helium would have to be thought of as being in equilibrium insofar 
they are connected through the film, whereas they show manifestly to be not 
when the exchange of helium vapour is considered. Indeed, in a general way 
it may well look surprising that we can speak at all of equilibrium between 
bodies having different temperatures. 

The following considerations should help clarify hie matter. 

The key to the paradox lies in the observation that we would presumably 
have a truly reversible process in a-transfer of superfluid alone. Indeed, it is 
conceivable that the change from rest to motion—as long as the flow remains 
subcritical—reduces to an adiabatic (in the quantum mechanical sense) alter- 
ation of the fundamental state eigenfunction, requiring no rearrangement in 
the occupation of states and hence no entropy change. The same however 
could not be said of any transfer of the normal fluid (let alone of the vapour). 

There would therefore be true equilibrium between two helium bodies con- 
nected only via a channel rigorously impervious to the normal fluid. However, 
apart from the fact that it would be often impossible to avoid other con- 
nections (as that via the vapour in the previous example), the above condition 
would always mean an ideal case. Indeed, through any real channel the normal 
fluid would have to move, too, though much more slowly, bringing eventually 
the system to « ordinary » equilibrium, in particular, to isothermal conditions. 

We therefore see that the treatment given in the previous sections is an 
approximation valid to the extent that the transfer of the normal fiuid (or 
the vapour) is negligible against that of the superfluid. In many instances, 
this may well represent a very good approximation. 


RIASSUNTO 


Questo lavoro concerne il comportamento del film (saturo) dell’elio in condizioni 
non statiche. Viene studiato principalmente l’effetto del movimento (flusso suberitico) 
e quello di inomogeneità di temperatura sulla forma del film. Si dimostra l’esistenza 
di una sortà di effetto termomeccanico del film. La teoria è essenzialmente fenomeno- 
logica e si basa sull’ipotesi della reversibilità dei fenomeni. La legittimità di questa 


ipotesi viene esaminata nella Sez. 4, mentre il caso del flusso irreversibile viene preso 
in considerazione nella Sez. 2. 
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A Decay of an Anti-hyperon, A°, in Nuclear Emulsion. 
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(ricevuto il 4 Agosto 1958) 


Summary. — An event which has been found in a stack of nuclear 
emulsion exposed to 105 4.6 GeV x--mesons per cm?, has been inter- 
preted as an example of the decay in flight of an anti-hyperon, Ae ihe 
event and the interpretation are discussed in detail. 


1. — Introduction. 


The concept of particle and anti-particle is now well established and it 
received experimental confirmation a few years ago with the discovery of the 
anti-proton and anti-neutron (12). Up until now however it has not been 
extended to the realm of the hyperons. As distinct from the other group of 
strange particles, the K+ and K~-mesons which are particle and anti-particle, 
the presently known hyperons are expected to have anti-particles which have 
not yet been discovered. They are however expected to exist and will fit 
naturally into the present scheme of elementary particles. The main experi- 


(*) On leave from the University of Bristol. 
(+) Partially supported by the United States Atomic Energy Commission. 

(1) O. CHAMBERLAIN, E. SEGRE, C. WIEGAND and T. YPSILANTIS: Phys. Rev., 100, 
947 (1955). 

(2) B. Cork, G. R. LAMBERTSON, O. PICCIONI and W. A. WENZEL: Phys. Rev. 
104, 1193 (1956). 
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mental difficulty in observing such particles lies in the fact that they must 
be produced in association with a baryon to conserve the baryon number or 
must be produced by an anti-baryon in such a way as to conserve strangeness, 
the threshold energy for production is therefore very high. The only hyperon- 
antihyperon pair whose threshold energy for direct production is within reach 
with accelerators presently working is that of the A°-hyperon. The reaction 
giving an anti A° with protons is expected to be of the form: 


DAS ee 
and with 7--mesons of the form: 
nm +tp>A°+A°+n. 


The threshold energy for production from a free nucleon by a proton is 
7.10 GeV and for a #--meson, 4.73 GeV in the laboratory system; from a 
bound nucleon it can be considerably less if the struck nucleon is moving in 
a favourable direction. If a composite nucleus is used as a target however 
a fraction of A° will be absorbed before they can escape from the parent nucleus. 
When intense beams of x--mesons of ~ 5.0 GeV became available at the Be- 
vatron, it was decided to attempt the detection of an anti-A°hyperon in 
nuclear emulsion despite the numerous difficulties involved, of which the last 
mentioned is one. Because the beam energy is so little above the threshold, 
the particles produced are ejected in a narrow forward cone in the laboratory 
system. It is therefore necessary to search for these particles against a large 
background of beam tracks; it is thus highly desirable to have a pencil beam 
with a sharp cut-off in one dimension at least. The anti-hyperons are expected 
to have a short lifetime, — 10-15, and it is therefore necessary to detect them 
close to the point of production. A convenient way of doing this in nuclear 
emulsion is to use it as a target and as a detector, particles with lifetimes as 
low as 10-15s can then be detected. If this technique is used however no 
separation of the secondary particles is possible either in momentum or mass, 
one is forced to detect all particles produced. There are hopes that very high 
fields which can be pulsed with the accelerator may eventually overcome these 
difficulties. 


2. — Exposure details. 


The stack size was 150 plates 8 in. x7 in. x 600 um G.5 Ilford emulsion. 
The top edge of the stack was exposed to 105 7--mesons/em? which is a density 
which does not allow one to follow through near-minimum tracks although 
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it is possible to search for double stars ete. The density of interactions is 
about 1 per field of view under a magnification of X600. The density of 
tracks falls off to 10%/em? in a vertical distance of 8in. The stack was ex- 
posed with the individual pellicles vertical as the fall off in this direction was 
much greater than in the horizontal direction. The mean beam momentum 
was 4.63 GeV/e at the centre line of the stack (plate 75). The circle of con- 
fusion, of radius 3in; was considerably greater than the dispersion of the 
system, 0.13 GeV/e per inch. The maximum momentum passed by the an- 
alysing system was 5.0 GeV and this momentum would just reach the first 
90 plates in the stack. The beam momentum in the first plates was therefore 
a rather homogeneous mixture between 5.0 and 4.47 GeV/e with a peak at 
4.86 GeV/c, or in terms of energies, between 4.85 and 4.32 GeV with a peak 
at 4.6 GeV. We thus have z~-mesons in the stack whose energies are above 
threshold for A°-production from a free nucleon. 


3. — Experimental procedure. 


The detection of neutral particles in emulsion is always difficult and in 
this case a method first used by FRIEDLANDER et al. (*) to find A°-hyperons, 
that of following back secondary z-mesons, was adopted. The A°-hyperons 
are expected to decay into positive 7-mesons and to antiprotons: 


A°>nt+p+ 37.45 MeV. 


The hyperons will be emitted fast and forward in the laboratory system, 
the x-mesons will therefore be emitted predominantly in the forward cone. 
The following method of scanning was therefore adopted: typical z-u-e decays 
were found by area scanning in a region of the plate where they could have 
more probably been produced by the decay of a A°-hyperon. They had to 
appear to originate in the intensely irradiated part of the stack and they had 
to make an angle of less than 60° with the primary beam direction. 392 7*- 
mesons were followed back until their parent stars were reached or until they 
were lost in the very high background of tracks. Unfortunately about half 
of the x-mesons met the latter fate after a mean range of 2.5 cm. The rest 
were traced to their parent stars and of these 150, 127 appeared to come 
from primary 4.6 GeV 7-meson stars as judged from the direction of a pos- 
sible primary. Three 7-mesons appeared to originate from the decays of 
K+-mesons. The kinematics of one event were consistent with a Ki, decay 


(3) M. W. FRIEDLANDER, D. KEEFE, M. G. K. MENON and M. MERLIN: Phil. Mag., 
45, 533 (1954). 
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in flight and the other events were consistent with the decay of a Ki, at rest, 
giving secondaries 2.01 and 1.5 cm long. The event of interest however is 
a «V » type event at which a z-meson of 32 MeV originated. In the following 
sections it will be shown that this event is the decay of an anti-lambda hyperon 
and probably represents the first example of its kind detected. 


4. — Description of the event. 


The position of the event in relation to the stack is shown in Fig. 1. The 
decay event was found 2.7 cm from the highly irradiated stack edge and 4.2 em 
from the front of the stack. The direction 


STACK EDGE of the A° can be inferred exactly from a 
oe Roe measurement of the direction and energies 


—T beam direction 31° — È 
ol È a_n of the secondaries. The line of flight makes 
43-50 Gev ati i, an angle of 21° with the 7--meson beam 
le E Ae SS direction. The potential path length of 


the A° in the stack is such that assuming 
the normal A° lifetime for the A°, it was 
probably produced in the stack. Although 
an attempt was made to locate the parent star, it was unsuccessful owing to the 
very large background and high star density. A drawing of the event is shown 
in Fig. 2. 

The x*-meson was followed for a distance of 1.70 cm from the x-u decay 
point to the « V » event. The z-u-e event was found 126 um from the glass 
interface and the x-meson traversed 5 plates from its origin in plate 46, 
290 um from the air interface to its decay point. The u-meson range was 
621 um and the decay electron was detected, the 7-u-e event is therefore per- 
fectly normal. The angle of dip of the x-meson was 37° at the « V » event, 
the meson travelling down in the stack. The other track from the « V » tra- 
velled up in the stack at quite a shallow angle of 6°. It traversed 5 plates 
before interacting in flight 2.71 cm from the «V». The interaction is very 
large having 11 prongs, three of which are shower tracks. The projected angle 
in emulsion between the 7 and the p is 51°, making the space angle (64+ Ly 
The ionization of the anti-proton track is about twice minimum at the decay 
point and it changes by 10% over the length of 2.71 em. As this change of 
ionization is very critical for the interpretation of the event, it is treated at 
length in the next section. To prove that the A° is the correct, interpretation 
of this event, it is necessary to show the identity and direction of both decay 
particles and to show that the Q-value is that expected from observations on 
normal A°-hyperons. Even then there is one possibility which cannot be com- 
pletely ruled out, namely that the decay event is the interaction in flight of 


Rive: 
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an anti-neutron producing a single z-meson and an antiproton with an ap- 
parent Q-value in agreement with that for a A°. Such a radical explanation 
is certainly unlikely and we ignore it. 


5. — The anti-proton track. 


As mentioned above the direction of this track is vital for the interpre- 
tation. Unfortunately only 2.7 em of range are available for measurement 
and the change in ionization is only 10% over this distance. Very careful 
measurements are therefore needed to detect the 10% unambiguously. The 
track is completely contained in 5 plates. It starts 290 um from the air inter- 
face in the first plate and ends 64 um from the air interface in the 5-th plate. 
In Table I are listed the corrected ranges in each plate and the portion of the 
total range on which ionization measurements are possible (excluding 50 um 
from the surface and glass interfaces and other regions of obvious distortion 
or bad local development). The effective ranges from the « V » event at which 
each ionization measurement was made are also listed together with the actual 
number of grains counted and the values of the minimum ionization in each 
plate. The value of minimum ionization was determined in the following 
way: the ionization of 4.6 GeV 7-mesons is « plateau » and these beam tracks 
therefore provide a very good means for normalizing each plate. There are 
a very large number present and it is therefore possible to normalize the 
exact portion of the emulsion which is actually required. In each plate 10 beam 
tracks were taken in the same region of the emulsion as the anti-proton track 
and 1000 grains were counted on each. 


TABLE I. 

Total cor- Range N ; Vini . Effective range 
Plate rected range counted SER EUR I from A° decay 
T(46) | 278mm | 2.27 718 15.00 2.104 1.13 mm 
2 (45) 5.88 4.72 1572 15:51 2.150 5.72 
3 (44) | 6.59 5.11 1780 15.62 2.230 11.98 
4 (43) 6.78 5.42 1890 15.39 2.260 18.67 
5 (42) 5.09 4.56 1601 15.19 «| 2.316 24.84 


The top and bottom 50 um were excluded. The results showed a gradient 
of 10% from the top to the bottom of the emulsion. For each plate a mean 
value was deduced by averaging measurements taken at equal intervals in 
depth and this is the value (after conversion to minimum from plateau) which 
appears in Table I. The accuracy of each value is + 1% and the deviations 
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from a common mean are therefore significant. ALEXANDER and JOHNSTON (4) 
give a plateau minimum ratio of 13%. Previous determinations (°°) have 
not given such a high value (about (6--10)%). However as we have adopted 
the calibration curves given by ALEXANDER and JOHNSTON, which are nor- 
malized to «minimum », we have used their value for this ratio. We have 
merely reduced our plateau value by 11.5% to give the minimum values. 
It is realized however that the lack of a precise knowledge of this ratio intro- 
duces a systematic error in the 
! determination of the energy of 
T-MESON x vs Range #K-MESON =the antiproton although it does 
"fa not affect the determination of 
its direction. On the experi- 
mental evidence it is unlikely 
hote that the ratio exceeds 13%, 
but it is possible that it is as 
low as 6%. . We have therefore 
included an asymmetric error 
in the ionization of ‘5%. 

The results obtained are 
shown in Fig. 3. It is clear 
that the track is moving away 
from the decay point and to- 
wards the large star. The lines 
shown in the figure are the 
best fits to the points for various 
particles. It is patent that the track cannot be that of a 7-meson. The points 
are bad a fit to the K-meson line, but probably equally good to the proton or 
hyperon line although the proton fit is slightly favoured. It is not possible to 
distinguish the proton and hyperon masses by these ionization measurements. 
The fit to a hyperon going in the reverse direction is exceedingly bad but if one 
admits the reverse direction, the track is probably a hyperon track. These 
measurements give a value of g* = 2.34 at the point of decay and assuming 
a protonic mass, correspond to a residual range of 10.0 em and an energy of 
230 MeV. The error without the normalizing error is +5 MeV; with, it is 
+7? MeV. 

The track is long, and flat enough to permit a measurement of the mul- 
tiple scattering. This does not give a result of sufficient precision to determine 
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(4) G. ALEXANDER and R. H. W. JOHNSTON, Nuovo Cimento, 5 263 (1957). 

(5) A. H. Morris: Phil. Mag., 43, 533 (1952); Phys. Rev., 91, 423 (1953). 

(9) R. R. DANIEL, J. H. Davies, J. H. MuLvey and D. H. PERKINS: Phil. Mag., 
43, 753 (1952). 
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the direction of the track but a better estimate of the mass can be obtained 
in this way. The value of pf obtained from 30 independent cells of 400 um 
each is (455 +80) MeV/e at the decay point. A proton of g* = 2.34 should 
have a pf of 410 and a hyperon (2), 540 MeV/c. The mass is thus 
(1020 + 150) MeV. Again it is not really possible to decide between a hyperon 
and proton mass. 


6. — The anti-proton star. 


Having proved that the particle from the decay event is entering the inter- 
action and not leaving, it is easy to show that the track cannot be caused by 
any known particle other than an anti-proton or anti-hyperon by examining 
the visible energy in the star. An attempt has been made to identify all the 
particles coming from the star. A list of the identities and the ranges is shown 
in Table II. The identity of the shower particles is especially important and 


TABLE II. 
; ‘ | 
ong Identity Range Ionisation Scattering measurements pi 
1 | Proton or | 256 um — — 6 
x-particle 
2 | Proton 1.56 cm = “= 67 
3 | n*-meson — g*=1.31- min at star pfatstar=(210+30)MeV/e 87 
4 | Proton (+) | 96 um se — 3 
5 | Proton (+) 1576 um —- — 10 
6 | Proton (+) |156 um — = 5 
7 | x--meson — g* =1.21-min at star, pB at star=190 MeV/c 108 
it interacts after |pf at intern.— 70 MeV/c 
19.8 em when g*=2.7 

8 | Proton (+) | 56 um — — 2 
9 | x-particle |106 yum — — 15 
10 | x*-meson 3.5 em \g*=1.70 at star pp at star= 86 MeV/c 50 
11 | «-particle 46 um — — 9 
Total visible energy 362 
Binding energy 70 
Neutron estimate 93 
Neutron B. E. 50 
nm -meson rest mass 420 
(+) We cannot exclude the possibility that these r°-meson K. E. estimate 122 
tracks may be due to deuterons. However, their energies x°-meson rest mass 210 

have been determined on the assumption that they 
are eae Total energy released | 1327 
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these tracks have had multiple scattering measurements made where possible 
together with ionization measurements. In one case the track was arrested 
in the emulsion and it produced a typical 7-1 decay. event. 


Track No. 2 has been arrested after 1.56 cm and is a proton of 67 MeV. 


Track No. 3 has an ionization of 1.31 xmin (+ .05).. This track could 
not be traced through after the first six pellicles and so no information on 
the range is available. Scattering measurements however give a pf of 
(210 + 30) MeV/e from 40 independent cells of 100 um. If the track were 
a proton it should have a pf of 1800, a z-meson should have 140 MeV/c. It 
has therefore been taken as a 7-meson of 87 MeV. 


Track No. 7 has an ionization of (1.21+ 0.03) x min and it interacts after 
9.8 cm when its grain density is 2.7 + 0.1. Thus by Ag*/AR it is a 7-meson. 
Confirmation from scattering gives pB =190 +40 from 25 independent cells 
of 100 um at the anti-proton star and (70 + 15) MeV/c from 25 cells at the 
interaction. As the energy at the point of interaction is only 16 MeV, it is 
probably of negative charge. The secondary prongs from the interaction are 
both protons of 55 and 5 MeV. 


Track No. 10 has been arrested after 3.51 cm. It decays into a normal 
u-meson and is therefore a 7-meson of 50 MeV. Scattering confirms this, giving 
a pp of 86 MeV/e with an ionization of g* =1.7 x minimum. 


The total visible energy as can be seen from the table is 363 MeV + 420 MeV 
in the rest masses of the shower tracks. The total is thus 783 MeV. .If we 
include 8 MeV binding energy per nucleon and include 4 the total x*-meson 
energy for neutral 7°-mesons and an equal amount of evaporation energy for 
neutrons as for the protons, one obtains a total energy of 1327 MeV—a 
little lower than that available for an anti-proton annihilation — 2108 MeV. 

With a total visible energy as high as 782 MeV all interactions involving 
hyperons, K-mesons or x-mesons are ruled out as explanations and the only 
tenable hypothesis is that it was caused by the annihilation of an anti-proton 
or anti-hyperon. The latter is virtually impossible, because of insufficient 
energy in the beam to create anti-Z hyperons. They can of course be produced 
by a secondary interaction of anti-proton: p+ p> dt+ K-+n. 


7. — The Q-value in the decay. 


The best value of the mass of the A° has recently been obtained by the 
Barkas group (7), it is 1115.3 MeV and leads to a Q-value of 37.45 MeV. 
From the opening angle of the «V» and the energies of the secondary par- 


(7) W. H. BARKAS: Report to the Padua-Venice Conference, (Sept. 1957). 
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ticles, the best Q-value obtained for this event is (35%) MeV. The main con- 
tribution to the error comes from the energy of the anti-proton from the ioni- 
zation measurement. Reducing the plateau minimum ratio in nuclear emul- 
sion would have the effect of increasing the Q-value. The dependence of the 
Q-value on 0, the opening angle, and on E. is shown in Figs. 4 and 5. The 


-A 
Q-VALUE 
(MeV) 


50 See 


p — 
100 200 300 350 


error in the opening angle + 1° contributes + 0.8 MeV to the error in the 
Q-value, while the statistical error in E_ contributes + 0.7 MeV. The Barkas 
range energy relation (5) has been used and the density of the emulsions at 
the time of exposure was determined by obtaining the mean range of flat 
u-mesons from 7-meson decays found during the course of the experiment. 


8. — Conclusion and discussion. 


The ionization measurements on the anti-proton track, the visible energy 
in the interaction and the apparent Q-value all point to the conclusion that 
the event represents the decay of an anti-lambda hyperon. The only tenable 
alternative, if one accepts the anti-proton direction, is that the decay event 
is an interaction produced by an anti-neutron. We can never exclude this 
possibility but it is surely very unlikely. 


(8) W. H. BARKA8, P. H. Barrett, P. Curr, H. H. HECKMAN, F. M. Suiru and 
H. K. Ticno: Nuovo Cimento, 8, 185 (1958). 
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If one does not accept the direction of the anti-proton track, one can postu- 
late that the large interaction is a neutral star which gives rise to à track 
which either decays or interacts to produce a single 7-meson of 32 MeV. The 
ionization measurements at the decay point and at the interaction enable 
one to exclude a number of decay possibilities: 


a) K},-decay in flight. The range of the z-meson defines the momentum 
of the K-meson at the decay point as 370 MeV/c. The grain density of the 
track should therefore be 2.1. This'is within the error of the measured value but 
the grain density at the interaction is 8 standard deviations from that expected. 
The pp measurement would also have to be wrong. 


b) X+-hyperon decay in flight. Again the momentum ot the hyperon 
is defined and the grain density should be 2.56, 8 standard deviations from 
the measured value. 


c) A positive cascade hyperon decay in flight. The momentum is de- 
fined as 460 MeV/e which would result in a grain density of 1.80, 6 standard 
deviations away from the measured value. 


The interaction of a positive or negative X-hyperon to give a single posi- 
tive z-meson is also very unlikely. Hence even if the direction of the anti- 
proton is not believed it is a very difficult event to explain away by an expla- 
nation involving known particles. 


We are greatly indebted to Dr. E. J. LOFGREN and his colleagues at the 
Berkeley Bevatron for making the exposure possible, and to Professor D. H. 
WILKINSON and Dr. WARREN CHUPP who carried out the exposure. Dr. ©. 
WALLER of Ilford Inc. is again to be thanked for the careful preparation of 
the stack. We thank Professor N. DALLAPORTA and Professor C. F. POWELL, 
Dr. D. Evans, Dr. H. Huzira, Dr. S. NATALI, Professor GERSON GOLDHABER, 
Dr. SULAMITH GOLDHABER, Professor D. H. SrorK and Professor H. K. TIcHO 
for many fruitful discussions. 

We appreciate the encouragement of Professor N. DALLAPORTA, Professor 
M. MERLIN, Professor C. F. POWELL and Professor J. R. RICHARDSON. D. J. P. 
thanks Professor C. F. PowELL for the loan of plates, Professor N. DALLA- 
PORTA and the University of Padua for hospitality in the fall of 1957 and 


Professor J. R. RICHARDSON for the kind hospitality of his laboratory during 
1958. 
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As an example of what we can expect kinematically we show in Fig. 6 
two curves which give the maximum angle of emission in the laboratory 
system plotted against the energy of the A°. In curve 1 we have assumed the 
momentum of the incoming +-mesons to be 5 GeV/c, the proton to be moving 


= Are Pa 
1 {2,25 GeV/c ROZ CMS. gp 


À 7 P. 
2 (P.2475GeVe = Bp70.16 CMS. À 
P 
p 
P 


a velocity of f = .23, and that 


in the opposite direction to the x-meson with 
opposite to both the mo- 


the momentum of the A® hyperon in the c.m.8. is 


menta of the A° and the proton. 
All values of 01» and Zx within the curve are permi 


example for an intermediate case. 


tted; curve 2 gives an 


RIASSUNTO 


leari esposte a 105 7- di 4.6 GeV/em? 


Un evento, trovato in un pacco di emulsioni nuc 
o di un anti-iperone, A°. 


è stato interpretato come un esempio del decadimento in vol 
L'evento e l’interpretazione sono discussi in dettaglio. 
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A Search for Anisotropy of Inertia. 


G. Cocconi and E. SALPETER 


Cornell University - Ithaca, New York 


(ricevuto il 6 Agosto 1958) 


Summary. — One of the possible consequences of the Mach principle is 
that the asymmetric position of the Galaxy with respect to the Earth 
induces an asymmetry in inertia. A test of this possibility is proposed. 
It consists in the microwave measurement of the Zeeman splitting of 
electron levels having different orientations relative to the direction 
Earth-Center of the Galaxy. The isotropy of inertia on the Earth can 
thus be checked with a precision of about 1 part in 101, 


One of the fundamental assumptions of mechanics and gravitational theory 
is that the equality of the inertial and gravitational masses of a body is exact, 
and independent of the position and motion of the body. This assumption 
is an extrapolation to infinite accuracy of observations which obviously have 
a finite precision. It is also an observational fact (of finite precision) that the 
so-called «fixed stars » do in fact represent an inertial reference frame. 

This last observation has led to the formulation of Mach’s principle, which 
asserts that the inertia of a body at any point in the Universe is determined 
by «the fixed stars », which must be understood as the total distribution of 
matter in the Universe. 

This distribution, however, is uniform and isotropic only to a certain degree 
of approximation. If Mach’s principle holds, we might then expect that the 
slight asymmetries in the distribution of matter at large would result in slight 
deviations from at least some of the laws of mechanics and gravitation which 
are commonly assumed to be exact. It is thus desirable to find tests that can 
demonstrate the validity of Mach’s principle and of its consequences. In this 
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article we shall concentrate on the question of the isotropy of inertia and 
describe a relevant experiment more sensitive than those performed thus far. 

To good accuracy, the acceleration produced on any body by a force points 
in the same direction as the force, and the ratio of force to acceleration, the 
inertial mass, is independent of the direction of the force; é.e., inertia is iso- 
tropic, and inertial mass is a scalar quantity. However, if Mach’s principle 
holds, these statement should be true only in first approximation and asym- 
metries in the matter distribution at large, say concentration of matter near 
the center of our Galaxy, could produce asymmetries in inerta. 

Mach’s principle alone does not specify the nature of the effect that matter 
at large has on the inertia of a body; however, a rather plausible possibility 
for the contribution to the inertia of a test-body, resulting from a mass M 
a vector distance r away, is the following: 


a) The contribution is proportional to M and to a negative power, 
— y, of the distance r. 


b) The contribution depends on the angle 9 between r and the direction 
of acceleration of the test-body (1). 


The most attractive possibility for this 6-dependence, which still preserves 
energy conservation in the usual sense is that the contribution to the inertia 
of the test-body has a maximum value for cos @|—1 (motion parallel or 
antiparallel to r) and is zero (?) for cos 0=0 (motion in a plane perpendicular 
to r (°)). 

Jith these assumptions, the inertial mass of a test body becomes a sym- 
metric tensor 977,; for any distribution of matter in the Universe. E.g., Newton’s 


(1) One cannot, of course, exclude the other possibility that the contribution to 
inertia is independent of 0. In this case inertia would still be isotropic for any distri- 
bution of matter in the Universe, although the numerical value of the inertial mass 
of a test-body would depend on its position in the Universe. 

(2) The more general case of a non-zero minimum for cos 9 = 0 gives qualitatively 
similar results, except that the expressions for AQT, to be discussed later, are 
decreased by a constant factor. 

(3) The assumption that inertia might be larger when the test body is, for instance, 
accelerated toward M and smaller when accelerated in the opposite direction, has the 
rather drastic consequence that ordinary forces would no longer be conservative in 
the usual sense. It would follow that, e.9., all atoms, molecules and, consequently, 
also bulk matter would be electrically polarized with the electric dipole moments pointing 
toward the Galactic Center; another consequence would be that any compound system 
with internal forces would be repelled by matter by amounts proportional to the 
strength of the internal forces. It would presumably be too fanciful to seek any con- 
nection. between this effect and the observed expansion of the Universe. 
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law is generalized to 


3 
F; — > Misa; , 


j=t 


where F and a are force and acceleration, respectively, and the kinetic energy 
T becomes $ > (977-);;p;P;, where p is the linear momentum. The mass tensor 


2 
is AEREA by the directions of the three principal axes and the corres- 
ponding diagonal elements 977;,. 

We shall be mainly interested in the effect of the mass concentration near 
the center of our own Galaxy. In this case; neglecting the mass distribution 
in the Galactic Disc, the mass tensor of a test-body has the direction to the 
Galactic Center as one principal axis with diagonal element 97 + A977, the 
diagonal elements for the other two axes being 97 —3A977 each. 

The order of magnitude of the expected anisotropy can be evaluated as 
follows: If o is the average density of matter in the Universe, the isotropic 
part 977 of the inertia of a body on the Earth is proportional to 


R 


, 470 450 
Inertial mass = 977 «| = € dr = —— R*” , 
r 3 —y 


ü 


where PR is the maximum distance from which matter still influences the inertia 
on the Earth. The anisotropic contribution A977 to the local inertia from a 
body of mass M at a distance 7 away is thus 


AQT a eu d 
GTA 4no R®-” : 


(1) 


To give a numerical example, assume that 9 = 105$ g em, and that R=cT = 


=10?8 em, the so-called radius of the Universe deduced from the Hubble 


constant. Then for »=1 (*) the extent of the anisotropy due to the matter 


near the center of the Galaxy would be of the order of 1 part in 107 and for 
y= 0.25 1 part in 10". 

It must be emphasized here that these numbers depend strongly on the 
values assumed for the density and the radius of the Universe, quantities 
which are not yet known with any accuracy. (The corresponding quantities 
for the Galaxy are instead rather sound.) The value of R is particularly un- 
certain since the value of the Hubble constant is derived only by extrapolation 


(4) Values of » larger than 1 can be excluded a priori since they would make the 
masses of the Sun and of the Earth prominent. 
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of the red shift from relatively small distances and the detailed relation between 
optical red shift and decrease of inertial influence of matter is not known. 
In fact, it is likely that our assumed value of R is an underestimate. 

In principle, this anisotropy would be observable in classical macroscopic 
experiments; e.g., there should be a diurnal variation in the period of a quartz 
crystal (or a pendulum) clock if one assumes that the anisotropy of electric 
forces (or gravitational mass) is either absent or at least different from that 
of the inertial mass. However, present-day, observations place an upper limit 
of about 1 part in 10° on a possible anisotropy observable with a quartz clock. 

We want to propose a more sensitive test that utilizes microwave measu- 
rements of the splitting of atomic energy levels. In principle it consists in 
the comparison of the periodic motion ot two atomic electrons, one moving 
in a line pointing towards our Galactic Center, the other moving in a plane 
perpendicular to this direction. The comparison is realized by measuring the 
energy difference of the two levels. 

Consider, for example, a single-electron atom in a static magnetic field 
with orbital quantum number / and magnetic quantum number m, referred 
to the field direction as axis (component of orbital angular momentum parallel 
to this axis). For m,= 0 (or small) the electron moves preponderately paral- 
lel to this axis; for |m,| equal (or close to) / there is preponderance of mo- 
tion in a plane perpendicular to this axis. 

What then must be measured is either the fine structure splitting between 
two single-electron levels with the same value of m and differing values of j, 
or the Zeeman splitting between two levels with the same value of j and dif- 
ferent values of m. What one needs is only the variation of such a fine structure 
or Zeeman splitting as the orientation of the magnetic field (which determines 
the m-axis) relative to the Galactic Center is changed. This should be ob- 
served either by rotating the experimental apparatus or merely by looking 
for sidereal variations (as the orientation of the laboratory co-ordinates relative 
to the Galactic Center is modified by the Earth’s rotation). 

For a given value of AQ77/977, estimated in eq. (1), the energy shift due to 
any anisotropy of inertia is given by 


AT — 
(2) AR =", TP, 


where T is the average kinetic energy of the single (or valence) electron and Pi 
is the expectation value of the Legendre polynomial of order 2 with the di- 
rection to the Galactic Center as polar axis. This perturbation is equivalent 
to that of a weak quadrupole field and P, changes sign when the magnetic 
field direction (with quantum numbers fixed) changes from parallel to per- 
pendicular to the direction of the Galactic Center. 


42 © Il Nuoro Cimento. 
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We quote a few examples for P, with the magnetic field pointing to the 
Galactic Center and for a single p-electron. Let 1, 9, I and F be the orbital, 
total electronic, nuclear spin and total angular momentum (electronic plus 
nuclear spin) quantum numbers; mu, #;, M, and m, the corresponding mag- 
netie quantum numbers. In the absence of hyperfine structure (1 = 0) Ps is 


zero for j= 4; P.=— for m=+$, P= +5. for mrs +4 if j= 3 (p3- 
electron). For non-zero nuclear spin J and very weak magnetic field (Zeeman 
effect < hyperfine structure) and F=I4+4+j we have P,=—} for m,= F and 


P,= (j —D/BG+D) for m,=F—1. 

Precision experiments of the kind described above have of course been 
performed, although they were not designed for investigating anisotropy of 
inertia. One example is the microwave experiments of LAMB and collabo- 
rators on fine structure splitting in hydrogen and deuterium (5), which measure 
the energy differences between 2p, and 2p;-states (about 104 MHz while the 
average energy corresponds to about 3-10° MHz). 

Still more sensitive tests for anisotropy of inertia are provided by precision 
microwave measurements of Zeeman transitions in very weak magnetic fields 
for atoms with ground states with j> 3. Relevant experiments which have 
already been carried out are those of KUSCH and collaborators (*) on Zeeman 
transitions in the p,-states of Ga and In. 

However, since the measurements conducted thus far were not intended 
for observing sidereal variations, the various hours in the sidereal day where 
not covered in a systematic manner, and the orientations of the apparatus 
(relative to the polar axis) were not chosen to give the maximum sidereal 
variation (7). Nevertheless sidereal variations in the Zeeman splitting energies 
as large as 0.2 MHz, say, would probably have been noticed by the exper- 
imenters and the negative results so far place an upper limit of about 107° 
on AQIT/W7. 

Further improvement in sensitivity could easily be achieved by similar 
experiments repeated with the purpose of studying sidereal variations or de- 
pendence on the orientation of the apparatus relative to the Galactic Center. 
Very weak magnetic fields giving Zeeman splitting energies less than 100 MHz, 
say, transition between F=I+j, m,=F and F—I+j, m,=#—1 in atoms 
with relatively small nuclear spins and p,-ground states would be particularly 
favorable. In this manner sidereal variations of the order of one kHz might 


(3) E. S. DAYHOFF, 8. TRIEBWASSER and W. E. Lams: Phys. Rev., 89, 106 (1953). 

(6) P. KuscH and H. M. Forex: Phys. Rev., 74, 250 (1948); A. K. Mann and 
P. KuscH: Phys. Rev., 77, 435 (1950). 

(7) At latitudes near 45° N the most favorable direction is the NS direction in a 
hôrizontal plane. This line points towards the center of our Galaxy (due 8) at sidereal 
time about 17.30 h; twelve hours later the Galactic Center is at the nadir. 
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be detectable and one might hope to reach accuracies for A977 /Q of the order 
of one part in 10" or more. With such accuracy any anisotropy due to the 
presence of our Galaxy would be detectable if the parameter v in Cd (1) 238 


around unity, even if the value of R is substantially larger than that used 
in the example discussed earlier. 


The authors are grateful to Professor W. Lams for helpful suggestions; 
they also wish to thank the Gravity Research Foundation, New Boston, New 
Hampshire, for the interest shown in the ideas expressed in this article. 


RIASSUNTO (*) 


Una delle possibili conseguenze del principio di Mach è che la posizione asimme- 
trica della Galassia rispetto alla Terra induca un’asimmetria nell’inerzia. Si propone 
una verifica di tale possibilità, consistente nella misura delle microonde della separazione 
per effetto Zeeman dei livelli elettronici aventi differenti orientazioni rispetto alla 
direzione Terra-centro della Galassia. È così possibile verificare l’isotropia dell’inerzia 
sulla terra con una precisione di circa 1 a 10". 


(*) Traduzione a cura della Redazione, 
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Decays of Hyperfragments in Nuclear Emulsion. 


C. GROTE 


Deutsche Akademie der Wissenschaften zu Berlin 
Kernphysikalisches Institut - Zeuthen bei Berlin 


(ricevuto 18 Agosto. 1958) 


Summary. — Seven decays of hyperfragments have been found in Ilford G-5 
emulsion. Five events are uniquely defined decays. Within the limits of 
error the binding energies agree with the results of other authors. 


170 double stars in nuclear emulsion were selected for decays of hyper- 
fragments. We used the criteria described by FILIPKOWSKT, GIERULA and 
ZIELINSKI in (1) for selection. Four decays have been found in plates which 
are part of the I-stack, three decays in plates exposed by the K~-beam of 
the Berkeley bevatron. È 


TABLE I. 
Primary star Secondary star 
| Time of flight 
Event PSE Kinetic energy ani of the 
no. Produced Type of all charged RER eae fragment 
by particles (MeV/c) (10-12 s) 
(MeV) 
| 
1 | Cosm. rays 6+2n 34.6 = 7.5 (if 5Hex) 
2 » 13-+1n 33.1 30 4.7 
3 » K -capt 35.1 10 1.1 
4 ue oes 17+0n 68.9 400 2.5 
5 K -beam K -capt 35.9 6 8 
6 » » 65.8 405 1.2 
7 » ) 38.1 9 2 (if *H,) 


(1) A. Frurexowskt, J. GreruLa and P. ZIELINSKI: Acta Phys. Polon., 16, 139 (1957). 
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TABLE II. — Range R, polar angle y, dip angle d, kinetic energy T and relative grain 
density g/qo of the secondary star tracks. 


(*) g is the angle between the track 
in the plane of the emulsion. # is the angle between the track 


emulsion. 
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Bvent Track KR œ (*) d (*) vil 
PE (um) (deg.) (deg.) 9190 ee (MeV) 
F 195 + 6 |349 22146643 45,.7He,|18 +1 (if Hea) 
1 1 IAE 6.6 | 121.5-¢1 |} — 25-43 1H LO 
2 = 2 sy 120 | w+45 3.45He | 0.63-+0.05 
3217515 +50 |301.54+1|—16+3|:2:0.4+0.15] n-|32.25+0.7 
F CSM 9 SEL) 2743 >iFlie (1200-05 
CATE RO 14.34.72 32 +6| ~+40 4He | 3.92+0.52 
2 dr CRUE Li AI 1H 6.56 +0.2 
3.19580 +100 |216 +1|—26+3| 2.5 40.25 | n- | 22.61 L0.54 
F 114+ 0.51316 +5|—30+5 5He,| 3.6 +0.2 
ey ARRE 10.4+ 0.5|118.5+3 0+5 4He | 3.1 0.28 
2 |10490 + 50 | 247.5+1|— 23+2| 2.55+0.1 n- | 23.85+0.54 
2. 400 CET 1327541) 4202 1H 8.18+0.11 
F CY Vga OO | 66 +2 0+5 Lig (13°. 40.5 
ACTA > 8500 8541-2300) 44 10.451 1H, | 59.2 22.8 
2 35.54 :6.5| 224 £2) +20415 3H 2.6 +£0.34 
3 BOG tk Bul 303, 41444712 2H 6.97 £0.15 
F 76. +01 |350 +3 | +1745 4H, | 4.8 +0.2 
Al 85.5+ 2 |359 +2 | +3943 3H 48 01 
2 24 + 4 |143.541|—37£42 1H 5.62+0.1 
3 111760 +140 |214 +2|—31+3|2.1 +0.15 n- |25.50+0.61 
F RSH iss Say ou | 104 +2|—28+5 TL tee sew 
6 | 1 | 9635 +170 26 +1) +42+2| 45 +0.2 1 91.10.75 
2 53.5+ 1 dg — 35-46 3H 3.47 £0.08 
3 493 + 6° 273 41) +4342 2H 11.16 £0.16 
F 320 + 4 |319.541| +10+3 23,44, | 11.3 +0.3 (if HA) 
tg | 982 + 3 |140.5+1|—25+2 1H 6.63 +0.1 
2 424 0.5| w 255 | ~ —10 1,2,3H 0.45+0.1 
3.116360 +120 |342.5+1 | +31+2| 1.8 +0.1 x- | 31.00+0.71 
| 
i 


and an arbitrarily chosen direction, both projected 
and the plane of the unshrinked 
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For the identification of singly charged particles we have used the range- 
energy-relation of BARKAS (?), for multiply charged particles the curves of 
VIGNERON (*) and PAPINEAU (*). In computing binding energies the tables of 
SEGRE (5) (isotopic masses) and ROSENFELD (5) (masses of elementary particles) 
have been employed. The following errors have been taken into account: 
a) Uneertainty of range-energy-relation (+ 0.60%, 8. (7)).b) Error of range 
measurement; €) Variation of stopping power (+ 0. 5%, 8. (5)); d) Hetero- 
geneity of emulsion (s. (°)); e) Range strag ggling (s. (*)); f) Uncertainty in 
Qi (8. (°)). 

In all events reported here we have assumed that the hyperfr agment decayed 
after coming to rest in the emulsion. The lower limit of the life time of the 
hyperfragments following from their time of flight is then 10° s. 

The interesting data of the events have been summarized in the Tables I 
and II. A drawing of each event is shown in Figs. 1-7 (A, primary star; 
B, secondary star; F, interconnecting track). 


Event No. 1. — The tracks of the secondary star are coplanar within the 
limits of error. The interconnecting 
track F has no visible d-rays or gaps. 
A few x before producing star B the 
fragment undergoes large-angle Coulomb 
scattering. Therefore it is supposed 
that the decay occurred at rest. Par- 
ticle 3 is a 7--meson producing à o-Star 
at the end of its range. Track 1 has 
four visible gaps suggesting a charge 
z<2. Assuming that particle 1 is a 
proton and particle 2 a helium nu- 
cleus, the total momentum will be 
zero. The decay scheme will therefore 
appear as follows: 


Bic. als 


Pee > 1H + 3.4.6He + x + (A) A 


Event No. 2. — The tracks of the secondary star are coplanar within the 
limits of error. Track F has no visible Ô-rays and undergoes wide angled 


2 


(3) W. H. BARKAS: private communication (March 1957). 
(3) L. VIGNERON: Journ. Phys. et Rad., 14, 145 (1953). 
(4) A. PAPINEAU: Rapport C.E.A. no. 543 (1956). 

(5) E. SEGRÈ: Experimental Nuclear Physics, vol. 1 (1953). 
(9) 

RE 

CR 

ay 


3 


is "i H. RosENFELD: private communication (April 1957). 


7) Genoa-Milan Report no. 5. 


Levi SETTI, W. E. Starter and V. L. TELEGDI: preprint April 1957. 
. H. Barxas, F. M. Suirx and W. BIRNBAUM: Phys. Rev., 98, 605 (1955). 


8 
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Coulomb scattering before producing star B. The end of track 3 seems to 
be a typical ox-event. The interpretation of this track as z~-meson has been 
supported by grain counting. Track 1 and 2 have 
been identified by momentum balance. The ioniz- 
ation of both tracks suggests a charge of 2 <2. If 
track 1 is interpreted as an «-particle and track 2 
as a proton, the residual momentum is 30 MeV/c 
corresponding to a neutron of 0.5 MeV energy. If 
we consider the emission of such a low-energetic 
neutron as improbable and the residual momentum 
as due to the error of measurement, we can only take 
one decay scheme into account: 


‘He, > ‘He + ‘1H +7 +0. 


For the binding energy B, we have found 


Fig. 2. 


B, = (3.8 £0.8) MeV. 


Event No. 3. — The tracks of the secondary star are coplanar within 
+4 degrees. Particle 2 is a 7--meson producing a o-star at the end of its 
range. Track 3 could be identified as 
a singly charged particle by measure- 
ment of the total gap length. If it is 
interpreted as a proton and particle 1 
as an «-particle, the residual momentum 
is only 10 MeV/c. The decay process 
therefore must be as follows: 


‘He, > ‘He +x-+1H4+Q; 


B, = (1.8 + 0.6) MeV. 


Event No. 4. — The tracks of the 
secondary star are not coplanar. A re- 
markable thinning-down of the frag- 
ment track suggests a charge of 2 > 3. 
Track 3 leaves the emulsoin. The par- 
ticle has been identified by measuring the grain density at different places 
of the track. The ionization of track 2 and 3 suggests a charge of 2 <2. 
Assuming this we computed all possible decay schemes. When the emission 


Fig. 3. 


3617 


656 C. GROTE 


of more than one neutron is excluded, only one reaction is in agreement with 
the hyperfragment interpretation, namely: 


"Li, > 7H 4° ns 0); B, = (3.7 + 5.2) MeV. 


ea As 
"a 
/ 
Fig. 4. Fig. 5. 
Event No. 5. — The tracks of the secondary star are coplanar within 


+2 degrees. The fragment undergoes wide angled Coulomb scattering a few 
u before producing star B. The end of track 3 seems to be a oz-event. This 
assumption has been supported by measurement of grain density. Assuming 
that particle 1 is a triton and particle 2 a proton, the total momentum will 
be zero. This interpretation agrees with the appearance of the tracks. Only 
one decay scheme can therefore be taken into account: 


SH, H+'H+r +0; B= 41.0, 490.7) Me: 


4 Event No. 6. — The tracks of the 
i secondary star are not coplanar. 
From an energetic analysis of the pri- 
mary star follows that, beyond reaso- 
nable doubt, the secondary star has 
been produced neither by a collision 
of a stable particle nor by the cap- 


É : ture of a negative meson or hype- 
us ron. Track 1 has been identified as 
f proton track by grain density. The 

Fig. 6. ionization of tracks 2 and 3 suggests 
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a charge of 2 <2. Assuming this we considered all possible decay schemes 
with the emission of only one neutron Only one case agrees with the hyper- 


fragment interpretation, namely: 


Li, > 4H +*H4+*H +n+Q; 


the binding energy is 
B, = (5.141.6) MeV. 
The decay is uniquely defined, if the emission of two or more neutrons is ex- 


cluded. 


Event No. 7. — The tracks of the secondary star are coplanar within the 
limits of error. Tracks F and 1 have been identified as singly charged par- 
ticles by measurement of the total gap length. Particle 3 is a x--meson pro- 


ducing a o-star at the end of its range. From momentum balance follows 
that track 3 corresponds to a hydrogen nucleus and track 1 to a proton. The 
decay process must therefore be as follows: 


234 > 1H + 4H +n +0 ; Q ~ 38 MeV. 
For the binding energy we find 
B,=(—1.1+0.8) MeV. 


* o >*k E 


I am greatly indebted to Dr. K. LANIUS for his encouragement and for 
the opportunity of many critical discussions. I further wish to thank Mr. D. 
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BeBeL, Mr. H. W. MEIER and the teams of the emulsion laboratories of Buda- 
pest, Praha and Warszawa for the many measurements on tracks. I would 
like to thank the assistants of our laboratory for their efficient help in scanning. 


RIASSUNTO (*) 


i Si sono trovati sette decadimenti di iperframmenti in emulsione Ilford G-5. Cinque 


eventi sono decadimenti univocamente accertati. Nei limiti d’errore, le energie di 


legame si accordano coi risultati di altri autori. ~. 


(*) Traduzione a cura della Redazione. 
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An Investigation on Plasmas in External Magnetic Fields (*). 


II. - Varying Fields. 


G. SCHMIDT 


Israel Institute of Technology - Haifa (+) 


(ricevuto 18 Agosto 1958) 


Summary. — The behaviour of a plasma placed in a varying external 
magnetic field is investigated considering two limiting cases: 1) the ideal 
hydromagnetic (co); 2) the collisionless one. Restricting the cal- 
culations to general cylindrical geometries (0/02=0) and quasi-stationary 
processes, a complete set of plasma equations are obtained for both cases. 
These include among others équations of state also for the tenuous 
plasma. There is a close resemblance between the two sets of equations 
referring to completely different physical conditions. The equations can 
be readily solved to express every macroscopic plasma quantity in terms 
of the external magnetic field strength. The extension of our consi- 
derations to non-quasi-stationary processes is briefly depicted. It leads 
to a phenomenon which suggests new ways for the increase of plasma 
temperatures. 


1. — Introduction. 


In paper I (:) the stationary behaviour of a quasi-neutral plasma placed in 
an external magnetic field was investigated. With the exception of the deri- 


(*) The present work was essentially finished by the end of 1957 and a full account 
of the basic ideas was given on the Ninth Meeting of the Israel Physical Society, 
April 2, 1958. (Bulletin of the Research Council of Israel, Section F, vol. 7F, p. 99). 
Due to the sharp expansion of literature on the subject several papers as refs. (?*) 
have appeared during the last year covering some parts of the present work. 

(+) Present address: Stevens Institute of Technology Hoboken, New Jersey. 

(1) G. Scumipt: Nuovo Cimento, 10, 55 (1958). 

(2) E. N. PARKER: Phys. Rev., 107, 924 (1957). 

(8) Ia. P. TERLETSKIJ: Journ. Exp. Theor. Phys., 5, 755 (1957). 

(4) Nucleonics, 16, no. 5, 122 (1958). 
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vation of the pressure balance equation—which has been performed in full 
generality—certain restricting assumptions were made: 


1) A tenuous plasma was considered neglecting collisions among particles. 


2) Only «2-dimensional geometries » were treated (straight and parallel 
B lines). 


In the present paper plasmas in time dependent fields are investigated. 
It is generally assumed here that the variation of the field is slow, therefore 
the process can be treated in a quasi-stationary way, 1.6. the system is moving 
through equilibrium states as described in paper I. The condition of 2-dimen- 
sional geometries is also retained. 

The physical problem to be investigated is the following: A plasma of 
finite dimensions in the x-y plane, is placed in an initially homogeneous magnetic 
field B (0,0, B.). Owing to plasma diamagnetism this field becomes distorted, 
so that in the steady state the pressure balance equation is satisfied. If now 
the external field is gradudally changed, the plasma characteristics must also 
change. It is shown that increasing the field, the plasma becomes gradually 
compressed and its energy (or kinetic temperature) raised. On the other hand, 
on decrease of the field, plasma expansion and cooling occur. Owing to 
the quasi-stationary condition these two processes take place in a symme- 
trical way. 

Two opposite limiting cases are considered. The first one is approximately 
valid for very dense plasmas, with frequent particle-particle interaction. In 
this case the plasma behaves similar to a fluid with very high conductivity, 
and therefore the magneto-hydrodynamic equations can be applied. For the 
sake of simplicity the conductivity is assumed to be infinite. 

The second case is that of tenuous plasmas as considered in paper I. In 
this article, however, a further restriction is also made; the particle motion 
is treated in the framework of the usual «first order theory ». It is hoped 
that a more general treatment on the lines of paper I will be forthcoming 
later. 

The calculation of these two extreme cases leads to very similar results. 
Validity of the assumptions made, and order of magnitudes of the parameters 
for realizing these processes under laboratory circumstances are estimated in 
Appendix A. 

Appendix B deals with considerations for non-quasi-stationary systems, 
where an interesting new phenomenon occurs. Calculations for a simplified 
model show that compression and expansion are no longer symmetrical pro- 


cesses, and a periodic operation leads to continually growing kinetic tempe- 
ratures. 


3622 


he "a 


AN INVESTIGATION ON PLASMAS IN EXTERNAL MAGNETIC FIELDS - II 661 


2. — Magneto-hydrodynamic approximation. 


We investigate the case of a fluid, with very high electrical conductivity o 
placed in a varying magnetic field. The variation of the magnetic field inside 
the fluid is governed by the familiar equation expressing the «constancy of 
magnetic flux». This can be obtained as follows: 

The current density i, in a fluid element moving with the velocity v, in 
the local magnetic field B, and electric field E, equals: 


(1) i=o(E+vxB). 


Taking the curl of eq. (1) and making use of Maxwell’s equations curl E = 
= — 0B/ct, we obtain: 


A 


B 
{2) curl: — 0 (- = + curl v <B) : 


With the help of the equation curl H = j where the displacement current is 
neglected, as o is high and ©/ct small, it follows: 


(3) - VB = | i + curl vx B) : 
py 7 


where use has been made of div B= 0. Integration of eq. (3) over a «ma- 
terial surface » S, whose boundary s is fixed to the moving fluid elements, 


yields: 
(4) Sa V*BdS (cn (vx B)| dS. 
uo ct 
Ss Ss 
The second term on the right hand side equals: 
(5) — prox Bas =p Berxds). 


This is the time rate of flux change within the material surface owing to 
the motion of the boundary, while the first term on the right hand side of 


eq. (4) signifies the change of flux caused by the time variation of the field. 
The total time rate of change of the flux within the moving boundary is 


finally : 


Oh See ei RR 
at Bas, 
ve de x 
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If o > co this reduces to: 


d@ 
(7) ate 0 


The connection between the magnetic field and the characteristics of the 
fluid is given by the equation of pressure balance. In case of the assumed 
two dimensional geometry (0/02= 0) and for a quasi-stationary process this 
can be written in the form (*): 


: pB? 
(8) grad (7 + r) = 0. 


To connect the material pressure with other quantities characterizing the 
fluid an equation of state is needed. For example if the process is adiabatic: 


f 1/x 
(9) v=% (22) 
p 


where V is the plasma volume, p, and V, are initial values of pressure and 


volume, respectively. x is closely connected with f, the number of degrees of 
freedom in the plasma: 


(10) 4 = 


The expression for the variation of temperature with pressure is corres- 
pondingly : 


(11) T=T, aie | 
Po 

Equations (7), (8), (9) and (11) give a complete description of the processes 
taking place in the plasma when the external magnetic field is varied. 

Let us apply these equations to the case where the plasma is initially con- 
fined within a material container and no magnetic field is present at t= 0. 
Due to eq. (7) the magnetic flux within any closed material curve remains 
constant, in this case p=0. Consequently the magnetic field inside the 
plasma remains B,,=0, while outside the fluid it is B,,— B(t). Along a 
line crossing the boundary the magnetic field is described by a step function. 

The «step» at the fluid boundary means that surface currents have been 
induced in order to prevent the penetration of the external field. These sur- 
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face currents interacting with the external fields produce a force density 


f=ixB directed into the fluid tending to compress it (Fig. 1).. This is ex- 
pressed by eq. (8) which yields after in- 
tegration: 
+ x noe XX x Rte 
he $ Rae oe. DEI 
(12) È + p =0 a Fluid' 
<p È à ae x x À 
Choosing the point €, inside and &, out- x x xx 
. . . x 5 Ce 
side the fluid, and observing that B(6,) = oe ie 
= B,,—0 and p(&)=0, eq. (12) becomes: Sl Sl 
3 = i mete ” ho 
(13) BTL por Ree he 
aU 
Fig. 1. — Conducting compres- 
As soon as B?(t)/2u >p, the fluid de- sible fluid in a growing magnetic 


taches itself from the material walls and is 


field. 


confined merely by the magnetic pressure. 

Since B.,(t) is known, p,,(t), V(t) and 7(t) are also known functions through 
eqs. (13), (9) and (11). As B(t) increases the pressure and temperature of the 
fluid increase as well, while the volume decreases according to the laws of an 
adiabatic compression, where the magnetic field assumes the role of a piston. 


3. — Independent particle model. 


When the plasma is so rare that particle collisions can be neglected, an 
entirely different physical picture must be adopted. Each particle moves in 
the electromagnetic field produced by external sources and by moving plasma 
particles. We therefore investigate the motion of a plasma particle under 
these conditions. Contrary to paper I where the stationary state of such a 
plasma was treated in general, we shall here restrict ourselves to the «first 
order theory », considering nearly circular particle paths, with the corresponding 
cyclotron radius superposed by a small drift. 

The drift velocity is composed of different parts. In the stationary case 
every particle drifts along the lines |B|= const with the velocity: 


Lmv2 
(14) Wy = an (Bx grad B?), 


where the index M denotes that this drift is due to the (inhomogeneous) mag- 
netic field. The value of B is to be taken-according to the first order theory, 
at the guiding center. 
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When the external magnetic field is to be varied, the magnetic field inside 


the plasma undergoes a change and corresponding electric fields are induced. 
The equation of motion of a particle is: 


(15) = TE + vxB). 


The electric field is «transformed away » in a co-ordinate system moving 
with the velocity of «electric drift »: 


Ex B 
(16) ea 
Indeed, putting 
(17) v= w, +, 
eq. (15) becomes: 
(18) ded ux Bye Le, 


by expanding the vector product (Ex B)xB and making use of (EB) — 0. 
Expressing w, as a series expansion in the neighbourhood of the guiding 
center eq. (16) yields: 


ExB 
(19) Ww, = ( | + (R grad))w,; , 


where the index 0 refers to the guiding center, and R has the magnitude of 
the cyclotron radius. Due to the assumption of small field variations in 
the range of a particle orbit the second term is small compared to the first. 
The quasi-stationary treatment demands on the other hand that |w,|<]|v|. 
‘The first term on the right hand side in eq. (19) is therefore of the first order 
while the second term is of second order and can be neglected in our approx- 
imation. 

Eq. (18) has the following meaning: In a co-ordinate system moving with 
the velocity w,= (w,) the particle moves with a velocity u, affected only 
by the magnetic field and the inertial forces due to the acceleration of our 
system of reference. In the moving system — dw,/dt is equivalent to a grav- 
tational acceleration g and gives rise to the drift velocity (5): 


mgxB m 1|d /ExB 
20 urta LEE 
ara di q A B | B|. 


(5) L. Spirzer: Physics of Fully Ionised Gases (New York, 1956). 
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with E and B taken again at the guiding center. In the quasi-stationary treat- 
ment of slow time variations w, is a second order term. 

It can be shown that, with Suitably chosen «material surfaces » eq. (7) 
is again valid. Attaching the «material curves» to the guiding centers of 
particles, corresponding at a given moment to the same |B |= const lines 
(Fig. 2), the rate of change of flux is: 


(21) Gi [pas = 4 [Bas + f Boyds) + f Bow, ds) 3 
s s 


where second order terms are negelcted. In consequence of our choice of 
material curves w,xds=0. The first term on the right hand side can be 
transformed into a line integral with 
the help of a Maxwell equation; 
further substituting eq. (16) into 
the last term and expanding the 
double vector product, we obtain: 


(22) i? =—PEds +p Eds, 


where we made use again of (EB)—0. 
Eq. (22) shows that the magnetic 
field between any two |B|= const 
lines is frozen into this region. 
Apart from this restriction regarding 
regions for which eq. (22) is valid, 
the analogy with eq. (7) is apparent. 

As it was shown in paper I, the 
pressure balance equation is valid for 
any stationary plasma irrespective 
of collision rates. Eq. (8) can be used therefore without restrictions also in 
this treatment of non-colliding particles. 

To complete the set of equations governing the plasma in varying fields, 
an equation of state is again needed. This can be obtained, by further inves- 
tigation of the particle motion. 

Let us assume for example that the external magnetic field is increasing. 
To satisfy condition (22), the regions S must shrink and the plasma becomes 
compressed. As the flux crossing decreasing material surfaces tends to be. 
constant , B is increasing even in the co-ordinate system moving with the 


Fig. 2. — Tenuous plasma in“ a varying 
magnetic field. 


43 - Il Nuovo Cimento. 


3627 


666 G. SCHMIDT 


guiding center of the particle. In this system the particle has a spiral path, 
gaining energy from the field by a betatron like mechanism (Fig. 3). 
The gain in energy € of a particle in this system is: 


4 
FA Las 


2 
8) dira 


(23) de — d | dI 


a) where p,, is the magnetic moment defined by 


(24) Um = =p 


It is known that 4,, is a constant of motion. It 
is worth-while to point to a physical meaning of this 
fact. Eq. (24) can be written in the form: 


A 


Fig. 3. — Path of a plasma 
particle in avaryingmag- (95) 1mu? _ Lg? (R:1B) 
netic field: a) in the 2 VB om M 

guiding center co-ordi- 
nate system; bd) in the 


laboratory system. The bracket expression is just the flux enclosed by 


the particle orbit. This means that the magnetic field 
is frozen into each particle orbit, providing again 
an equation of the form of eq. (7). The «material curve » is in this case the 
particle path. 
It follows from eq. (23) making use of the constancy QE us: 


(26) de =p, 08, 


(27) 


that is the particle energy is proportional to the magnetic field measured 
always at the guiding center. From this follows furthermore that picking up 
a bunch of particles in a small region of the plasma, the energy of each of 
them will be multiplied by the same factor B/By owing to the same drift com- 
ponent in the direction of erad|B|. (Only w, possesses such a component in 
the first order theory). Consequently a local Maxwellian distribution is main- 
tained during the process, for the “particles concerned, and eq. (27) can be 
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written in terms of two dimensional kinetic temperatures, 
(28) — = 


The plasma pressure can be expressed in terms of the particle density and 
of the kinetic temperatures (see for example eqs. (34) and (25) of paper I}: 


i; My A es 0 
a ae he (7) a 


For each material region for which eq. (22) is valid we may rewrite it in 
the form: 


V No By 
(30) D ne 


By combination of eqs. (28), (29) and (30), we obtain directly the equa- 
tions of state: 


(31) Vay; (di 
and 
Cer fo 
(32) r=r,(£) 
Considering that the number of degrees of freedom is f == 2, eq. (10) yields 


x — 2 and it can be readily seen, that eqs. (31) and (32) agree with eqs. (9) 
and (11) respectively. With this, the formal analogy between the behaviour 
of plasmas in the two entirely different limiting cases in completed. 

Otherwise it can be shown in a more general—although more formal—way 
that a tenuous plasma in arbitrary electromagnetic fields obeys always adia- 
batic equations of state (*). Namely for such a plasma the Boltzmann equa- 
tion reduces to df/dt — 0. This means that the density of any bunch of par- 
ticles in the six dimensional (r, v) space remains constant in the course of 
motion. As the probability, of finding a particle in the six dimensional unit 
volume is uniform in the phase space the thermodynamie probability of any 
bunch of particles remains constant as well. This is equivalent to the con- 
servation of entropy of any part of the system, proving our statement that 
the thermodynamic transformations concerned are adiabatic. 


(5) This was derived recently on different lines by I. B. BERNSTEIN: Phys. Rev., 
109, 10 (1958). | 
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We turn now anew to the problem of the behaviour of an initially field 
free plasma in à growing magnetic field, treated in Sect. 2 in the magneto- 
hydrodynamic approximation. In the first stages of the process, when the 
applied field is yet small, the particle paths become slightly bent. Only in the 
later stages are the previous considerations applicable, in case that the fol- 
lowing restrictions are satisfied: 


1) The cyclotron radius of any plasma particle is much smaller than the 
linear dimensions of the vessel. In.this case, the great majority of particles 
are confined by the field and not by the material container. 


2) The usual condition for the first order treatment is satisfied. This 
requires that the variation of the magnetic field should be small inside a 
cyclotron orbit. 


Tf the magnetic field is sufficiently large (and È small) the treatment above 
becomes applicable, provided that the time rates are sufficiently small for the 
quasi-stationary treatment. We may inquire now about the variation of various 
plasma characteristics in terms of the external magnetic field B(00, t). 

The combination of eqs. (28), (29) and (30) leads to: 


B? 
(33) DATE Be ) 
which can be rewritten in the form: 
B?*(€, t) 
34 | Acs 1) Pee 
( ) p(E, ) P(E0; Bet, 0) b) 


where the new notations express clearly the fact that quantities belonging to 
different instants refer also to different space co-ordinates. The co-ordinate § 
equals & at t= 0, and follows the motion of the guiding centres of the par- 
ticles. 

To express the variables concerned in terms of the external magnetic field 
B(co, t), we use the integrated form of the pressure balance equation (12), 
which yields for any point § at the instant t: 


BE, t i 
(35) PE) sea 
For t= 0 eq. (35) reduces to: 
pB? ae 0 2 
(36) ae + p(B, 9) = oe . 
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By combining eqs. (34), (35) and (36) we readily obtain: 


(37) B(E, t) _ Bloo, t) 
BE, 0 Blog, 0)” 


This is a very important result. It means that, increasing the external field 
to its » fold value during the time t, each particle experiences a » fold increase 
of the magnetic field in its own moving co-ordinate System. If the initial 
distributions at &— 0 are known, each quantity may be readily obtained as 
functions of B(oo,t) only. From eqs. (34) and (37) we have: 


B*( co, t) 
(38) P(E, t) = p(Eo, 0) ) Bi(00,0) | 
Combining this with eq. (32) we obtain: 
fo (00, t) 
39 #4 = 2 
( ) (E, t) TED B( , 0) 


or with eq. (31) the volume occupied by any bunch of particles is: 


B(co, 0) 


(40) Hie, t) a V (Eo, 0 B(co, t) . 


APPENDIX A 


Supposing we produced somehow a partially or totally ionized « cool » 
plasma in a cylindrical container, then both temperature and pressure are 
enhanced by application of a growing magnetic filed. 

Assuming an initial temperature of 7,~10° °K and pressure of p,10-3 atm 
and a magnetic field of about B,,; —16000 gauss peak value (actually a 
rather conservative estimate), the corresponding material pressure after the 
compression is roughly: 


Bou N 
(41) Pmax = & 105 mnt w10 atm. 


24 


Considering a «two dimensional » compression (x = 2), the corresponding tem- 
perature is from eq. (32): 


(42) dI cd (==) DS 10° °K : 
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The plasma volume is during the compression compressed to 
\4 
(43) LR ( — | & 


A «three dimensional » compression yields with x=3 
Ta es A108 °K and Vaan  4°10-* Vz: 


Which x is to be chosen, depends upon the compression time T. 
The «self collision time», the average time required to reduce any lack 
of anisotropy in the velocity distribution is (5): 


ae 11.4 Ai T# 
(44) 12 in 4° 
where A = ™/Mbprotony % = particle/cm?, and In A is a slowly varying function 
of T and n, and it may be taken In A & 10. In a plasma of deuterium ions 
Z=1 and A=2. Assuming a linear increase of B with time, and x= 2: 


t t 
B = ee Deo ; fis EME 7 re s M 
T À T T 


The probability of a considerable change of the distribution function by col- 
lisions during the compression, is: 


T T 


di ae In An o, (5 fo In A2 Nmax 
be RE 11.4:v2T%, a Fe 11.47%,. 


Do, 
0 0 


(45) 


The assumption of a «two dimensional compression » is justified until this 
quantity is much less than 1, or 


11.47%, 


46 ca 
x = In A 1/2 Ninax 


& 4-10 68. 


For electrons however, this quantity is much less because A #&1/2000. In 
this case we obtain the condition: 


TOUR ASE 


Therefore it is obvious that in case of t< 10-78 the compression is two di- 
mensional, when t>4-:10-*s, it is three dimensional. In the intermediate 
range, the situation is rather complicated. It is generally impossible for the 
two types of particles to undergo a different type of compression, because 
this would result in net space charges and immense electrical field strengths. 
The two types of particles are furthermore coupled by the common magnetic 


field determined by currents of both types of particles. For this « mixed » com- 
pression §<x%< 2. 
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For t>4-10- there is a Strong collision interaction, and the magneto- 
hydrodynamical model describes the phenomena. There exists an upper limit 
for t, because of the diffusion of magnetic field into the plasma, due to the 
finite conductivity. We can make an estimation for this with the help of 
eq. (6). If we approximate V2B ~ B/L* (where L is of the order of mag- 
nitude of linear dimensions) and substitute L2 for S, we obtain: Z 


(47) Si B Tr? & ag Lo 
po L? di 
Since BL? = D we get 
(48) ca RI S Ê 
D © uoL? 


The time constant for diffusion is therefore 9 ~ uoL?. The conductivity 
of a plasma is o = T?/653 (5). In case of a cylindrical geometry L? is pro- 
portional to V. With 


3 3 
= =" Lar | Li | and L? = far (=) . 
p 


We obtain for #: 


Tt da 
49 9 = yw — max min 


This shows that # is the same constant in every phase of the compression. 
An initial value of L, — 0.1m results ¢.g. Lin & 4-10-5 m? and 9 & 6:10-4 5. 

In the collisionless case the size is limited by the condition L>R. With 
x—=2, and the same cylindrical geometry: L? = Diin(Bmax/B). It can be easily 
shown that L/R is a constant during the compression, since combining 
R=mv/qb with eq. (24) we obtain by means of an elementary calculation: 


/ 
JE 


2UmM 


16 / 
50 —— L in 
(50) R = a] 


For deuterons with B,,,. = 16000 gauss, Rui, & 4 mm. A final plasma diameter 
of some centimeters yields an initial size of some 10 cm, which is quite rea- 
sonable. 

The assumption of an adiabatic compression is justified only when energy, 
losses are negligible. Indeed, assuming a compression time of about 1 us, 
the loss by recombination radiation per particle is of the order of 10? eV, 
the bremsstrahlung radiation is 10-*eV whereas the energy gain during the 
compression is (10?—10%)eV. The neglect of particle loss by recombination 
is also well justified, the rate during the contraction amounting to some 
10~ recomb/particle. 


ni 
ri 


a UM Sn, tee NG Ke Fa Les 


LR. sd LE Um We. 
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APPENDIX B 


Finally we investigate the case of high compression velocities, when the 
velocity of the « magnetic wall» reaches the order of magnitude of the par- 
ticle velocities in the plasma. We excluded this case in our former treatment 
in order to avoid «non-quasi-stationary » processes. 

We investigate in the following the magneto-hydrodynamic model in the 
ideal case of no magnetic field penetration. The phenomena in this case do 
not differ from the simple process of compressing à gas in à cylinder with help 
of a moving piston, only the piston is substituted by the « magnetic wall ». 

Assuming now a piston in periodic operation of compressing and expanding 
a gas, it is easy to see that this process can be regarded as quasi-stationary 
only as long as the piston velocity V is much lower than the molecular velo- 
cities of the gas v. When however V reaches the order of magnitude of è, the 
energy gain of the gas in the compression period—when the piston is running 
towards the gas—is larger, while the energy loss by expansion—when the 
piston is running away—is smaller than in the quasi-stationary Case. The net 
effect arising from this irreversible phenomenon is, that in a cyclic operation 
the gas is heated to ever increasing temperatures (7). 

A statistical explanation of this fact is the following: Given an ensemble 
of microscopic particles interacting with each other and with the piston. The 
average kinetic energy of the piston, {4MV?) corresponds to a very high 
temperature, due to the high mass M. The particles tend to reach the tem- 
perature equilibrium with the practically infinitely « hot » piston, by obtaining 
kinetic energy from it (*). 

We investigate on a simple example how this process works. Assuming 
a one dimensional problem, with not interacting particles, and two synchro- 
‘nously moving walls with a momentary distance L(t), the wall velocity is 


4 1 dL 
51 $2 ees 
(51) | V ee 


If the particle velocity v is small compared to V, we obtain for the time of 
flight between two reflections: 


(52) LS L-EVa 
v 
or 
(53) ess edie: 3 
Dh 


(7) The «magnetic pumping » mentioned in ref. (4) seems to be similar to this 
method. 


(8) G. Karman, L. Pocs, G. SCHMIDT and K. SimonyI: Periodica Politechnica 
Budapest (Electr. Eng. Series), 1, 53 (1957). 
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The number of reflections in the time dt equals dt/t while the gain in ve- 
locity is 


(54) SNS LP 


In case of v/V + co denoting v in this asymptotic case by «, we have 


du V dL/dt 
5 = 9 no di, 
(55) = zu 7 di 


After integration we obtain 


L 
(56) u = UT : 


where « and Z, correspond to the same # =. As u?=u5(T/T,) and Lp/T = L,p, (ae 
from eq. (56) follows 


p $ 
57 T=T (2) È 
(57) "ln, 
As for a one dimensional case x =3, eq. (57) equals eq. (11). For the quasi- 
stationary case we put v= u-y in eq. (54) and obtain 


” 2 pe 
(58) y du + u dy =—2V + dt + dt; 


with help of eqs. (55) and (56) it follows: 


9 
5 — —_ |/V?dt. 
(59) {LO nil 


From this we obtain finally 
t 
Uo Ly 2 ri 
0 


We see that in the irreversible case v is always greater than in eq. (56). In 


our case of « Magnetic piston » we control the pressure instead of L. Expres- 


sing again L and v by p and 7 we obtain from eq. (60) 


py p vat So 
(61) a (5) +240 (5) + ¢ == 
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applying the following notations: 


TÉ 2 


0 


AS — 2 (va. 
Po Ly 


As A and C are both positive, it is obvious, that the same T corresponds now 
to a smaller p/T and therefore to a smaller p then in the reversible case 
when C=0. 


RIASSUNTO (*) 


Si esamina il comportamento di un plasma posto in un campo magnetico esterno 
considerando due easi limite: 1) il caso idromagnetico ideale (0 > co); 2) il caso privo 
di collisioni. Restringendo i calcoli a geometrie cilindriche generali (0/0: = 0) e a 
processi quasi stazionari si ottiene per i due casi un sistema completo di equazioni 
del plasma. Vi sono incluse tra altre equazioni di stato anche equazioni per plasma 
tenuo. Esiste una notevole somiglianza tra i due sistemi di equazioni riferentisi a con- 
dizioni fisiche totalmente differenti. Le equazioni si lasciano prontamente risolvere 
per esprimere ogni quantità macroscopica di plasma in termini dell’intensità del campo 
magnetico esterno. Si descrive in breve l'estensione delle nostre considerazioni a processi 
quasi non stazionari. Ciò conduce a un fenomeno che suggerisce nuovi metodi per 
aumentare le temperature del plasma. 


(*) Traduzione a cura della Redazione. 
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On the Inversion Problem for a Klein-Gordon Wave Equation. 


V. DE ALFARO 


Istituto Nazionale di Fisica Nucleare - Sezione di Torino 


(ricevuto l’11 Agosto 1958) 


Summary. — We treat in this paper the problem of obtaining the static 
potential which corresponds to given resonance parameters in the case 
of a Klein-Gordon wave equation. Our results are an extension of those 
already found in the non relativistie case. We restrict the class of our 
potentials to those vanishing outside a finite range a. The knowledge 
of the resonance parameters allows to obtain an approximate solution 
of the so called Gel'fand-Levitan equation which provides the desired 
potential. 


1. — We consider a Klein-Gordon equation for S-waves: 
(1) p'(E, x) + [K? — 2E V(x) + V'(x)]q(E, x) = 0, 


with the purpose of deriving the potential V(x) when its spectrum is known; 
the most direct way is provided by the so-called Gel'fand and Levitan equation 
corresponding to eq. (1). 

This reads (*): 

x 
(2) K g(x, Y) = Tas(æ, y) a yO; t) Let, y) dé ) 
0 

(*) In this equation x, 6 have two values 1, 2. Its first derivation is due to CoRI- 
NALDESI (1), who has limited himself to «=1 and B=1, 2. We have followed here the 
formulation of VERDE (?), by which Ky, is a 2x2 matrix. The reason for the appearance 
of it is due, as we will precise below, to the fact that the Klein-Gordon operator is 
not self-adjoint and the completeness relations require two-component eigenwaves. 

(1) E. CORINALDESI: Nuovo Cimento, 11, 468 (1954). 

(3) M. VERDE: The inversion problem in wave mechanics and dispersion relations, 
to be of published on Nuclear Physics. 
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PAU) contains the whole information, being defined by 


aL 


(3) fapl@, Y) = + [ran Ê = aH DE, 01 PE, Ye - 


Fig. 1. — Path of integration. The com- 

plex plane represented here is the &-plane. 

E and k are single valued functions of £: 

E = m cosh (£/m); & =m sinh (£/m). The 

whole Æ plane is mapped in the region 

— 3n/2 <Im £<x/2. @ Zeros of f(K); 
x Zeros of f(— K). 


G(E) 


Here D,(E, x) is the wave matrix (de- 
fined later) whose first component is the 
eigensolution of (1) in absence of the 
potential normalized so that for æ — 0, 
D (E, a) ~ &. G(E) is related to the 
so-called spectral function as we shall 
specify later on. | means that the path 


of integration has to be taken as in 
Fig. L. 

From eq. (2) the potential V (2) fol- 
lows, being 

K,(0,%) le 

MR ee IA È, i y(0)]. 
Clear e oe Oe 
We wish to approximate the f_;, once a 
set of resonance parameters is experi- 
mentally obtained, so that it yields a 
simple solution for eq. (3). 


2. — In agreement with Jost and KOHN (?) we define the two independent 
solutions of (1) f(k, v) and f(—k, 2) with the asymptotic behaviour 


f(k, ©), ~ EXP [— ke], 


f(—k, x) XP [ike]. 


The solution of (1) which for æ +0 goes like æ can be expressed as: 


1 
Be) = SEE) CR; wv) — f(—k) f(h, @),] 


f(k) = f(k, 0). The S matrix and the phase shift are defined by 


S(k) = exp [216] = 


{(k) 


{(—k) 


(3) R. Josr and W. Koun: Dan. Mat. Fys. Medd., 27, no. 9 (1953). 
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Two solutions of (1) corresponding to energies 7 and £’ Satisfy the quasi- 
orthogonality relation 


(ce) 


(5) | LE + B27) 9H, 0) gE’, à) dx = sE — PI) GE, 
ò 
where G(E)= f(k)f(—k) and |E|>m. 

We know that g(E, x) is an entire function of Æ of type æ (*). f(k) however 
is a two valued function of k, having two branching points at k= + im. 
The same happens to E considered as function of %. This Suggests that G(H) 
might be an one valued function of Æ because it is even in & and it has the 
Same branching points. A close examination of G(E) shows that this is in- 
deed true. 

The set of functions g(E, x) is two times complete, i.e. it can be used to 
expand two functions. This is clear when one considers the classical limit, 
in which one retains only the positive energy solutions of (1) which form a 
complete set in the interval 0 -:: oo. It is customary to define two com- 
ponent wave functions ®(E, x) defined by 


| g(E; 2) | 
(6) DE, a i 


E — V)9(#, x) 


On this set the scalar product is defined as 
(7) (OE, 2), DB, 2) = | D (E, o, DEI, a] de, 
0 
where 0, = : 0): The set (6) is orthogonal in the sense of (7). It is also 


complete because of the expansion 


+ © 


(8) dpt — a) = | do(E) D,(E, «)[o, DE, x’): 


— oo 


(*) Let F(k) be an entire function of k and M(h) the maximum modulus of F(k) 
in the circle of center k=O and radius h. The order of F(k) is defined by 


li log log M(h) be 
Due log h ee 


The limit lim sup h-@ lg M(h) is the type of F(k). See R. P. Boas: Entire functions 
ho 
(New York, 1954). 
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Here the spectral function is defined as 


do(#) __F if |E|>m; 
dE a) 


do(E) ui 12 ei, - a rt , 
— = 5, |2| (Hj =" Pdxl ÔEŒ—E) for |E|<m 
oF) _ x, bfe-nialde) à : 


E, and y, denote the energy and eigenfunction of the 1-th bound state. 


A more compact form of (8) can be obtained by going into the plane of 
the complex energy 


kd Mt 
date = [gy Pa Bs olo: DES) 


e 


3. — Many of the results of the non-relativistic analysis can be generalized 
without any difficult to the Klein-Gordon equation. It is obvious that G(E) 
is an entire function of E of type 2a, where a is the range of the potential. 
The same approach as in (*) yields quite similar results concerning the distri- 
bution of the zeros of G(E) and the corresponding expansion of this function 
into an infinite product. In dealing with these properties we observe that it. 
is by far convenient to use as independent variable the energy E instead of 
the wave number k, just because q(E, x) is not a uniform function of k, while 
it is uniform in FH. 

We are now interested in the classification of the zeros of G(E). The real 
zeros of G(E) belonging to the interval — m --- m imply the existence of par- 
ticle or antiparticle bound states. The remaining zeros either belong to some 
zero norm pathological states (*) or they are the energies of the so-called re- 
sonance states (in infinite number) or finally a few of them define the virtual 
states. Apart from the bound states and the zero norm states the solutions 
which correspond to the other roots are not acceptable from a strict quantum 
mechanical point of view because of their non summability; yet in many 
practical cases (e.g. in nuclear physics) the complex energy resonance states 
are treated as true states with some individuality. We want to investigate. 
the completeness of the set S of ®(E, «) which correspond to the zeros of G(E). 


(4) T. REGGE: Nuovo Cimento, 8, 671 (1958); Nuovo Cimento, 9, 491 (1958). 

(*) If the width and depth of the potential are large enough some states of vanishing 
norm may occur (5). Although their spatial behaviour is quite correct their energy is. 
complex. We shall not consider this case as it leads to unwanted difficulties. 

(5) L. I. ScHIrr, H. SNYDER and J. WEINBERG: Phys. Rev., 57, 315 (1940). 
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We shall find that this set S is complete in the interval 0 --- 2a. This is the 

Starting point for an expansion of the functions fig Of (3) in terms of the 

set S, in order to obtain a form for them which is practical for computation. 
Let ®(E,,x) be any function of the set S. If 


H 
5 n 


> | D, (En, x)[0, E(x)] CHANT) 


a 
0 


for every n implies £(x) — 0 in 0---H, then $ is closed in 0 --- H. 
We introduce the entire function I(E) of type H defined by 


H 


IE) = | D, (E, æ)[o, E(e)], dr, 


I(E) can be factorized as follows: Z(E) = G(E)R(E), where R(E) is an entire 
function. Since the type of G(E) is 2a the type of I(E) is at least 2a. There- 
fore unless H > 2a we must take I(E)= 0, R(E)=0. The conclusion is that 
the set S is complete in the interval 0 --- 2a. The independence of the set 
follows from the existence of another auxiliary set of functions 1,(a#) de- 
fined by 


(9) SOR =| 0.08, x)[0, Tn(@) la de . 


0 


Eq. (9) is a generalized Paley-Wiener theorem and its detailed proof requires 
some care and some additional considerations. From (9) we deduce 


It is evident from (9) that r,(x) is orthogonal in the sense of (7) to all ®,(E, x) 
but ®,(E, x), and it vanishes for x > 2a. Any relation of the kind 


[col 


Y_ ADE, à) =0 


0 


holding in the interval 0,; - - 2a implies therefore that all 2, =0. 
The completeness and independence of the set S in 0 --- 2a imply the sa- 
me for the set ®,(H,, 7) of the unperturbed solutions. The set 


G(E) 
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satisfies orthogonality properties : 
2a 
| 9,.(æ)Co, Po(Em de = Om (Em) - 
ò 
4, — The completeness of the set D,(E,, 1) allows to obtain an approximate 
formula for the kernel f,, of the Gel'fand and Levitan integral equation. We 
can expand f,, using the set of ®(L£,, 2): 


; PACE Y) me + DAT) [01 0 (LE, Y)] ? 
; where Wnj(@) is given by: 
% 2a 


ea 
na (x) = ct) | Poi dy i 


ul 1 
= a] Ing(4) ay | var [ - DAE, x) 0, DE, y)le = 


0 


1 rl 
== al k dé È œil Dos (É, x) NED: A 


Taking into account that 0,(x) satisfies the equation 


G(E) — 


| a 4 
Ÿ'a(x) + fo(— Kn x) “si =| ban oi E. DE: x) ’ 


€ 


‘ley 
exp [ine] | we obtain for 


where fo(— kn, 2) is defined by: fo(£n nX) = 
(E,— V) exp [ik,®] 


fat, y) the expansion 


Al 
fap (x, y) = > G'(E,) [d,(2) +- fol Kn, x)],[0, P(E, y) |p : 
But 
1 
2 aH, On l#)[O, PolHns Y) Ia = Ô,pte—y}= 0 if a>Y 
sì so that the expansion assumes the simple form 


ho 1 
x fap(&> Y) er > G'(E ) fox kr; æ)[o DE»: Y)le 5, 
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If æ < y we can exchange # and y in the formula in virtue of the symmetry. 

With this expansion we can construct PACA y) from the experimental data. 
This form of f,,, as it is evident, is especially suited for the integration of 
eq. (2), since our approximation of f,, consists in a finite sum of products of 
functions of the variable x or y only. 


* *% * 


I wish to thank Prof. M. VERDE for many helpful discussions and Prof. 
T. REGGE who suggested this problem. 


RIASSUNTO) 


In questo lavoro trattiamo il problema di ottenere il potenziale statico corrispon- 
dente a dati parametri di risonanza nel caso di un’equazione d’onda di Klein-Gordon. 
I nostri risultati sono un ampliamento di quelli già trovati per il caso non relativistico. 
Restringiamo la classe dei nostri potenziali a quelli che si annullano oltre un range 
finito a. La conoscenza dei parametri di risonanza ci consente di ottenere una soluzione 
approssimata della cosidetta equazione di Gel'fand-Levitan che fornisce il potenziale 
richiesto. 


(*) Traduzione a cura dela Redazione. 


44 - TI Nuovo Cimento. 
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The Polarization of Knock-on Electrons Arising 
from Scattering by Polarized |.-Mesons. 


D. FiveL and T. FULTON 
The Johns Hopkins University - Baltimore, Md. (*) 


(ricevuto l’11 Agosto 1958) 


Summary. — An experiment is proposed for the direct measurement of 
the helicity of u-mesons arising from x decay by studying the polarization 
of the high energy electrons produced in the collision of the polarized 
u-mesons with unpolarized electrons. Polarization effects become signif- 
icant only for high energy incident u's (of the order of a few GeV). 
The separation of a polarized beam of u-mesons from the decays in flight 
of high energy x’s is in principle made easier by the fact that the labo- 
ratory momenta of u-mesons which are emitted in the forward and 
backward directions in the 7 rest system differ by about 42% as a con- 
sequence of the relativistic transformation from the center of mass to 
the laboratory system. Expressions in the laboratory frame are also 
given of the cross section for scattering polarized electrons by polarized 
u-mesons. In the C.M. frame, the results agree with those calculated 
by Bincer. 


1. — Introduction. 


It is well known that knowledge of the helicity of u mesons produced in 
the decay of the 7 would provide considerable information concerning the 
coupling involved in the decay interactions. Techniques used to measure the 
helicity of electrons in f-decay, which utilize spin effects in Rutherford scat- 
tering and bremsstrahlung, can probably not be successfully applied to  me- 
sons because of their relatively large mass (7). Suggestions haver ecently been 


(*) Supported in part by the U.S. Air Force through the OSR. 
(1) J. FRANKLIN and B. MARGOLIS (Phys. Rev., 109, 525 (1958)) discuss the pos- 
sibility of utilizing Mott scattering of u--mesons by nuclei to determine the y helicity. 


3641 


THE POLARIZATION OF KNOCK-ON ELECTRONS ARISING FROM SCATTERING ETC. 683 


made to study the w helicity by observing the asymmetry of the electron 
angular distribution in the $-decay following the capture of polarized uy mesons 
in the 2C present in nuclear emulsion (?). Depolarization effects of various 
kinds, and the low energy of the recoiling !*B nucleus (used in identifying the 
event), make this experiment a difficult one to perform. Another possibility 
which has been discussed (?) involves the knock-on scattering of polarized 
u mesons on aligned electrons. It is suggested that u mesons, appropriately 
selected to be only from forward decays in the rest system of the x, be made 
incident on a block of iron magnetized so that a fraction of the electron spins 
are successively parallel and antiparallel to the » beam. The cross-section for 
knock-on electrons is then to be compared for the two cases. 

= We are also concerned with the knock-on process in the present note. If 
one scatters unpolarized electrons by polarized u mesons, the final electrons 
are polarized. We obtain expressions for this polarization in Sect. 2. Since 
spin dependent effects arise from a spin spin interaction which is of much 
shorter range than the Coulomb interaction between electron and w meson, 
the spin dependent terms can be expected to become important only for 
u mesons which have laboratory energies large compared to the rest mass of 
the meson. The high energy electrons arising from a knock-on scattering by 
u mesons of a few GeV are sufficiently polarized that the polarization of the 
electrons in the resulting shower (4) may possibly be detected. The low cross- 
section for u-e scattering at the high incident meson energies required may 
make detection of the polarization difficult. 

A partially polarized beam of high energy y mesons can be obtained by 
selecting a fraction of those which have appeared from the decay of a beam 
of energetic 7 mesons of uniform momentum. For example, 7 mesons with 
a total energy of 3 GeV result in 3.003 GeV/c momentum for the 4 mesons 
which were emitted in the forward direction in the x center of mass system, 
and 1.726 GeV/c for those emitted backward. 

In the Appendix, we list the results in the laboratory frame for p-e knock-on 
cross-sections when both particles are initially polarized. These results, when 
transformed to the center of mass system, agree with those of BINCER (*). 


2. — Calculations and results. 


We consider the final state density matrix for the process in which a 
u meson, of four-momentum Yi = (Pi; 41), kinetic energy £ and mass 
M (ce, = E+M), scatters an electron of mass m at rest (pe: — (0, im)). The 


(2) J. D. Jackson, S. B. TREIMAN and H. W. Wy Lp Jr.: Phys. Rev., 107, 327 (1957) 
(3) A. M. BINCER: Phys. Rev., 107, 1434 (1957). 
(4) K. W. McVoy and F. J. Dyson: Phys. Rev., 106, 1360 (1957). 
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a / x è È A 
final state is a p meson of four-momentum p, and an electron of four-mo- 
= ! ! ! ER A 
mentum p,, and polarization L: An energy E' (E'=e,—m) is imparted to 
the electron. We define the four-vectors (°) 


E A 
(1) a= Lay ela 
and 
ES 1 (Li Ps) ! sharp 
(2) È (1 RP m(e, +m)??? Um | i 


P, is the percentage polarization of the y and is positive or negative, depending 
on whether it is parallel or antiparallel to its direction of incidence; pisa 
unit vector in the direction of incidence, and is v/e for the incoming vu. 
A simple calculation, which utilizes the positive energy projection operator 


(3) A+(s, p) = (4m)1(1 + èysy :8)(m — ty-p) 


for taking into account the specific spin orientation, leads to a final state 
density matrix 
! ! é 1 
Ole, + & — Ei — E2) dps Y 
[ (pi pa)? — m?M?} ep PD) 


(4) es = rim! 


In the expression above, 7, is the classical electron radius, and Y is given by: 


(5) Y=M+Y, 

where 

(6) Y, = (m2 + M*)(m? + psp.) + (pipe)? + (pi Po)? ; 
(7) Vi — mM[s1°8,(p, — DE + 81° (Pr — pi) 82° (Pe — P,)) - 


The differential cross-section and polarization at final electron energy EF’ 
are given by (5) 


(8) Co == ARR OF ’ 
and 
ei 
(9) Pe AL 
Leo 


(5) H. A. ToLHorK: Rev. Mod. Phys., 28, 277 (1956). 
(5) L. WoLFENSTEIN and J. ASHKIN: Phys. Rev., 85, 947 (1952). 
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where it is understood that we let A => and that the traces involve sum- 
mations or integrations over all final variables with the exception of the 
electron energy. 

If we specialize to laboratory co-ordinates, we find that the differential 
cross-section has the well-known form (°) 


27nmrs |, ° Hi bs) ae CARA, 
10 E, BH’) dH’ = |1 ae LE’ 
(10) CE, B')¢ 3° | Pe tara) | ‘ 


where E, is the maximum energy which can be imparted to the electron. This 
value of E’ is reached for electrons scattered in the forward direction. Equa- 
tion (10) is the cross-section calculated for unpolarized particles. 

The component ot polarization along the incident direction, expressed in 
laboratory co-ordinates, is given by 


P,( E' La 2m—H 1, ut 4m 
\NE+ M ro om CT En 2(E + M)} E'+2m 


MEI au: de 
E, 2\E+M 


(11) P(E,E)= 


. 


The average polarization is of course entirely along the incident direction. 
Tf only knock-on electrons between energies of Æ,, and nE', are detected, this 
polarization is defined by: 


Em 
Pi|oE, E')P(E, E')dE' 
6 | D,,|P1 io. nEm 
de) SRE nee Em 
[ou E')dE' 
nEm 


The behavior of J"; and 2, as functions of E is presented in the Table. 


We would like to thank L. MADANSKY for suggesting this experiment as 
well as the means of obtaining polarized beams of high energy p mesons. 


(7) See, for example, B. Rossi: High Energy Particles (New York, 1952) D. 16. 
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APPENDIX 


We. summarize here the results obtained for the differential cross-sections 
when polarized u-mesons scatter electrons polarized along or opposite to the 
direction of incident u-meson momentum. In the expression for the density 
matrix, Eq. (4), Y is replaced by Y’, where 


(A.1) = + I, 
and 
(A.2) Yi = — mM[3s,-8,(p, — pi)? + 81° (pi — pi) 82° (Pa — Pa) - 


In Eq. (A.2) above, s, is the four vector 


(A.3) Sat P.(P, 0) , 


and P, is the percentage polarization of the initial electron, positive if parallel, 
negative if antiparallel to the direction of the incident u. The differential 
cross-section in the laboratory co-ordinates has the form 


(A.4) o(E, E') dE’ =[o,(E, E') + o.(E, E')]dE", 
where 

2a mrs E' E E' 
A.5 (CITI > Pods al 
pe) oe BE" aul Etra) 


In order to get a measure of the spin CREER part of the cross-section, we 
define A,, where 


NEm 


The behavior of À, as a function of Æ is presented in the table. 
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Tagce I. — Table of results. 


| | ye ae A 5 tai A 
FA. ne | 2 .5 5 5 9 .9 9 
D PTE 8 (mb) | (%) | (% |(10-® mb) (% | (%) 
| 10 | OEE 4992 13 6.4 6.6 1.87 13.7 14.4 
20 | 729 .9977 2% 12.4 12.5 .46 31.5 34.2 
30 | 1480 .9996 | LE 18.1 18.3 .25 46.9 48.6 


RIASSUNTO (*) 


Si propone un esperimento per la misura diretta dell’elicità dei mesoni y derivanti 
dal decadimento x per mezzo dello studio della polarizzazione degli elettroni di alta 
energia prodotti nella collisione dei mesoni y polarizzati con elettroni non polarizzati. 
Gli effetti della polarizzazione appaiono significativi solo per y incidenti di alta energia 
(dell’ordine di qualche GeV). La separazione di un fascio di mesoni p polarizzati dai 
decadimenti in volo di x di alta energia è facilitato in linea di principio dal fatto che 
gli impulsi nel sistema del laboratorio dei mesoni p emessi in avanti e all'indietro nel 
sistema di riposo dei x differiscono di circa il 42% come conseguenza della trasforma- 
zione relativistica dal sistema del centro di massa a quello del laboratorio. Si danno 
anche espressioni valevoli nel sistema del laboratorio per la sezione d’urto degli elet- 
troni polarizzati dallo scattering su mesoni p polarizzati. Nel sistema del centro di 
massa i risultati si accordano con quelli calcolati da Bincer. 


(*) Traduzione a cura della Redazione. 
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(ricevuto il 18 Agosto 1958) 


Summary. — The radiative pion-nucleon scattering is studied in statie 
theory treating one meson exactly, and the other one in the one-meson 
approximation. Explicit expressions for the rescattering corrections are 
obtained. 


1. — Introduction. 


In this paper we shall apply the static meson theory of Chew in order to 
evaluate the matrix elements for radiative pion nucleon scattering, i.e. for 
the process: | 

T+HN>T+N4+4yY. 


This problem has already been treated by CUTKOSKY (!) who has been able 
to treat the incoming and outgoing pions in a symmetrical way. His method, 
however, leads to a very difficult evaluation of the rescattering corrections. 

In our paper we are able to obtain explicit expressions for these corrections. 
The method treats one pion exactly and the other one in a one-meson ap- 
proximation. It is especially adequate for situations in which only one of the 
two pions is in the (33) resonant state. The case in which both pions are on 


()-R. E. Curkosky: Phys. Rev., 109, 209 (1958). 
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the (33) resonance will be treated separately in a forthcoming publication. 
In this case indeed the photon energy k is quite small and therefore, following 
Low (2), one can derive a rigorous low energy theorem. 


2. — Radiative pion-nucleon seattering matrix element. 


The matrix element of our process is 


(1) A;(q; p) = (fesso doti) , 


where p and q represent the initial and the final meson respectively, & is the 
emitted photon and the other symbols are the same as in CHEW-Low's paper (°). 
As in the C.L. paper, we put: 
(2) FIA a Pigi 
with the condition: 
Cas 7,474) — [as , fr + js] Fi 0, 
(Po; Liv +is1?%) = (PTE) = 0. 


(3) 


Then by the same argument of C.L. we get for the vector current matrix 
element: 


RES 2 ari sd F,(k?) …, fimo/k ea) _ 
& ape eo ee |) = 
TA Uy Un (5 | F(k?) (wo vow?) 
2 ela ole ee ees 
where V° is given by 
Vt Ts De: v(w) 
1/20, 


and corresponds to a meson of momentum u=k /\e and isotopic variable DA 
We have thus to study the matrix element: 


po, VP), 


(2) F. E. Low: Phys. Rev., 110, 974 (1958). 
(3) G. F. CHEW and F. E. Low: Phys. Rev., 101, 1579 (1955). In the following 
extensive use will be made of the results and of the notations of this work, that will 


be quoted as C.L. 
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It can be expressed in the following way: 


70+ 1 7 2A 
(5) (Ea VETTA = CPs Ag vee) + CA) ine rs y 122 veg 


LH —o,—t 


0 (Wo, Vatu dba: Vets yD D (Pei Va PONE, Vase 
= Ome (Por Vato) > LL hp Op ve E,— «,— 1e ; 
where the identity: 


(6) a Pi = 0, Pi (EO, — = ie) Fe > gt 


has been used. 
We now introduce the new function 


(PO, 2 PSP) 8, (Boy POD) = (0 By. 


From (5) we obtain for (q-, p*), the following equation in the one-meson 
approximation: 


(Py, PAP, VER) (D, PO), VIP) 


Wa Wy — 1€ a 


pe x 


n (Way Ve Lo) Loy Ve) (Po, TPE SUR, Va) (Pos VR) PH) _ 


Op Wn — Wa — VE 7 O+— D, — de 
e x Fe Vecio P* )à > Ue, Ve RR, Va ote a 


Or — Ma — 1€ ro Or + Wy 


ee 


Pos re Ero (Ee, Vi) 


Or + Op — Wa — VE 


We point out that in order to obtain in (7) all the terms with the zero and 
one meson in the intermediate states, we need to consider in (5) also the two- 
mesons terms. In fact the two-mesons matrix element (Y,V,W,.) contains a 
ò-funetion term, as one easily recognize by writing for this matrix element 
an equation similar to the (5). When this 6-function is substituted in (5) it 
leads to a one-meson term (1). 

We should now study the matrix element FK;(q, p) due to the scalar current 
density. However, since we have the same reason as in photoproduction of 
neglecting this term in the final result, it will not have any consideration in 
our paper. 


(4) L. RoDBERG: Phys. Rev., 106, 1090 (1957). 
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We turn out to the matrix element 


(8) 2040) = (WO, [Fe AD aes”) = (un. far Cia) Fe 


Hi (M, P) 


+ (, [ar (— jz - AF) at”) a= > CA Vo* po) Ù 


n E, Wir ie 


Using the relation (6) to express a, P° the second term can be written 


i (—) 70+ gst 
(ir far (de Ana") =—3(¥., dv i ane) Var hs) _ 


E, + wr" Oo — 2 i 
since from (3) follows: 
(9) RAR a = 0 . 


Finally Eq. (8) takes the form: 


(10) 59,7) = C2 D | A (<i ag] 5”) = 


(POVHED) 


n + Da — Wp — ie 


x, VEE wn) Ae 


_ Pido(— j- AP 
| | (dada) )z On — D — €" 
where there are no contributions to the sum for n — 0 because of (9). 

We shall define 


(11) He (4; P) = (Po, [an (in At] TI 


The commutator appearing in this expression is the same as the one en- 
tering in the C.L. work on photoproduction problem; therefore we can write 


(12) 7° (4, p) = Hou (4 D) + Hin (4 P) è 
where 


Cr 5 + 
Feld; P) = an (Yo, [T2001 — 7% dw]S -eP, iF 


2ie,q ° 
AGP) = Fi (On de Des dde slo Y)} + 
DaDg+k 
20q+k 7 7 \ (+) 
deal, "a 7. (WoL , Bath dav Seek Ôw] Pr ) 
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3°. is coming from the interaction current and Æ,, from the conventional meson 
current defined by Eq. (3) of C.L.’s work. 

If we limit ourselves to the one-meson approximation, Eq. (10) becomes 
an integral equation for the matrix element Æ#7(q, p), since #6; (9, p) is known 
and the remaining term is given by the product of the solution of Eq. (5) times 
a matrix element of photoproduction. 

We observe here that Æ7(q, p) does not contribute any more to the re- 
scattering correction in the integral equation, since it is a S-wave in the b 
variable. It can be expressed in a suitable way by means of the scattering 


matrix. 
We have 
Ani 
(13) FE? (; P) = (Ske wo wh 287 Pr a) P, (€, P)hi(@,) ? 
€ Dal 5) Ty 


if x and pf are the isotopic index of the p and q mesons respectively and where 
P,(e, p) = CEC 2p, 
P,(e, p) = 9 €"p—-0'E0'p,; 


(14) 
Pi(y, à) = $4, t 


Moreover, we can neglect the term 


ii FAL (+) 
> (Vi farine) CH VE Pa) 


= + — ©, — te’ 


because of the smallness of such a matrix element in comparison with other 
inhomogeneous terms, and since unless the photon energy is bigger than the 
meson rest mass, the deniminator cannot go to zero. 

To recapitulate, we have: 


(15) FE (Gs P) = Hr (4,P) + FEE (4), 
where € (9, p) is defined by the integral equation 


(16) FE (q, p) = FAG, pi D (vet Woy Fa P) 


Wr — y — de 


Therefore the total radiative pion nucleon Scattering matrix-element is 
given by: 


Fi (G, P) = FE (4, p) + FEL(G, p) + FEM (4, p). 


3654 


THE RADIATIVE PION-NUCLEON SCATTERING IN STATIC THEORY 693 


3. — Solutions of the integral equations. 


We study in this section the solutions of Eq. (7) and (16). Starting wtih (7) 
we observe that the last three terms can be neglected since two of them have 
a denominator which never vanishes, the denominator of the remaining one 

can only go to zero if the photon energy is larger than the meson rest-mass, 
and this case will not be considered as will be discussed in Sect. 4. Using the 
projectors ,P, and ;P, defined in Appendix A we can write: 


a7) ee ae ee 


, 
iT (8040,Mu)? 


where R satisfy the integral e uation: 
ti è =) 


AT sd 1 hi, o R Or, © È 
(18) RB (@., ©) =— + —2— + = (0) Bs (Orr On) (1) dor. 
Va Op — Wa N WO, — Wa — ie 


The matrices AY and Bi are also defined in Appendix A. For the solution 
of such an equation, we use the RODBERG (5) result which in our case reads: 


AT (© ) BU (M ) Og x 
TJ , i i, j DI a0 7 D i hi; 7 
(19) Ei (04) Wp) | À: ' hy;(03) Wp Og — Wo ts (09) 
where for 
1 a 9 CA 
t; = |1 3] 5 hee al oa 
e, ez: 


We now turn our attention to Eq. (16). Since the equation is linear, for the 
sake of convenience we write: 


(20) dei (4, p) == Hei (4, PD) = Fei (4; p) ; 
where pes 
(PV ) 
(20°) Hi = AGP) — LT, = o — ie) FEM; D) è 
(P ee D ) 
(20") HE = Hi (q; P) — om CSP) Hi P) 
de 
(21°) FEY (yp) TI Opa e bapa Per Po) è 


(Ska, Ogre 


SERIA A) 
TA 2 = © apy Ys HE 
(Aka)? © — Oqik + 


(21") ME (ID) = 


(5) See equation (28) of ref. (*). 


fi 
À 


wi 
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Here « and B stay for the isotopic index of p and q ,mesons. Starting with 0) 
in order to deal with integral equation without 6-function it is convenient 
to introduce a new function C;,(q, p) defined as: 


(22) On (4, p) = Hi (14, p) — Hr (4, PD) - 
The resulting equation for C,(q, p) is: 


Die €: ER ot be Se 
3), 04g, D 0 = Se, sel 


(8ko,~,)* @,— my — ie 5 


704+YS(-) 
(Contatta a pr). 


On — Ma — LE 


where r= p—k. 
It can be easily seen that (23) is satisfied by: 


(24) C.(g p, r) = CH(q, p, r) = 


= OP © pce Py, Po) XE]. la) Pg, 7) Cogley 0) 


(80,0,k) ©, 


iJ 
if C,,(@ , ©) is solution of the equation 
1 1/ v(n)n? do, 
20) ICE (0 — ht (@,)+- — | — RE (n) Cra, ©). 
( ) CF Day € >) Or =a) oy — ite 7g à, T x | On eu Wa ea pal | n) rah n°, r 


This equation is of same type as (18); the explicit expression for C,, is: 


i hE (@r) Wa 1 | 

(26) Cibi ani Laz kw). 

Finally considering (20"), we have to observe that the kernel will rescatter 
only the p-waves. Then we can consider in the 3; only that part which con- 
tains p-waves in q variable. 

We separate our amplitude into magnetic dipole and electric quadrupole 
parts. Furthermore, we can label our magnetic and electric amplitudes with 
the quantum number of initial and final pion-nucleon System. The electric 
quadrupole transition of initial j = 4 state can only lead to final j = 3 states. 
All the remaining transitions can lead from one of the four initial states to 
one of the four final states of the pion-nucleon System. 

We can therefore set 


(27) Rg, p) = non p-wave part or Æ?7(q, p) + 


8 0(q) v(p) I 
Zoni DOG Pg k Pye, p) + q-ePy(k, p)] + 


+ Oi, P) aL P;(g, k) Pyle, p) + P;(q, €) P;(k,p)] + ME(q, p)ilP,(q, k/\ e)o-p]+ 
+ Mi:(4, p) 3[P.(k, q)P;(e, P)=2£,(q, e) P,(k, p)] <p |\P.t*P,|a). 
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The indices have the following meaning: the upper ones are the total angular 
momentum (at left) and the total isotopie spin (at right) of initial state; the 
lower indices are in the same order the angular momentum and isospin of 
pion-nucleon system in final state. The Q and the M are the electric quadru- 
pole and the magnetic dipole amplitudes. 

The equations for these amplitudes become: 
Bey una) | Me 0) 


A | Mn — Ma — 1€ 


1 3y2(n) do, 
(28°) Mi;i(4 p= MGu(Q, D; k)h,y(0y) + se ci RFO) M PAU p); 


Wn — Wa — VE 


J=1 J=3 
be 3 154 j=t 
er Re 
The explicit expressions for Go Gy are 
| samy | DI TO 22 — nile — 1) for |æ|< 1, 
Gy= . 
| reg] netti 20 for |æ|> 1, 
vi] Ma 
(29) 3 A 
Premiers arr) 6s for jal 
Go= 
| 4 
leo ia +3) for |e|> 1, 
and 
na == Pn 2. 
eae i È 
2kq 


The solution of equations (28) was given by OMNES ‘(Bie 
(30') Qii(4, P) = Gala py k) hop) C08 dis(0a) + 
+= exp [oulon)]P | AE FUER exp [— 04(0)] dg] exp lidia] ; 
(30”) MPG 2) = similar expression with Gg replaced by AjGx 


and 


(6) R. OMNES: Nuovo Cimento, 8, 316 (1958). 
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4. — Discussion and conclusions. 


The results of our investigations show that the scattering amplitudes for 
cr+p—>p+zx-+y are most important, and therefore could have some chance 
to be experimentally observed, when either one or both mesons are on the 
33 resonance. Our method, which treats one of the two mesons exactly, is 
particularly appropriate for the first situation: the pion which is to be treated 
exactly is the one in the resonant state. 

In principle our method applies equally well when either the initial or the 


final pion is the resonant one. 

However we think the more interesting case is when the initial meson is 
a resonant state: this makes the total energy smaller and therefore the static 
approximation and the one meson approximation more reliable. 

Concerning the shape of the photon spectrum our method is appropriate 
for the study of its intermediate energy part. As already pointed out in the 
introduction, the low energy part can be accounted exactly. 

On the other hand, it is not easy to evaluate the error for k~w. The 
high energy part needs further investigations. 


MCE 


The greater part of this work has been done while one of the authors (B.B.) 
was fellow at Theoretical Group of CERN in Geneva. 

We are deeply indebted to Prof. B. FERRETTI, Dr. S. FUBINI and Dr. R. 
OMNES for very enlightening discussions and useful suggestions. 


APPENDIX A 


In the final state we have to couple three angular momenta (4 of nucleon 
spin, 1 of meson q of the vector u — k/\e) in order to get 4 or è of the 
initial pion nucleon system (the initial meson is D). 

It is easy to verify that if we choose a coupling scheme in which the q meson 
is coupled to the nucleon to give the partial angular momentum j the four 
operators of total angular momentum J and partial j which we shall indicate 
as ;P,(u, q, p) are: 


| iP, = A;Pi(q, k\e)c-p, 


(A.1 
a |P = i[P(g;e)P,(k, p)— P(g, E) P,(e, pl, 


wih tie ~>, 
~ 


2 
where À,— | | 
— 1 


P; are the usual operators as defined by CHEW. 
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The corresponding ones for the isotopic spin, which we shall indicate as Py, 
are: 


(A.2) «Py = <P|P,t*P,|a), 


where ti, = ie. 


In (A.1) we introduce the matrix B; in order to preserve the symmetry 
between angular momentum and isotopic spin. The matrix A7 and B? can 
then be written in the form (all the terms are omitted but the 33 resonant one): 


» Q 
eli AT jis 4 pe 
ell! eae 97 33 ? DCR iF 381 sasso rata 
20 8 4 
1 Fe A: LE — I 192 » 
AT = "OF hx, È e = — ne hs ; AS _ 97 ha , 
10 25 10 
13) = > ARS ree ‘ 5 she * 
An= gola; A= on ha — 3h33; AGS = = hss 
27 a 27 
33 * 
A5s 97 hs ’ 
11 211 13 713 gore 
Bii BY Br BY Ba 0 ’ 
4 2 4 
5 eee FE ASE. ne SI 199. JE} 
By = 3 Mss ; By = ae Kiss 3 Bi = Mss 5 
LI 
9 hs h 
= 33 33 
B® == h Be = 8” 
13 3 38? 33 3 


The remaining elements are given by the relations: 
Fe . BE SR 
Le ca A ) Ba vw Bs È 
The h,; are function of ©, and are defined as: 


hj(@,) = exp [70,;] sin 6,;[prvr(p) . 


R TA SES UNTO m7) 


Si studia nel quadro della teoria statica lo scattering radiativo pione-nucleone 
trattando un mesone esattamente, e l’altro nell’approssimazione di un solo mesone. 
Si ottengono espressioni esplicite per le correzioni di riscattering. 


(*) Traduzione @ cura della Redazione, 


a 45 - Il Nuovo Cimento. 
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An Alternative Method of Quantization: 
the Existénce of Classical Fields. 


J. A. CRAWFORD 


Berkeley Astronomical Department, University of California - Berkeley, Cal. 


(ricevuto il 18 Agosto 1958) 


Summary. — A new method of quantization, applicable to an arbitrary 
physical system, is presented. For systems with a classical analog, the 
introduction of quantized canonical variables proceeds without the elimi- 
nation of the corresponding unquantized variables (although the Liou- 
ville equation is in general modified). These remain as legitimate physical 
observables having the same expectation value as and no greater dis- 
persion than their quantized counterparts. Their (non-negative definite) 
joint distribution in phase space then determines the state of the system. 
In the case of bosons interacting only with external fields or sources, 
the unquantized fields continue to satisfy the deterministic equations of 
motion of the unquantized theory. It seems significant that classical 
relativistic fields with.external sources are not eliminated by the process 
of quantization. All the observable predictions of quantum theory are 
identically reproduced and no new results are obtained. This method 
of quantization is therefore an alternative to the usual one. A natural 
extension of the theory yields a deterministic model for direct boson- 
boson interactions (e.g. pion electrodynamics), free of mathematical diver- 
gences. This will be developed in a following paper. 


1. — Introduction. 


While a future physical theory will have to yield the verified predictions 
of quantum mechanics, it may present these in a new mathematical form. 
This possibility makes it worthwhile to search for alternative methods of quan- 
tization. 

The aim of this paper is to introduce a new method of quantization which 
is an alternative to the usual one. Both methods are of comparable mathe- 
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matical simplicity and lead to the same observable predictions. They are not 
however simply equivalent in the sense that, say, matrix mechanics is equi- 
valent to wave mechanics. The new method leads, for instance, to the intro- 
duction of additional physical observables, in particular classical canonical 
variables for those systems possessing a Classical analog. 

The usual method of canonical quantization eliminates the canonical va- 
riables of a classical theory by replacing them with observables which satisfy 
uncertainty relations. As a result, the quantized theory is indeterministic. 
This situation has not unnaturally led to the impression that the existence 
of complementary (i.e. not simultaneously measurable) observables, by itself 
contradicts the possibility of a deterministic physical theory (we shall call a 
theory deterministic when a) the statistical state of the system it describes 
can be completely specified by the joint distribution of certain parameters 
and b) these parameters satisfy an equation of motion of first order in the time. 
An example is classical statistical mechanics; the parameters in this case are 
the canonical variables). 

Boum (1) has proposed a deterministic re-interpretation of quantum theory 
in which, accordingly, complementary observables play no part. Bohm’s theory 
is therefore classical in the usual sense. Quantum effects are accounted for 
by special «quantum forces », of random character, introduced for this pur- 
pose. As we shall see, however, it is not in fact necessary in order that a 
physical theory be deterministic, that it should exclude complementary obser- 
vables. Thus, even if one were to adopt the arbitrary philosophical position 
that the ultimate goal of physics is the construction of a completely objective 
model of reality (2), there would be no reason for ruling out complementary 
observables. 

Much confusion has grown around this question because of the difficulty 
of interpreting a theorem due to VON NEUMANN (3), concerning the possible 
role of hidden parameters in physics. Contrary to what seems to be generally 
assumed, this theorem does not disprove the possibility of the deterministic 
description of a physical system possessing complementary observables. 

The theorem assumes in effect that a physical system must indeed own 
some complementary observables and goes on to draw the conclusion that if 
hidden parameters exist, not all observables can be functions of these para- 
meters: since, even if the values of all hidden parameters are given, not all 
observables can be simultaneously known. However the theorem does not 


(1) D. Boum: Phys. Rev., 85. 166 (1952); D. Boum and J. P. ViGier: Phys. Rev., 


) 
96, 208 (1954). gt i ; ; ? 
(2) P. A. SomrpP: Albert Einstein: Philosopher-Scientist (New York, 1951), p. 675: 
(3) J. V. NEUMANN: Mathematische Grundlagen der Quantenmechanik (New York, 


1943), pp. 170-173. 
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prove (it does not actually claim to) that the state (i.e. the total information 
required to determine the distribution of every observable) of the system cannot 
be completely specified by the values of the hidden parameters or, in general, 
by their joint distribution (*). : 

There exist in fact several familiar examples of deterministic theories which 
are also capable of describing complementary observables. To consider only 
one: the pitch v and the epoch t of a musical note may be classed as comple- 
mentary acoustical observables since they satisfy the relation (°) 


AvAt > 1/42. 


And yet acoustical theory is deterministic. 

Thus, contrary to a commonly held impression, a deterministic model for 
a physical system is not a priori inconsistent with Bour’s fundamental idea 
of a complementary description of physical observables. While this conclusion 
in no way suggests that a convincingly simple deterministic model of all phy- 
sical phenomena can be found, it does imply that the physical world may be 
indeterministic only in part. This is certainly possible at least if one is willing 
to describe the world as composed of separate parts, i.e. if he does not insist 
on a unified theory. If we agree to adopt this viewpoint, at least temporarily, 
we are then naturally led to reconsider the reasons for giving up the « success- 
ful» (relativistic) classical models, namely the classical field models of gravi- 
tation and radiation in the presence of external sources. 

The fact that electromagnetic. energy is quantized implies that the conti- 
nuous classical expression 3 dt(e2--h?) cannot be interpreted as energy. It 
was essentially the difficulty of finding a quantized expression for electro- 
magnetic energy within a theoretical framework which would also include the 
classical field that led PLANCK finally to abandon his attempt to preserve this 
field and Drrac to replace it by its quantized counterpart (°). Although this 
method of quantization was completely suecessful, the foregoing considerations 
make it worthwhile to reconsider the question of the necessity of giving up 
the classical field. The problem may be stated as follows: «Can an alternative 
method of quantization be found which does not destroy the classical field? » 


If this is possible, the quantized and classical fields may clearly not be iden- 


tified with each other: for instance, the components of the former are joined 
by the Heisenberg uncertainty relations whereas those of the latter are not, 
Once this necessary distinction between the two fields is realized, certain para- 
doxes are more easily avoided. 


(4) F. Bopp: Ann. Inst. Henri Poincaré, 15, no. 2, 81 (1956). 

(5) A», At are not in this case statistical dispersions; they simply give the accu acy 
to which » and ¢ are defined. 

(5) P. A. M. Dirac: Proc. Roy. Soc., 114, 2434, 7100027) 
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In the following Section a new general method of quantization is presented. 
When applied to a system with a classical analog, this method introduces 
quantized canonical variables without however eliminating the corresponding 
unquantized variables: their (non-negative definite) distribution in phase space 
then determines the state of the system. The physical meaning of these clas- 
sical (unquantized) variables is made clear by the proof that they always have 
the same expectation value as, and no greater dispersion than, their quantized 
counterparts. The measurement of a quantized variable thus produces a value 
for the corresponding classical observable. 

Classical boson fields in particular are retained as legitimate physical obser- 
vables. When bosons interact only with external fields or sources, the clas- 
sical boson fields are shown to satisfy the deterministic equations of motion 
ot the unquantized theory. It is worth noting especially that the classical 
electromagnetic and gravitational (and also mesic) fields with external sources 
are not eliminated by the new process of quantization (°7). All the predictions 
of quantum theory are identically verified and no new results are obtained. 


2. — Quantization. 


Our starting point is the Hilbert space of quantum theory. 
2°1. States. — We postulate that the state of a system is represented by a 
time-dependent self-adjoint linear operator ®. This D is not the statistical 
operator U of von NEUMANN (5), but we shall find later that 
as D2, 
If M, N are two linear operators, we define their scalar product (M, N) by 
the relation 
(M,N) = Tr MN). 
It follows that, if (M, N) exists, 
(M, N)= (N, M), 
and if (M, M) exists (M is then said to be normalizable), 
(M, M) > 0. 
(7) S. N. Gupra: Rev. Mod. Phys., 29, 334 (1957). 


(8) J. V. NEUMANN: Mathematische Grundlagen der Quantenmechanik (New York, 
1943) p. 168. 
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Normalizable linear operators form a manifold which is itself a Hilbert space. 
We shall refer to it as the associated Hilbert space (AHS). 


2°2. Observables. — Physical observables will be represented by self-adjoint 
linear superoperators defined as follows: a superoperator is àn operator in AHS. 
The theory of operators in Hilbert space (*) is therefore directly applicable 
to superoperators. For example, a superoperator s transforms every linear 
operator M of its domain in AHS into some normalizable linear operator M': 


sM = M'. 


Superoperators s, s* are said to be adjoint if they have the same domain, and 
in that domain 


(st M, N) = (M, 8). 


The linear character of s means that its domain is a linear manifold and, if 
a, …, d are complex numbers, 


s(aM +... + dP) = asM +... + dsP. 


2°3. Expectation values. — It is postulated that the expectation value of 
the observable s in the state ® is (1) 


<8) = (®, 59). 
This requires 


(1) (Op) ans, 


It is further assumed that the observable f(s), where f(s) is a real power series 
in x, is represented by the superoperator f(s), and the observable s+¢ by the 
superoperator s+. 

If @ is multiplied by an arbitrary phase factor, expectation values are 
evidently unchanged. However the requirement that ® is self-adjoint restricts 
this factor to +1. Thus +@ and — © represent the same physical state. 


2°4. Special observables. 


a) Twin quantized observables. — Let A be a self-adjoint linear 
operator in Hilbert space (A — At). Then the superoperator a’ will be de- 


(9) J. V. Neumann: Mathematische Grundlagen der Quantenmechanik (New York, 
1943) p. 47-48. 


(19) For convenience, we use the same symbol s for the observable as for its 
associated superoperator. 
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fined by (1) 


(2) Col =A, for every M. 
a! is obviously linear, and we have 

(@M, N) = Tr ((AM)'N) = Tr (M'AN) = (M, aN). 
Thus a! is self-adjoint. We find 


(3) <f(a')) = (®, {(a)D) = Tr (8 (A)®) = Tr (Uf(A)), 
where 

Ca 
We may also define (1) 


(4) at M = MA, for every M. 
a’ is also linear and self-adjoint, and we find 
(5) <f(a") = Tr (U f(A)) = <f{a')» . 


Thus the observables represented by a! and a’ are indistinguishable, in 
essentially the same sense already familiar from the quantum theory of iden- 
tical particles (e.g. the co-ordinates 7, 2, of two electrons are indistinguishable 
observables in quantum theory). 

Tf U is now identfied with the statistical operator of von Neumann, the 
expression 

Tr (Uj(A)) 


is the expectation value of the quantum-mechanical observable f(A). It is clear 
then that, according to the theory proposed here, what is called in quantum 
mechanics the observable A actually represents two indistinguishable « twin » 
observables a! and a”. This doubling of all quantum observables suggests that 
this method of quantization might be called double quantization. 

When a! (or a’) is precisely known, 4.6. 


Cal) — <a» =0, 
we have, from (3) 


Tr (UVA?) — [Tr (UA)F = 0. 


(11) This definition will hold also when A is not self-adjoint. 
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This is well known to imply that Tr (UA) is an eigenvalue of the operator A. 
Thus, from (3), when a! is precisely known, its value is an eigenvalue of A. 


b) Classical observables. — It is of interest to examine the pro- 
perties of the observable ae 


(6)  @2=4(a'+ a‘). 
A real function f(a) of such an observable, and also a sum 
f(a) + g(b) +... 
of such functions, may conveniently be, called a classical observable. Classical 
observables do not necessarily commute with each other, but we shall see in 
| the next Section that classical canonical variables, defined according to (6), 
do mutually commute. 
It follows from (3), (5) that 
<a> = <aly= <at> = Tr (UA). 
Consider now the dispersion of a. We have 
<(Aa)®) = <a?) — Ka)? 
= H<a") + Kata”) + Kara!) + <a>] — <a? 
= Ka) + Kala”) + Kara) — Cab. 


. It is readily verified that «1» observables commute with «r» observables: 


[a', b']M = A(MB)— (AM)B =0, for every M. 


[ai, bh] = 0. 


In particular [a!, a] == 0. We can then define a joint Dobe bite distribution 
P(a', a’) as follows: 


Ple, y) = <6(a'— x) d(a"— y), 


dla' — x) = (1/27) fas = [iz(a'— 2)], ete. , 
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and exp[iz(a'—x)] may be written as a power series in e(a—x). It is 
easily proved that this distribution is non-negative definite. We can then 
use the Schwarz inequality 

3/23 


: = FE N PEER ae 
Lea = tae = <a" ><a 3° = Ka! >. 


We then obtain 
<(Aa)?> < Ca) — <at> = <(Aat)*) = <(Aa”)?) . 


Thus a has the same expectation value as, and no greater dispersion than, 
a (or a"). The measurement of the « quantized » observable a' (or a") (which 
in quantum theory is interpreted as a measurement of A) produces therefore 
a value of the classical observable a with at least as good accuracy. 
It can also be shown that the uncertainty relations connecting any two 
classical observables o, + take the trivial form 
<(Ao)?>#<(Arz)2>? > 0. 
2°5. Equation of motion. — The mathematical apparatus presented above 
closely resembles that of quantum theory, the state operator ® playing the 
role of the pure state vector in quantum mechanics. Pursuing this analogy 
we write down the equation of motion as follows 


(7) ih 0D|/ct = k®, 

where Æ is a self-adjoint superoperator which we shall call the superhamil- 
tonian of the system. This guarantees that the normalization condition (1) 
is conserved. % must also be so chosen that ® remains self-adjoint. The time- 
derivative § of an observable s may then be defined by the equation 


ihs = th os/et +[s, k]. 


In this paper we shall limit ourselves to a choice of k which guarantees that 
the results of quantum theory are identically reproduced. We choose 


(8) k—=h—k, 
H being the quantum-mechanical Hamiltonian ordinarily associated with the 


system. Thus k does not represent energy. To see that this choice of k does 
in fact lead to the results of quantum theory, we note that (7) may be written 


ih e@jot = HO— OH, 
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which conserves automatically the self-adjoint character of D. Von Neumann’s 
equation of motion ('*) 


ih cU/ct = HU — UH 


immediately follows. 
Later (Sect. 6) we shall have occasion to consider the possibility of choosing 


k differently. 


8. — Classical canonical variables. 
Consider the linear operators Q, P satisfying 

[OE 

It follows readily from (2) and (4) that 
[a,pl= th, 
[g,p]=—h. 

If we define the classical canonical variables according to (6): 
debe), 


D= Hip Ep), 
we find 
[p,q] = 0. 


In general, if we have several linear operators, 7.e. 
[Q:, Pi] = hous, 


we find that the set of classical variables q, p mutually commute. 
If we define the linear expressions 


| A=YcQ+4P,, 
| B= Dei +hPs, 


(c, d, e, f complex numbers), it follows that 


(10) [a,b] =0. 


(12) J. V. NEUMANN: Mathematische Grundlagen der Quantenmechanik (New York 
1943) p. 186. 
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4. — Quantum electrodynamics. 


We shall now apply these considerations to quantum electrodynamics. The 
Maxwell potentials A, in quantum theory may be assumed to satisfy the 
equations 
(11) LIA, =—-(1/0)J,. 


it 


It follows that 


and therefore 


The usual quantized field components 


F,, = 0A,/ox, — 0A,,/ 0x, 


i 
are expressions of the type (9), and therefore the classical Maxwell field com- 
ponents 
fay = 04,102, — Ca, / Cx, 
mutually commute, according to (10). According to the discussion of classical 


observables in Sect. 2°4, a measurement of F’,,, produces a value for f,,. 
If we neglect the sources, the equation 


shows clearly that the propagation of the classical potentials is deterministic. 
The same is true if the sources are external, so that Me j, and j, are replaced 
by the same explicitly given function of the space-time co-ordinates. 

The supplementary condition needed to insure that on fi ODEY Maxwell’s 
equations may be written, according to Fermi’s method (1) 


xD = Dy =—0, 
7D = DY = 0; y = 0A,/0®,. 


(13) G. WENTZEL: Quantum theory of fields (New York, 1949), p. 112. 
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These conditions are conserved by the equation (11), and 


5. — Boson fields. 


If we neglect all interactions except those with external sources or fields, 
the quantum-mechanical Hamiltonian H(Q, P) of a boson field is a quadratic 
function of the Q,, P;, the coefficients being possibly time-dependent (see 
equation (19)). The canonical equations of motion 


Q, = 0H(Q, P)]oP;;. “P,=— 0H(Q, P)/00, 


are then linear. It follows from (2), (4) and (6) that the «/» and «7» cano- 
nical variables as well as the classical canonical variables satisfy the same 
equation. We have 


(12) di = CH (q, p)/0p: ; D: = — OH (q, p)/ eq: . 


Thus the classical canonical variables obey the deterministic equations of 
motion of the unquantized theory. It follows that the probability distribution 
9(q, p; t) in classical phase space must satisfy the Liouville equation (4). In 
the following Section we shall construct an expression for g(q, p; t) for an ar- 
bitrary physical system and verify the Liouville equation for the case dis- 
cussed in the present Section. We will at the same time show, for this case, 
that the theory of bosons is deterministic and furnishes a model identical 
with the one provided by the unquantized theory. This will enable us to for- 
mulate a suggestion for a new treatment of direct boson-boson interactions. 


6. — General canonical representation. 


The state. operator ® of an arbitrary physical system may be represented 
as a matrix if an orthogonal basis in Hilbert space is chosen. More generally 
we shall consider a Dirac representation based upon the usual canonical co- 
ordinates Q; possessing a classical analog and other variables N, without clas- 
sical analog (for instance fermion occupation numbers). If a boson field y 


(14) For charged boson fields, the canonical variables as usually defined are complex. 
They may however be replaced by real variables (cfr. beginning of the following Section) 


. 


AN ALTERNATIVE METHOD OF QUANTIZATION: THE EXISTENCE OF CLASSICAL FIELDS 709 


is complex it may be decomposed into real fields y,, y, according to the 
relations (1) 

y = (1/V2)(y1+ ivo) 

y' na (1/V2) (y: — ty) 
and we may therefore assume without loss of generality that the Q; are real. 
The «matrix element» of ® is then 
(13) <Q'N'|D|QO"N"} 


in the usual Dirac notation. 

We shall not for the moment make use of the fact that @ is self-adjoint, 
and for brevity we omit reference to N', N" in the matrix element, writing 
<Q'|®\Q"> in place of (13). 

In order to represent the superoperators qj, dj; Pi, p; we consider for 
instance the matrix element 


<Q! \q'B|Q") = <Q" |Q.|Q"> =|d9"K'12s10"><Q" | BIQ"> = @;c0'|D|0"). 
Similarly we find 
<Q! |g, ®|Q"> = 00" |D]Q") - 
If the operator P; in this representation is given by 


ba = (hi) 0/00; 
we obtain 
<Q’ (pid |Q" > = (fili) o<Q'|D|Q" » EQ, i 


Likewise, 


| 


<Q' |pi®|Q"> =[agrca'ian<on| P10” = |dQ”"<Q" | P,|Q">* <Q" | P19") = 
= — (h/t) 0Q' | B|Q">/2Q; . 


fr l an N Fr ri " 
We see that the superoperators Qi, dr Pr Dj are represented by the following 
1 Che f/f)! ; i 
hermitian functional operators acting upon the matrix element <Q'|?|Q"> 


(15) G. WENTZEL: Quantum theory of fields (New York, 1949), p. 12. 
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considered as a function of the variables Q,, 9%, t: 


Rada (1), = 9 
te | . An! ex PANE 
| (Pap = (Ji) 0/00, ; (Dios pS: (Ai) 0/00; . 


For instance <Q'|p'®|Q"> is a function of 0°, 0°, t which is obtained by ope- 
rating upon <Q'|®|Q"> with the hermitian operator (#/î) 0/ oq). 

Starting from this representation we may then make use of the technique 
of transformation theory to obtain other representations. We consider first 
a simple change of variables: 


| 
(15) | 


We write 
<Q'| BQ") = x(q‘, k';t). 

From (15) follows 
ae (Ji) 0/0Q, = (h/t) 0/ek; + (k]2i) 2/04, , 

— (hfi) 0/6Q; = (Ji) d/dk, — (h/2t) 0/04, . 
Reealling (14), the relations (15) and (16) give us the hermitian operators re- 
presenting q°, 9", pî, p' in the q’, k’ representation. We can now very simply 
obtain a canonical representation by performing the Fourier transformation 


(17) P(q', p'; t) ', k'3 t) exp 


in rap 


This yields 


piP(q', p's t) (fi 


uu fai =]. pri 


[(R/i) Ox (q', k’; t)/0k;] exp 


in! sai 


c;x(g', k'; t)] exp È AE La 


à i ‘ i © pa o 1 7 l » 
Thus, in the canonical representation, q;, 9; p;, Pj are given by 


Doe aed, ih |2) 0] Op., ; r Za his : 
(18) | (don = 9; + (0h/2) 0/0p, (Dh. = 4 — (14/2) 0/0p, , 


(Pion = Pi (4/2) 0/0g, 3 (ph, = pi + (64/2) 2/04, . 

It follows in particular that 
(Gon cati q, ’ 
(Plon = Pi 
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The probability distribution function g(q', p'; t) in the classical phase space 
may Dow be obtained. Let f(q, p) be a function of the canonical variables 
expressible as a power series. We obtain 


Ha, p> = | | TI 40:40; <Q'N"|8|Q'N">*-<Q'N |, 9) DION") = 


= 2 > -{T dg dp’ f(q', p')}\P(a', p's N', N E 
-|--|T] dgidp! gly’, p's Did) - 


where we have written in the explicit dependence on N,, N°. 
s 


Thus the (non-negative definite) classical phase space distribution is (drop- 
ping primes) 


g(g, Bs 1) = X DM pi NN"). 


Ni N; 


It may readily be verified that the requirement that is self-adjoint, 
N09 = NOON, 
is equivalent to the condition 
Pg, p; N°, N"; t) = P*(q, ps N", N;1). 


In particular, if the physical system has a classical analog, the variables N SUN 2 
do not occur and ¥ is real. We have then 


qa Pr 


If in that case g is a given analytic function of the classical canonical variables 
diy Pi, its Square root Ÿ is given by 


wW—=- +Vg, analytic, 


i.e. is determined except for the sign. Inversion of the transformations (17) 
and (15) then determines ® except for a factor + 1. We have seen however 
(Section 2°3) that two @’s differing in sign determine the same physical state. 
Thus, whenever the system has a classical analog, the classical phase space 
distribution (assumed to be analytic) determines the state of the system. When, 
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in addition, the q;, p; satisfy first order equations of motion, as in (12), the 
theory of the system is deterministic. 

In this case (bosons interacting only with external sources or fields), the 
quantum-mechanical Hamiltonian is a quadratic expression 


(19) H(Q, P) = > dis PP; + b:5Q:Q; +36(P;0; + Q;P,) à + a;P,+ b;Q; + €. 


The equation of motion for Ÿ may then be obtained by writing (7) in the 
canonical representation. It is readily verified that 
h'= H(q', p') ; k° = H(q', p’)- 
Then 


k = H(q', p') — H(q’, p’) . 


If the operators (18) are substituted for q', qj, p°, p! in this expression, 
we find after some cancellations: 


À r ! La! cH Tr ! y 
(20) By = Ma Bh OE RIE 
2 cq; CP; Cp; Gi 


The equation of motion 


ih OV jot = kV 


0? 
is then simply the Liouville equation for Y. Since it is linear in the diffe- 
rential operators €/0q,, 0/0p;, 0/0, it follows that g — Y2 also satisfies this 
equation (15), as is required, since in this case the classical canonical variables 
obey the classical equations of motion (12). 
| The canonical representation is of some interest as it suggests a possible 
‘ay of modifying the physical theory. ; 

The fact that the deterministic equations (12) are satisfied when bosons 
interact only with external fields or sources, strongly suggests that the theory 
should be so modified that they remain satisfied also when bosons interact 
directly with each other. Y would then, also in that case, obey the Liouville 
equation (17). If this were indeed the case, Ÿ would remain real, and ® there- 
fore self-adjoint as required, Also the superhamiltonian would have the form (20) 
which, as one can easily verify, guarantees that it is self-adjoint. However k 


(9) See G. A. BAKER jr.: Phys. Rev., 109, 2198 (1958), Appendix. 
(7) For instance in the electrodynamics of charged pions (the electrodynamics of 
‘deuterons is of course not a direct boson-boson interaction). 
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would no longer separate the «7» and «7» variables as in (8), and departures 
from the predictions of quantum theory would be met. 

This program will be carried out for charged boson electrodynamics in a 
following paper. This deterministic model is also interesting since, being a 
classical field model, it is manifestly free of mathematical divergences. 


* KOK 
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RIASSUNTO (*) 


Si presenta un nuovo metodo di quantizzazione applicabile ad un sistema fisico 
arbitrario. Per sistemi aventi un analogo classico l’introduzione di variabili canoniche 
quantizzate si esegue senza eliminare le corrispondenti variabili non quantizzate (benchè 
l’equazione di Liouville sia in genere modificata), Queste restano come osservabili fisiche 
legittime avendo lo stesso valore probabile e dispersione non maggiore dei loro equi- 
valenti quantizzati. La loro distribuzione (non negativa, definita) nello spazio delle 
fasi determina allora lo stato del sistema. Nel caso di bosoni interagenti solo con campi 
o sorgenti esterni i campi non quantizzati continuano a soddisfare le equazioni deter- 
ministiche di moto della teoria non quantizzata. Appare significativo che campi rela- 
tivistici con sorgenti esterne non siano eliminati dal processo di quantizzazione. Tutte 
le predizioni osservabili della teoria quantistica si riproducono identicamente e non 
si ottengono nuovi risultati. Il presente metodo di quantizzazione è pertanto un’alter- 
nativa per quello usuale. Una naturale estensione della teoria fornisce un modello deter- 
ministico per le interazioni dirette bosone-bosone (p.e. l’elettrodinamica dei pioni) 
esente da divergenze matematiche. Questo sarà sviluppato in un prossimo lavoro. 


(*) Traduzione a cura della Redazione. 
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Magnetic Perturbation of Cathode Rays. 


S. YAMAGUCHI 


Institute of Physical and Chemical» Research - Tokyo 


(ricevuto il 29 Agosto 1958) 


summary. — When an electron beam penetrates a magnetized body, 
splitting of the electron beam occurs. This phenomenon is explained in 
terms of polarization of the cathode rays. 


1. — Introduction. 


An electron beam is perturbed by magnetic fields, and spreading and de- 
flection of the beam takes place here. The behavior of an electron beam that 
penetrates a magnet was studied in the present work. Splitting of the beam, 
distinguishable from the above-mentioned phenomena, was. observed after the 
beam had tunneled through the magnet body. The author has elucidated this 
phenomenon as the polarization of free electrons. 


2. — Experimental. 


The present experiment was carried out with an electron diffraction ap- 
paratus. A single crystal of iron (length: 10 mm, diameter: 1 mm) was edged 
mechanically and then the edge thus prepared 
was etched with hydrochloric acid. One piece of 
permanent magnet (coercive force, about 700 Oe) 
=Fe was attracted to the rear side of the test piece, 
Single crystal à illustrated in Fig. 1. In this way the single 

crystal was magnetized. An electron beam (wave- 

7 length, 0.0308 A) grazed the edge of the single 

Uj Mosne crystal. The diffraction pattern so obtained is 

Fig. 1. - Arrangement of shown in Fig. 2. In Fig. 2 the incident beam 
experiment. runs nearly perpendicular to the (111) plane of 


Electrons 
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Fig. 2. — Diffraction pattern obtained from Fig. 3. — The central spot found in 

the magnetized single crystal of iron. Fig. 2, enlarged optically 60 times. 

Wavelength, 0.0308 A. Camera length, One sees splitting of the electron 
495mm. Positive enlarged 2.3 times. beam. 


the iron crystal. The central spot found in Fig. 2 was enlarged optically 
60 times, as shown in Fig. 3. It should be noticed in Fig. 3 that the electron 
beam is split into discrete beams. 
Fig. 4 shows the incident beam which 
grazed the non-magnetized test piece. 
There is no splitting of the beam. 


Fig. 4. — The electron beam which grazed 
the non-magnetized specimen. 60 times 
enlarged. There is no splitting. 


3. — Discussion I. 


When the electrons pass through a magnetized body, interaction occurs 
between the spins of the incident electrons and the steep gradient of the 
magnetic field. This gradient of the field exists at the border between two 
Weiss magnetic domains and the Bloch magnetic wall. In Fig. 5, the relation 

between the direction of incidence and the field is 


agit illustrated. The deflection AZ of the spinning 

electrons by the magnetic field is given by the 
& : Fa equation 

ta GE | (A) 

Z=B\|—\m LI LS 1), 

7 /\ gq) A2-8 inf), >), 
L | 

Fig. 5. — Relation between the incident electrons and 

the magnetic field. L, camera length. 1, dimension of 

Screen : ES 
le-A 2 ke A 2 >| Weiss domain. AZ, observed splitting. 
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where B: the Bohr magneton (0.93-107?° erg/gauss), 0H/0Z: the gradient of 
the magnetic field, m: electron mass (9.1-10-28 g), À: the wavelength of the 
electrons (0.0308 A), J: the distance travelled by the electrons in the region 
of magnetic field gradient, L: the camera length (495 mm), and h: Planck’s 
constant (6.6:10-?’ erg:s). 

The spacing measured between the split spots in Fig. 3 (AZ) is 0.01 cm. 
The linear dimension of the Weiss domain (10-3? cm) was employed as the 
value of | in Eq. (1). Therefore, from Eq. (1) we obtain 


It is expected that there is a steep gradient of the magnetic field at the border 
between the Weiss domains and the Bloch wall, especially at the surface of 
the magnet body. It is known that the thickness (AC) of the Bloch wall is 
about 100 A. If we assume 


| AR 
+ NO 
we obtain 
AH ~ 10° G, 


where AH means the strength of the Weiss magnetic field. This observed 
value is approximately equal to the known value of the Weiss field. 

Here the author has considered the Weiss field as a real magnetic field. 
This treatment is not quantal, but classical. This is quite similar to Bragg’s 
relation regarding net planes in electron diffraction, which regards interplanar 
spacing as real. In quantal discussion, the present electron polarization should 
be treated as the effect of the exchange force between the spins of the incident 
electrons and the oriented spins in the magnet. 


4. — Discussion II. 


According to the uncertainty principle, it is impossible to polarize a beam 
of free electrons into two beams of separate spins by means of the Stern-Gerlach 
process (1). Spreading of a flux of electrons by a magnetic field with gradient 
masks completely the splitting of the electron spins by the same field. The 


(1) N. F. Morr: Proc. Roy. Soc., A 124, 440 (1929); N. F. Morr and H. S. W. MASSEY : 


The theory of atomic collisions, 2nd ed. (Oxford, 1950), p. 61; E. van DER SPUY: Nuovo 
Oimento, 4, 1349 (1956). 
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force required to split the spinning electrons (F,) is negligible as compared 
with that to spread the electron beam (F;). That is 


F, > IBS . 
This is because 
h 
Ao:0>-, 
m 
or 
Az À, 


where v is the velocity of the free electrons, Ax is the cross-section of the 
electron flux, and 2 is the wavelength of the free electrons. It is impossible 
to observe the electron spins in the Stern-Gerlach process, insofar as Ax, that 
is, a slit for the electron beam, is not comparable with the wavelength applied. 
Here, the intensity of the electron beam passing through this narrow slit is 
not only too weak to be observed, but also it is perturbed physical-optically 
at the slit. 

In the present work, it is shown that there is a natural slit of the order 
of the de Broglie wavelength, whose periphery does not screen the incident 
beam. This slit does never perturb the incident beam physical-optically. There 
is a very steep gradient of the magnetic field at this slit, i.e. at the border 
between the Weiss domains and the Bloch wall. On grazing the surface of 
a magnet body, the electron beam is perturbed by the steep gradient of the field. 


RIASSUNTO (*) 


Quando un elettrone penetra in un corpo magnetizzato si verifica una separazione 
del fascio di elettroni. Si spiega tale fenomeno in termini della polarizzazione dei 
raggi catodici. 


(*) Traduzione a cura della Redazione. 
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The Use of Methyl Iodide-Propane Mixtures 
and Tri-Fluoro Bromo Methane as Bubble Chamber Fluids. 


F. W. Buttock, C. Dopp and G. E. KALMUS 
Physics Department, University College - London 


(ricevuto il 1° Agosto 1958) 


Summary. — The suitability of a methyl iodide-propane system as à 
«heavy liquid » bubble chamber medium has been investigated. Operating 
conditions have been found for various mixtures, and also for tri-fluoro 
bromo methane, another possible heavy liquid medium. 


1. — Introduction. 


. Attempts have been made in various laboratories (**) to find a bubble 
chamber medium having high density and atomic number, and therefore good 
+ conversion properties. One of the most promising liquids from the point 
of view of radiation length alone, would appear to be CH,I, but unfortunately 
iodine is liberated from this material under the influence of light and heat, 
and the liquid rapidly becomes photographically unsuitable. Several substances 
have been suggested as possible inhibitors of the disintegration, the most 
popular of which is propane since it is itself a well-tried bubble chamber liquid 
of fairly low operating temperature, and therefore may be expected to lower 
the inconveniently high operating temperature (210 °C) of pure CH,I. 


2. — Chamber. 


We have constructed a small chamber 14in. dia. X1 in. entirely of monel 
metal in which heavy liquids can be investigated. Expansion is by means 


(1) J. L. Brown, D. A. GraseRr and M. L. PERL: Phys. Rev., 102, 586 (1956). 


(2) E. D. ALYFA jr., L. R. GALLAGHER, J. H. Murrins and J. M. TEEM: Nuovo 
Cimento, 6, 1480 (1957). 


(3) I. PLESS: private communication. 
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of a fluon diaphragm backed by rubber, placed at the end of a 3’in. dia. x8 in. 
pipe along which a temperature gradient is maintained. The chamber tempe- 
rature is controlled by closely wrapped electrical heating tape, and measured 
by thermocouples placed in the metal body + in. from the top and bottom 
of the liquid. 


3. — Methyl iodide-propane system. 


Various percentages by volume of methyl iodide in propane were intro- 
duced into the chamber by distillation and the operating range of sensitivity 
determined. A 10 millicurie Co source was used for most of the work but 
in a few cases gammas from a 
29 MeV microtron provided the 
ionizing particles. The results 600 
are shown in Fig. 1 as a plot of 
operating temperature against 550 
equilibrium pressure, the solid 
lines representing the limits of 
sensitivity for the system. 

Propane was found to be a 
successful inhibitor only in large 


— 


Ibs/sgr in(147 Ibs/sqr=\ atm) 
gauge ) 


concentrations. For instance a dr 

50%-50% mixture by volume —400-2 & SI 
remained completely colourless g ; fi 
for only a short period once the 350 È si 


operating temperature (115 °C) 
was reached. Continuous oper- 300 
ation was possible only for about 
an hour by which time the liquid 
was pink enough seriously to 
limit the photographic repro- 
duction. However, after subse- 
quently leaving the liquid to 
cool to room temperature it was 
found to clear, the free iodine 
being deposited as a film around the chamber walls. Higher CHI concen- 
trations could be operated for correspondingly shorter times, whilst the pure 
liquid was already a deep red colour by the time operating temperature had 
been attained. 


Lure 


4. — Tri-fluoro bromo methane (CF;Br). 


Because of their low critical temperatures and pressures certain freons 
offer attractive possibilities as bubble chamber liquids (*). Further, it is pos- 
sible to obtain brominated freons, which, whilst retaining the convenient oper- 


(4) G. G. Brinov, Iv. S. KRESTNIKOV, M. A. Lomanov and Ia. IA. SHALAMOV: 
Zu. Exp. Teor. Phys., 5, 1281 (1957). 
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ating conditions also have fairly short radiation lengths (5). We have inves- 
tigated in particular tri-fluoro bromo methane (trade name Isceon B.1) and 
have determined its critical temperature and pressure to be (67+1) °C and 
(595+ 8) Ibs/sq. in. respectively. Sensitivity of this liquid to ionizing radiations 
was found to commence at 22 °C whilst the foam limit occured at 34 °C. 


5. — Conclusion. 


The radiation length of a 50%-50% by volume mixture of CH,I and C,H, 
at its operating temperature is about 10 cm. However, a 25% CH,I concen- 
tration, which is almost the highest concentration that can be operated for 
any reasonable length of time without obvious discolouration, has a radiation 
length of about 17 em. This compares with a value of 11.0 cm for CF,Br 
at 30°C. The latter substance is completely stable at its operating temper- 
ature and prolonged exposure to light has produced no observable discolour- 
ation. In addition, it is non-toxic, non-inflammable, and non-corrosive. AS a 
bubble chamber liquid it is therefore ideal, and would seem to be much pref- 
erable to the CH,I—C,H, system which, besides its less convenient operating 
conditions, also presents difficulties in mixing and in retaining uniform con- 
centration, disadvantages inherent in gas-liquid systems. 


(5) D. V. Buca: private communication. 


RIASSUNTO (*) 


Si é esaminata la convenienza di un sistema ioduro di metile-propano come riem- 
pitivo di una camera a bolle a «liquido pesante ». Si sono trovate le condizioni di 
funzionamento per diverse miscele ed anche per il trifluoro bromo metano, altro pos- 
sibile liquido pesante di riempimento. 


(*) Traduzione a cura della Redazione. 
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Natuurkundig Laboratorium, Rijksuniversiteit Gent - Gent 


(ricevuto l’11 Agosto 1958) 


Summary. — Phenomena, due to the multipactor effect and observed in 
the linear electron accelerator of the Ghent University are reported. An 
explanation is proposed for the variation of the anomalous attenuation 
due to the multipactor effect as a function of the R.F. power fed in 
the accelerator. 


While operating the linear electron accelerator at the Ghent University 
the following interesting phenomena, were observed. 


1) As soon as the magnetron is switched on and R.F. power is fed into 
the machine, the gas pressure in the apparatus increases by a factor of about 
three. Then it decreases gradually and returns to its original value after 
several tens of minutes. Often discharges occurred in some places in the 
accelerator and waveguide structures, especially at the input « door-knob » 
coupling. If the R.F. power is turned off for a few minutes and then turned 
on again, the highest value of the pressure is less than before and returns 
faster to its original value. Finally after about ten cycles of turning the R.F. 
power on and off, the pressure is no longer influenced by the R.F. power. 


2) If the focusing magnetic field H, is switched on while the magnetron 
is running, the pressure rises again. A variation of the current in the focusing 
coils changes the gas pressure in the accelerator tube. The time needed to 
return to a pressure slightly higher than the original pressure, decreases for 
every next run, but never becomes zero. 


3) H, causes an anomalous attenuation of the R.F. power if this is at 
the normally used level of 0.6 mW. If very low power levels are used, no 
such effects are observed. For power levels of ~ 100 till 600 KW and with 
a constant focusing field, this anomalous attenuation is not a monotone 


function of the input power. 
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4) At the normal R.F. power level and with 4, on, a rather high amount 
of X-rays is observed over the whole lenght of the accelerator tube. pa 


5) After the accelerator has been in operation for many hours, the discs 
of the cavities show several concentric rings, differently, coloured. The dia- 
meters of these rings are different from cavity to cavity, but on all discs, the 
rings cease to appear at a few millimeter from the edge near the outside wall 
of the tube. The rings are most pronounced near the R.F. power input. It 
must be noted that all the cavities have different dimensions and are made 
to have a 7/2 mode travelling wave with a constant effective longitudinal 
electric field of 20 kV/cm. 


6) The accelerator operates best at a frequency slightly different from 
the frequency for which it was designed and matched. 


All these phenomena were reproducible as far as the effects themselves were 
concerned. Doing some rough measurements, however, lead nearly each time 
to different numerical results. 

Phenomena, similar to these described here, were also observed in some 
accelerators constructed in England (5). It was pointed out by MULLETT, 
CLAY and HADDEN (3) and by HADDEN (*) that the mentioned troubles were 
worse for dielectric loaded linear accelerators than for all-metal machines, 
and that they were less serious if mercury diffusion pumps were used instead 
of oil pumps. In the latter case, the kind of oil influenced the effect. In this 
regard it should be mentioned that the linear accelerator in Ghent has an 
all-metal corrugated wave-guide, that we had silicon oil in the diffusion pumps 
and that no cold traps were used. 

MULLETT et al. (3) explain the phenomena observed in accelerators on the 
basis of the multipactor effect. It is an effect associated with surfaces whose 
secondary emission coefficient 7 is greater than one, under the conditions 
that the effect occurs. The multipactor effect was first described by FARNS- 
WORTH (%7), was further developed by several authors (*-1°) and especially by 
KREBS and MEERBACH (1415). None of these papers, however, contains any- 
thing with relation to linear accelerators. They all consider mainly the case 


_B. R. SHERSBY-HaRrvie and L. B. MuLLETT : Proc. Phys. Soc., B 62, 270 (1949). 
W. Mircer and G. Saxon: Nature, 172, 1163 (1953). 

. B. MuLLetT, R. E. Cray and R. J. B. Happen: AERE GR/R 1076 (1953). 
J. B. Happen: AERE G/R 1161 (1953). 

B. WaLKER and E. L. Lewis: Nature, 181, 38 (1958). 

T. Farnswortu: Journ. Franklin Inst., 2, 411 (1934). 

T. FARNSWORTH: Zeits. Phys., 129, 491 (1951). 

W. B. Girt and A. von ENGEL: Proc. Roy. Soc., A 197, 107 (1949). 

W. B. GiLL and A. von ENGEL: Proc. Roy. Soc., A 192, 446 (1948). 

. J. Haron and H. S. WirLiams: Journ. Appl. Phys., 25, 417 (1954). 

. Lax, W. P. Aris and S. C. Brown: Journ. Appl. Phys., 21, 1297 (1950). 
. G. ABRAHAM: Stanford Microwave Laboratory Thesis (1950). 

. KosseL and K. Kreps: Zeits. Phys., 139, 189 (1954). 

. Kreps: Zeits. Angew. Phys., 2, 400 (1950). 

. Kress and F. MeERBACH: Ann. Phys., 15, 189 (1955). 

. Kress and F. MeeRBACH: Ann. Phys., 18, 146 (1956). 
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of two parallel flat plates with a R.F. electric field, normal to the plates. 
A simple representation can be given as follows. An electron comes out of 
one of the plates (e.g. by field-emission) at a time ¢ and is accelerated towards 
the other plate so that it arrives there at a time t+(2n-+1)/2f, f being the 
frequency of the R.F. field and n being an integer. The secondary electrons 
find the same situation at t+(2n-+1)/2f, as the first electron at time t, (if there 
is no time-delay between the arrival of an electron and the emission of the 
secondaries). If 7<1, there will be a rapid build-up. | 

The phenomenon of the coloured rings must be due to thin oil-films on 
the iris-walls. The multipactor effect itself produces cracking products, which 
can also give rise to a large increase of 7. The mechanism of the phenomenon 
has been studied and described by AITKIN (1). 

Since no higher power source was available, we have not been able to feed 
more than 0.6 mW in the linear accelerator. We should have liked to do so, 
in view of an experiment reported by MILLER and Saxon (?). They measured 
in a linear accelerator an attenuation of 2.5 db for all power levels, in the 
absence of the magnetic focusing field. In the presence of such fields, the 
attenuation was still found to be of the order of 2.5 db for power fluxes of 
a few mW and several MW. There was however a variation of attenuation 
with power in an intermediate range. They noted a maximum of 12 db for 
an input power of 250 kW. If the power level is increased, the attenuation 
falls, approaching 2.5 db for an input of 3 MW. Since, in our mind, no satis- 
factory explanation has been given for it, we want to turn our attenuation to 
this phenomenon of the dependence of attenuation on power level. 

Let us consider the condenser type electric field E sin(w@t+q), in which 
“© = 27f, p=phase angle at which the electron is emitted. Furthermore m 
and e are the mass and charge of the electron and d is the distance between 
the plates. In order to have multipactoring, the condition. 


| by 


e : 
d i: [En + 1)x cos y + 2 sin ¢], 


S 


must be satisfied. The requirement of phase stability or bunching, (then the 
multipactoring can occur in a stable way), leads to 


0<tgy <2/n(n +1), 


so that for each n, E has to lie between two extreme values: Émin and Emax: Es, 
corresponds tO @min = arctg [2/x(2n +1)] and KE, corresponds to y = 0. Table I 
gives these values for the case À = c/f — 10.00 cm, d = 2.50 cm and n = 0.1 ..., 5. 
It gives also @min and the extreme values of the energy e with which the elec- 
trons can arrive at the opposite plate. We have 


ere ron, |, 
+ [2n+1)x cos y + 28ing]? 


(1?) D. AITKIN: To be published in Proc. Inst. Electr. Eng., 105 B (1958). 
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TABLE. L. 
n P min Ein (kV/cm) Bo (kV/em) | Emin (keV) Emax (keV) 
0 32° 19! 134.3 159.2 128.1 253.1 
1 11° 59’ 51.90 53.05 25.65 28:13 
2 (Oy 31.575 31.835 | 9.81 10.13 
3 6° 127 22.643 22.735 5.094 5.162 
4 4° 03' 17.640 17.685 3.094 3.125 
5 319% | 14.445 14.468 | 2.075 2.094 | 


The regions in which the multipactoring is possible in a stable way, are 
rather narrow. For increasing n, these regions become smaller and smaller 
and the fields needed decrease. e decreases also and becomes so small that 
the multipactoring will be impossible because of the fact that the secondary 
emission coefficient becomes less than one, for low enough values of e. It is 
clear that in reality the limit values given in Table I are not very sharp. 
Among other things, this is due to the fact that the initial velocities v, of 
the electrons have a statistical distribution. (In the derivation of the above 
formulae, v was supposed to be Zero). 

The electric field E, parallel with the axis of the accelerator is, for a <r <b, 
given by 
gy Jolkr) ¥ (kb) — ¥o(kr)J (kb) 


E.) = BaF (ka) Y (kb) — Yo(ka)J (kb) ’ 


in which J, and Y, are Bessel functions, k = 2n/A and H, = (E.)--a The meaning 
of the other symbols is clear from Fig. 1. This function Æ, has been drawn 
on Fig. 2 for À = 10.00 em, a = 1.30 em and b=4.00 cm. E, is approximately 

a decreasing linear function of r. For a given value 


ir of n the multipactor effect can occur only in a ring 

whose radii r, and r, are defined by EH, = Ex and 

E, = En. It is easily seen that such a ring is very 
I Ï ] | | ! 2 narrow, especially for n #0. 

Sr a RE ani The fact that the anomalous attenuation is mea- 

d 126 sured only in the presence of the focusing field can 


be explained satisfactorily. In the corrugations, there 

is not only an electric field E,, but there is also a 

Fig. 1. component Ho, causing a radial force on the elec- 
trons, which pushes the electrons away from the 

regions where the multipactoring is possible in a stable 

way. The axial magnetic focusing field tends to prevent electrons from moving 
out of these regions. È 
If the power flux through the accelerator is increased, the values of E, 
and E, are unchanged. The interval of r(r,, 7.) where the effect can occur, 
moves however toward the outside and falls off. This is illustrated by Fig. 2. 


Let us suppose that the curves of this figure are straight lines, given by 
the equation i 


E=f0—7), (B> 0): 
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Then, the area S of the ring where the multipactor effect is possible, is 
0 (b — a) <E, 


_Pa = ely eae E, <Bp(b—a) <E, 
(1) S= B 


EG] sen 


It is obvious that, for any particular n-value, 8 is directly proportional to 
the energy loss. Fig. 3 gives S(B) for a= 1.30 cm, b — 4.00 cm, À — 10 em, 
d — 2.5 cm and n=1. S increases 
fast with , in the interval [ E,/(b— a), 
E,/(b— a)]. This interval is very 10 
small. For f > E./(b — a), S is first 
still increasing but very much slower, 
reaches a maximum for f=(£,+£E.)/b, 
and then falls off, going to 0 for 
B — co. 

For very low f, only high values 
of n could give rise to multipactoring, 
no high fields being required in that 
case. The effect will not occur, how- 
ever, because of the fact that € is 
too low to have n>1. If fis ino £)-------4-----,---- 
creased, lower n values will enter the £,,------- 
picture, for which also e is big enough 
to have 7 > 1. Increasing the power 
flux still more, will bring in more 
rings, the area of which is each time 


larger than the area of the one be- oe 
fore. The intervals between the pow- 0 NOTE TES D 
er needed for having successive Fig. 2. 


rings increase very rapidly. Even 

if the power were raised so much 

that the very low n-values could be reached, no more multipactor rings would 
come in because e would be so high that 7 would be less than one. 


i8+ units S(8) for a=1.30cm,A=1000cm, 


b=400cm,d=250cm, 
n=l 
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The most part of the energy loss, caused by multipactoring, must be due 
to the inner ring, corresponding to the lowest n-value reached in any par- 
ticular accelerator and which can give rise to the effect. This is because of 
the larger area involved and also because of the larger value of €. 

Let us now consider the results given by MILLER and SAXON. For 250 kW 
input 10 db of attenuation is due to the multipactor effect. By increasing 
the power till 3 MW, the electric field will become 4/12 times larger. From 
Table I one can see that no such large intervals exist between two successive 
multipactoring E-regions. However, the next ring corresponding to a lower 
n-value, possible as far as the needed electric field is concerned, will not give 
rise to multipactoring because of the high e value (y <1). Taking into account 
only the ring corresponding to the lowest n, that gives rise to multipactoring, 
and supposing that there is only one cavity, in which all losses occur, this 
attenuation occurs for B = (E, +E.)/b. One finds then P, =10P, (P, = input 
power, P, = output power), and S —2b?(H, — E,)/(E:+E;). For 3 MW one 
has P’=12P,, B'=V126 and S'= S(2\/12 —1)/12. From this one finds easily 
that the attenuation is now only 0.157 db. 

One comes somewhat nearer to reality if an accelerator with a purely 
periodic structure, without other losses than those due to multipactoring, is 
considered. Because of these losses each cavity has a different value for f. 
One can write 6 = A(z). Again we consider only one n-value. Using 


dP ~ pdp=— C'S dz, (C'= constant), 


and the third expression for S, one finds 


E,—EÆ, 2b 
db — 0 (Aa +) dz, (C = constant), 
so that the dependance of B on + along the accelerator is given by 
1 (60) — 6? | +E 
9 cla e 1 2 rat 
Gare otiras Ou LUS 


(E + E)? (EZ, + B,)* , B(O) — (Ai + E.)/2bX 
Nee eT HU) = Ni Rai A 


Maximum attenuation will occur approximately for 
BO =3(E + B,)/40 . 


(4 =lenght of the accelerator). Taking into account the value of a/b for most 
accelerators, this is roughly the value for which the multipactor ring is just 
in the last cavity. P,— 250 kW in the considered case and for an attenu- 
ation of 10 db one has B(0) = 34/10(#,+ F,)/4b. Using (2), this leads to 


(3) 201 (Fia) = 6.229 
E, + E, M 
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If now P,=3 MW then B/(0) = 4/12 B(0). Using (2) and (3) one finds an 
equation for A(l) leading to 


BU) = 8.126 a 


14 
: and log = = 2 log —__ = 0.022. 


Thus we have 0.22 db attenuation for 3 MW input. This agrees well with 
the experiment. 

We want to point out that the decrease in attenuation is mostly all due 
to the increase of P,. Even if the absolute power losses were constant for 
different power levels, we would find a low attenuation (0.34 db) for 3 MW 
if one supposed to have 10 db for 250 kW input. 

Finally it is interesting to remark that the increase of the input power 
will also increase the radial forces which try to push the electrons away from 
the stable region, while E remains constant in the considered ring. This will 
decrease the losses due to the multipactor effect even more. 


KOK Ok 


We wish to thank the Interuniversity Institute for Nuclear Sciences. 
(Belgium) that sponsored this work. 


RIASSUNTO (*) 


Si riferiscono fenomeni dovuti all’effetto multipactor e osservati nell’acceleratore: 
lineare dell’università di Gent. Si propone una spiegazione per la variazione dell’atte- 
nuazione anomala dovuta all’effetto multipactor in funzione dell’energia a r.f. con la 
quale si alimenta l’acceleratore. 


(*) Traduzione a cura della Redazione. 
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esponsabilita scientifica degli scritti inseriti in questa rubrica è completamente lasciata 


(La r 
dalla Direzione del periodico ai singoli autori) 


Self-Energy oî the Longitudinal Neutrino. 


W. W. WADA 
Nucleonies Division, U.S. Naval Research Laboratory - Washington, D.C. 


(ricevuto il 18 Settembre 1958) 


As is well-known the longitudinal neutrino theory is based on the assumption 
that the physical mass of the neutrino and the anti-neutrino with the inclusion 
of all interactions is strictly zero (1). The question then arises, in comparison to 
the photon self-energy, if it is possible to demonstrate that the self-energy of the 
longitudinal neutrino (and the anti-neutrino) vanishes by virtue of the so-called 
neutrino gauge. We wish to show by direct calculation that the self-energy in 
question vanishes in the lowest order perturbation calculation, in which the inter- 
actions involving neutrinos are taken to be those containing four Fermions repre- 


sented by the Puppi triangle (?). 
First we take the B-decay interactions. The interaction Lagrangian density 


. may be written (1) 


5 
T(x) a 210 (wo (x) O; Wn (a)) (We (x)O; (1 En) Vs) Wy (2)) SE 
+ OF (#,(2)0; ws (@)) (Py (@) (1 + 75) Oi pe (x))] , 


where 0,;= 1, O.= yx» O3= Gn, Or = Vas» and O;= 75. 
The S-matrix corresponding to the lowest order self-energy diagram (v= p +n+e-) 
is then found to be 


S = (2x)-*5(p'—p) > > O: ota ky) (d* ks) 


‘(U,(1 + y;)O; (yp — vhs + yka + im,)O; (1 — ys) Uy: 
trace [(yky + imy)O;(yk, + im,)O;|-[{(p — ky + ka)? + mg} (Ki + mp) (kes + mn) |, 


where m,, m,, and m, are the masses of the electron, proton, and neutron. p and 
p' are the energy-momentum 4-vectors of the neutrino before and after the self- 
energy interaction, and U, is the Dirac spinor for the neutrino. 


(2) L. LANDAU: Nucl. Phys., 8, 127 (1957); A. SALAM: Nuovo Cimento, 5, 299 (1957); T. D. LEE 
and C. N. YANG: Phys. Rev., 105, 1671 (1957). 

@) G. Pupri: Nuovo Cimento, 5, 505 (1948); also for the general discussion and other refer- 
ences, M. GELL-MANN and A. H. ROSENFELD: Ann. Rev. Nucl. Sci., 7, 407 (1957). 
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; For the purpose of analysing the above expression we will first transform the 
integrand in following steps. First making use of the identities 


1 
(ab) = [tas ui), 
0 


(ab = 2 | dr (1-2) [ar + b(1— a) 
ò 


we may write (g = k,—k,) 


[((p — 9)? + me) (ki + me) (kz + mi) = 


1 


= 2fea --2)d+[(a—2) de magma x2(1 — z)m? + 
0 0 


+ {(ke + 92)? + mal — 2) + m5 2}(1—-a)|- 


Then introducing variables a and 6 defined by k,+qz=a and q — px=b, we obtain 


1 1 
S = (22)-16(p’— p EE 0,05fannasfe mas fam. 
ae 0 0 


(D, + y5)0; {(1 — x)yp + im, — yb} 0; (1 — y;)U,)- 
‘trace [ (ya) O; (ya)O; + {yb + (pb) (1 —2)a+ im,} O; {— 2(yb) — ex(yp)+ ima} 0;]: 
-[2(1 — 2) b? + 22(1—2)m? + (1 — x) (a? + m3 (1 — 2) + emd)|3 


Let us first take i=j=1, or 0;=0;=1. Then, one may set 
(U,(1 + 5) {(1 — x)yp + im} (1 =" 5)U,) 0, 


immediately, since (yp) U,=0 because of the zero rest mass of the neutrino. If the 
neutrino were not longitudinal, then its wave function would not be in the form 
of (1—y;) U,, and therefore the above expression would not vanish even for the 
vanishing neutrino rest mass. For the longitudinal neutrino we then find that the 
integrand is proportional to (U,(1+y5) (vb) (pb1—y:)U,)., which, when averaged 
over directions of the 4-vector d, becomes a(Q(1 ++y;)b2(yp)(1 — y; )U,): Thus, we 
find that the above particular term (i=j=1) is of the form of the product of the 
rest mass of the neutrino and acertain divergent integral, vhich may be regarded 
as zero. It can be easily shown that the remaining terms, i.e. i=j71 and the cross 
terms, each vanish in similar manner. This situation persists with regard to other 
interactions in the Puppi triangle, namely vZp+n+u and væu +et+v. Thus, 
the self-energy of the longitudinal neutrino arising from the interactions involving 
four fermions vanishes in the lowest order perturbation calculation by virtue of the 


neutrino gauge. 


47 - Il Nuovo Cimento. 
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Some Remarks on the Quantum Mechanies of a Lattice Gas. 


J. CZERWONKO = 


Physics Department, Bolestaw Bierut University - Wroctaw 


(ricevuto il 27 Settembre 1958) 


It is shown that there exist two different Hamiltonians for a lattice gas which 
correspond to the same Hamiltonian of a system of N particles in the theory of 
second quantization. Calculations based on these two Hamiltonians lead obviously 
to different results. It seems a priori impossible to find theoretical arguments in 
favour of one of them. 

Following MatsuBara, Matsupa (1) and Morita (?) we consider the lattice 
gas in which N molecules are on L lattice sites (7,...7,). We use here arbitrary 
lattices satisfying the following conditions: 


a) there exist vectors 7,;, à = (— 2/2, ..., — 1, 1, ..., 2/2) where 2 is the co-ordi- 
nation number of the lattice; such that: 


nyt Th = Tra; and |n;|= d = min|r,— rm). 
sm 
b) n= — 7-3 n £ N; if à #j. We define an operator of the ratio of differences: 
1 
(1) Ani (= (Ts Gera 1D \p UE SACS. Tr 25) = CAES Ve Ni 1x) . 


We have (because of b)) 
(2) Tia Ti= 1. 


From (1) we have: 


1 
(3) Aia = FD (Ti. lai 27,3) . 
6 


(1) T. MATSUBARA and H. MATSUDA: Progr. Theor. Phys., 16, 596 (1956) and 17, 19 (1957). 


(3) T. Mortra: Progr. Theor. Phys., 18, 462 (1957). 
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We take into consideration only state functions y sara on the set of. lattice 
points rx. It may be easily seen that by taking n;=d(d;1. 0:35 Og) for i=1;.2, 3 and 


n-i= —; we get vectors possessing the properties a) and b) for a simple cubic 
lattice. By virtue of (3): 


1 
dipl Se ae [ y(... Ly + 2d...) — QWyp(... M+ a...) + yl... My cette 


If we compare the Morita definition of the operator A2/Ax% with this result, we 
see that: 


A? 


> 2 
Ara a Late - 
AT 


Taking into account the correspondence with the continuous case we see that 
the kinetic energy ought to be proportional to the sum of squares of displacement 
operators changing the argument of the state function by 7,. Thus for the operator 
A/Ax, such that (4/Ax,)2? = A2/Axz we have: 


where 


Tt sy... dpr) == Vidi) 


As we see the operator /1//4x, operating on the state functions creates functions 
which are not state functions since they depend on points which are not lattice 
points. Thus it is seen that the Morita operator of kinetic energy cannot be treated 
as the sum of squares of operators. in this space. 

We shall now show that to the two equivalent hamiltonians 


2. h,2 


he. 2 a a h 
pre VC, ANG, op 


[vs V?y(dr), 


correspond two essentially different hamiltonians for the lattice gas. 
2/2 
We replace here the expression Vy* Vy by = Ani v*4x;w, the expression w* V?y 
2/2 


by > p*/A%w, and the integral by the sum over the lattice points. We obtain: 


h2 212 
— | Vy* Vy(dr) +  — An: P* (rs) Ara P(r) » 
(4) =f y* Vy(dr) 3 > ni P(e) Ax (1%) 
h2 h DI ape 5 
el *V2uy(dr) = W* (1) Ari P(r) 
(5) =v V2 y(dr) > AZZ (1%) Ax. (T7) 


The expressions (4) and (5) go over into hamiltonians for the continuous position 
spectrum when d tends to zero, only for simple cubic lattices. This is caused by 
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a linear dependence of vectors 7;; j=1,...,2/2. On the other hand in expressions (4) 
and (5) we must take into account summations over all displacements 7;, 1°", 2/2; 
Let us consider the expression (4). We take after Morita: 


L 
wre) = >. dry Wi» 
j=1 


where the operators a; and a* fulfill the commutations relations given by Matsubara 


and Matsuda. We obtain: 


L L 
LE ones emi  — dry MI LX Brgere tne Om — Org re Am] = 


2md* pal Dent dal m=1 


h2 L ‘212 à 4 he k 
NES 2 > (aan; — Ar) (Ae 4m; — ay) = md? 2 — a;)(a; — G;) . 


Summation over <i, j) denotes summation over all the nearest pairs of points. For 
simple cubic lattices the hamiltonian (4’) equals that given by Matsubara and Matsuda. 
Let us consider the expression (5). We obtain: 


h2 L  af2 L L 
, 3 * a —- 9 — 
Mg ) P 2 12 i? a Or s.r} [> (dm Or mary +216 24m Orn etn: + Um One) 
MA = TL tal g=1 m=1 
hg £ h2 h2 L 212 
* * * 
= cs Od Gi apy Oy, On st, * 
amd? 21 md? 2 OM roe ae + 


where 7,, denotes the index of the vector 7+ 2n;. As we see the expression obtained 
in this way is quite different from (4’). It is interesting that the expression contains 
interactions with lattice points with the distance 2d as well, and it does not contain 
the interaction terms describing interactions with lattice points at distances smaller 
than 2d and larger than d. i 

There appears the question why these two lattice point hamiltonians are not 
equivalent. This is connected with the fact, that in difference calculus: 


A(ab) = (Aa)b + (Ta)(Ab). 


Hence: 
212 L 
(6) O= D Dd Aus (w* (7) Ax, v(1)) = 
i=1k=1 
a2 L 2/2 L 
= > > Age WS (re) Ari Pr) + > > (Tra y* (rx)) Ari (Tr) - 
is ped SZ i 
Hence: 
h2 2/2. ZX h2 2i2 L 
(7) Da > bi Ana D (Tr) Ana Pre) = — ome > > Tri (W* (Te) Tri Ati p(x) « 
421 bal Mi 1 K=1 


The last equality is the consequence of (6), (2), and of the connexion T(ab)=TaTb. 
The expression (7) up to the operators a; and aj at the band-points of the lattice 
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is identical with: 


h2 zl2 L 


o a > (y* (1a) Tes Aine y(Tr)) - 


2m 2121 


This expression for the simple cubic lattice equals the hamiltonian given by Morita. 

By elementary computations we obtain the result that hamiltonians (4’) and (5’) 
commute with the total number of particles. Thus we have no reasons for pre- 
ferring one of them. | 

It may be noted that the transition from the quantum mechanics with con- 
tinuous position spectrum to quantum mechanics of lattice gas is very difficult, 
because of two reasons: (i) the classical analogue does not exist and (ii) it is im- 
possible to use the canonical formalism, because of: 


[Axa» xi] ni din Tra o] 


On the other hand, as shown in this note, a procedure based on the correspondence 
with the case of continuous position spectrum is not unique. 
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Pion-Proton Scattering Dispersion Relations. 


L. BertoccHI and L. LENDINARA 


Istituto di Fisica dell Universita - Bologna 
Istituto Nazionale di Fisica Nucleare - Sezione di Bologna 


(ricevuto il 20 Ottobre 1958) 


We have re-examined the Puppi-Stanghellini (P.S.) discrepancy in the x-p scat- 
tering dispersion relations (1). 

The main difference between our analysis and that of P.S. is the use of dif- 
ferent values of the low energy 7--p total cross sections for the evaluation of the 
integrals. The data which P.S. used for drawing the low energy total cross section 
were: (12.9-+1.7) mb at 37 MeV (?), (17.6+2.2) mb at 58 MeV (°) and the two values 
(21-8) mb at 110 MeV and (31+9) mb at 113 MeV (*), the last two with large errors; 
P.S. did not use the cross section (14.7-+1.8) mb at 65 MeV (5), obtained by inte- 
gration, which was in disagreement with the other data. 

Since the P.S. analysis, other experiments have been finished. These new total 
cross section data yield (7.8+0.7) mb at 41.5 MeV (6) and (21.1-+0.6) mb at 98 MeV (7). 
Taking these new data into account, the best fit to the ensemble of total eross sections 
below the resonance is substantially lower than the cross sections used by P.S. 
It follows that also the 65 MeV point is in good agreement with the other data, 
and the only point remaining in disagreement is the 58 MeV value. This change 
has a negligible effect on the D* dispersion relation, owing to the fact that the 
integral for D* involves only o- divided by w'+ ©. But the modification may be 


noticeable in the D dispersion relation, especially in the region in which the slope 
of o- is increased. 


(1) G. Puppi and A. STANGHELLINI: Nuovo Cimento, 5, 1305 (1957). 

(3) C. E. ANGELL and J. P. PERRY: Phys. Rev., 92, 835 (1953). 

(3) P. J. Isaacs, A. M. SACHS and J. STEINBERGER: Phys. Rev., 85, 803 (1952). 
(ELSA 


. ANDERSON, E. FERMI, E. A. LONG, R. MARTIN and P. E. NAGLE: Phys. Rev., 85, 
934 (1952). 
6) D. BODANSKY, A. M. SACH3 and J. STEINBERGER: Phys. Rev., 93, 1367 (1954). 


(9) S. BARNES, B. Rose, G. Gracomenni, J. RING and K. MryAKE: University of Rochester 
Report NYO-2170. 


() J. R. Hour: private communication to prof. PUPPI. 
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The result of the calculation of D with these new data is shown in Fig. 1, 
in which the two points obtained from the angular distributions at 41.5 MeV and 
98 MeV have been added. In Fig. 1 the curve calculated independently by SALZMAN (8) 
using the Anderson-Piccioni cross-sections (*!°) is also shown. The good agreement 
between the two curves is due to the 
fact that the new total cross sections ob. oP + sy ob-0d) 
used by us are consistent with the ve ia 
Anderson phase. shifts. 


----- Salzman 
Bologna 


id 
w (pi) 2) 
o 15 20 D die 
-0 È 
ga (MeV) 
-02 Dee Le è 
0 100 200 300 


Fig. 1. — D?" versus the total energy in Fig. 2.—(D°_+D}.)/2r,+(1/2Mo):(D°_ — D? )/2r0 
units of 4 with the subtraction at o = w. versus the total energy. 


Comparing with the P.S. results, an improvement in the region below and around 
the resonance may be noticed. No appreciable improvement is noticed in the region 
above the resonance. But one must take into account that all points have been 
obtained with an analysis of the angular distribution made only with S and P waves. 
As Hamilton and Chiu (!!) have shown, the points around 300 MeV, analyzed with 
S, P and D waves, have large errors in their real part of the scattering amplitude, 
and in addition give values rather different from those obtained an alyzing only with 
S and P waves. Therefore these points cannot discriminate between the various 
D> curves above the resonance. Thus, the only points which remain in disagreement 
with a value of f? around 0.08 are the two Ashkin points (12) at 150 and 170 MeV. 

In spite of the fact that Hamilton has suggested that one can resolve this dis- 
agreement by modifying the total cross section used in the optical theorem at these 
two energies, we have mantained the P.S. values of D° because at the moment 
no new experimental data are available in this region. Fig. 2, in which is plotted 


(8) G. SALZMAN: private communication to prof. PUPPI. 

(®) H. L. ANDERSON and N. METROPOLIS: Proc. of the Sixth Annual Rochester Conference (1956). 
(2) R. Coot, 0. Piccioni and D. CLARK: Phys. Rev., 103, 1082 (1956). 

(11) J. HAMILTON and H. Y. Caru: Phys. Rev. Lett., 1, 4 (1958). 

(12) J. ASHKIN, J. P. BLASER, F. FEINER and M. O. STERN: Phys. Rev., 101, 1149 (1956). 

(3) J. HAMILTON: Phys. Rev., 110, 1134 (1958). 
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the relationship 
DID are Ike 


275 2M 215 


in which the f? term disappears, shows in a clear way the discrepancy between the 


150 and 170 MeV points and the total ensemble of experimental data. 

In this analysis, as well as in that of P.S.. the subtraction is made at the total 
energy @= u. This way of handling the dispersion relations may introduce at 
least two sources of trouble: 


1) As CHEW et al. (14) have pointed out, the P.S. analysis was not completely 
from free the hypothesis of charge independence, because they used Orear’s scattering 
lengths, which were obtained from elastic scattering experiments as well as from 
charge exchange experiments. 

2) Because of the mass differences between charged and neutral pions, and 
between protons and neutrons, the total o- cross section is different from zero at 
©= n, so that the integrals appearing in the dispersion relations are meaningless, 
except if one makes particular assumptions about the behaviour of the imaginary 
part of the scattering amplitude in the region between wo and x, assumptions which 
we prefer to avoid. 


For these reasons, it is more convenient to make the subtraction at a total 
energy different from 4, at which the real parts of the scattering amplitudes are 


TR b 
0e g $ D w )=0.104 * 0.014 


OPC w 6):0102: 00% 
w6=41.5 MeV 


-O45- 


07 
HA ed 


Fig. 3. — D, jr) versus the total energy with the subtraction at © = 41.5 MeV. 


(14) G. F. Carew and H. P. Noyes: Phys. Rev., 109, 566 (1958). 
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1 


experimentally known. Such an analysis has already been made by Cuiu (5), who has 
chosen the kinetic energy 150 MeV as the subtraction energy. However, this par- 
ticular energy is not very useful for the study of the region around the resonance, 
since all the curves are forced to pass through the subtraction point, independently 
of the value of f?. 

So we have chosen as subtraction point the Rochester 41.5 MeV point (f), and 
we have calculated the two relations: 


foe} 


Dh < 2}? (1— 1/42) 1 | n do' | ot g- + 


2 
ae he a à | et 
lire n Gi QE o + 1/2M (@ì — 1/4M?) la 4n°r9,) 0° — wo o Fo: ow +o 


The results are shown in Figs. 3 and 4. 
For the D* relations the situation is the same as P.S. found, although the curves 
are less sensitive to the value of f?. The D relation, in comparison with the P.S. 


DÒ (w)=0.104* 0.014 


p> 
14 02( w,)=0102 + 0.034 


0.2- 


-01 


-02 


Fig. 4. — D! Ir versus the total energy with the subtraction at w = 41.5 MeV. 


(25) H. Y. Curu: Phys. Rev., 110, 1140 (1958). 
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curves, shows some improvements, there being essential agreement with our result 
obtained subtracting at © = y. Unfortunately, in this way, the continuous curves 
have errors which arise from the real parts at fixed energy which are much larger 
than those errors arising from the scattering lengths in the w = y subtraction case. 
Thus one cannot obtain definite conclusions from dispersion relations used in this 
manner, until more precise experiments are finished in the low energy region. 

In conclusion, from the present values of the experimental results one can say 
that the general behaviour of the data agrees with the prediction of dispersion relations 
and is compatible with a coupling constant f? around 0.08. There is a remaining 
discrepancy in the region between. 150 and 200 MeV, which cannot be solved with 
further manipulations of the dispersion relations, but only with new experimental 
results. Thus, until these new experiments are finished, one cannot worry about 
the origin of this discrepancy or look for a violation of some fundamental principle. 


* KOK 


We wish to thank prof. G. Puppi and Drs. A. MINGUZZI and A. STANGHELLINI 
for stimulating discussions. 
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LIBRI RICEVUTI E RECENSIONI 


R. J. BLIN-STOYLE — Theories of 
Nuclear Moments. Pagg. 89, fig. 12, 
Oxford University Press, 1957 


Come dice il titolo stesso, e come è 
espresso nella prefazione, scopo di questo 
volumetto è l’inquadramento teorico dei 
risultati sperimentali concernenti i mo- 
menti nucleari in base ai vari modelli 
del nucleo. 

In questo senso esso potrebbe essere 
un utile complemento ai testi di argo- 
mento affine, come ad esempio quelli di 
Ramsey e di Kopfermann, nei quali i 
momenti nucleari sono considerati pre- 
valentemente dal punto di vista della 
loro interazione con campi interni (ato- 
mici, molecolari, cristaliini), od esterni. 

Purtroppo però l'argomento viene 
condensato in un numero di pagine 
esiguo rispetto all’ampiezza delle que- 
stioni trattate, per cui c'è da dubitare 
che un lettore non specializzato possa 
trarre un effettivo vantaggio dalla sola 
lettura del libro. 

L'utilità del volumetto si può rico- 
noscere se lo si considera piuttosto come 
una guida alla consultazione della biblio- 
grafia in esso citata. 

G. BOLOGNA 


Reports on Progress in Physies, Vo- 
lume XXI (1958), a cura di A. C. 
STICKLAND; Ed. The Physical So- 
ciety, London, pagg. 383. 


La mole di questo ventunesimo vo- 
lume della raccolta di lavori monografici 
pubblicata dalla Physical Society di 
Londra è alquanto minore di quella del 
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volume precedente, segnando così una 
battuta di arresto nella serie rapidamente 
crescente di quelli pubblicati negli ultimi 
anni. Ciò non significa tuttavia che 
l'interesse della pubblicazione possa con- 
siderarsi minore, perchè essa conserva 
l’alto livello nel contenuto e l’accura- 
tezza nella presentazione che sempre 
l’hanno distinta. i 

La scelta degli argomenti, che vanno 
dalla biofisica alla fisica nucleare, dalle 
tecniche di laboratorio agli sviluppi della 
teoria, è fatta come sempre dando la 
preferenza a problemi moderni e di largo 
interesse. 

È notevole in questo volume il posto 
occupato dalla Fisica dei Solidi alla 
quale sono connesse ben tre delle nove 
monografie; in una di queste dedicata 
ai calori specifici dei metalli a bassa 
temperatura D. H. PARKINSON, dopo 
aver separatamente esaminato il con- 
tributo reticolare e quello elettronico, 
si occupa brevemente dei casi specifici 
dei superconduttori e delle terre rare, 
concludendo con una esposizione delle 
più recenti tecniche sperimentali. Uno 
studio sugli spettri di assorbimento e di 
fluorescenza degli ioni nei cristalli dovuto 
a W. A. RUNCIMAN riunisce poi i risul- 
tati più recenti ottenuti per i vari tipi 
di ioni con particolare riguardo allo 
studio dei livelli energetici. Sempre in 
questo gruppo si può poi considerare 
l’accuratarassegna critica sull'analisi delle 
strutture cristalline a mezzo della diffra- 
zione degli elettroni di J. M. CowLEY e 
A. L. G. Regs, in cui si cerca di coordi- 
nare i metodi e la nomenclatura illu- 
strando in particolare l'applicazione dei 
metodi di Fourier. 

Passando agli altri argomenti ricor- 
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diamo che M. Brrnx in un articolo di 
carattere essenzialmente teorico si oc- 
cupa dei recenti progressi relativi alla 
interpretazione della struttura del nucleo, 
con particolare riguardo agli aspetti 
« collettivi », mentre R. HuBy tratta 
della eccitazione elettromagnetica dei 
nuclei a mezzo di proiettili nucleari o di 
elettroni, con una esposizione equamente 
bilanciata della parte tr rica e di quella 
sperimentale. i 

Un lavoro sugli atomi mesonici 
(D. WEST) ne passa in rassegna le carat- 
teristiche principali insieme alla infor- 
mazione che se ne può trarre sulle pro- 
prietà dei mesoni e dei nuclei, tra loro 
interagenti; anche in questo lavoro è 
fatto il debito luogo sia alla discussione 
dei dati sperimentali che alla teoria. 

J. M. VALENTINE ed S. C. CURRAN 
si occupano poi dell'energia mediamente 
spesa per ogni coppia di ioni prodotti in 
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un gas da vari tipi di radiazione descri- 
vendo gli orientamenti più recenti di 
questi studi e fornendo interessanti ta- 
belle di risultati. 

La produzione-e la misura di altis- 
sime temperature hanno fatto impor- 
tanti progressi connessi anche con i 
modernissimi studi sulla fusione control- 
lata; W. LOCHTE-HOLTGREVEN ci offre 
qui un quadro, notevole per la sua com- 
pletezza, dei metodi che possono essere 
impiegati e dei risultati ottenuti; la 
produzione di alte temperature con 
mezzi chimici è esposta in appendice 
da R. ScHALL e F. WECKEN. 

Finalmente la fisica del sistema ner- 
voso, un tipico esempio di collabora- 
zione fra le scienze fisiche e quelle bio- 
logiche, è trattata da P. FATT, uno spe- 
cialista di biofisica dello University 
College di Londra. 

F. A. LEVI 
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